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ABSTRACT 


The  Output  Processor  Module  of  the  Department  of  Defense  Land  Fall¬ 
out  Prediction  System  is  described  and  instructions  are  given  for  its  use. 

Working  in  close  liaison  with  the  Particle  Activity  Module  (Volume  V),  the 
Output  Processor  converts  the  output  of  the  Transport  Module  into  a  vari¬ 
ety  of  displays  in  a  directly  contourable  numerical  (map)  form  by  means 
of  the  off-line  primer.  It  requires  only  two  sets  of  input  data  in  addition 
to  the  inputs  called  for  by  the  Particle  Activity  Module:  (1)  a  magnet!"  tr.pe 
containing  descriptions  of  sets  of  grounded  fallout  particles  —  an  input  from 
the  Transport  Module,  and  (2)  card  inputs  by  vhich  the  user  may  request 
any  number  of  processing  tasks  to  be  carried  out  cn  the  grounded  fallout 
particle  data.  In  each  request  any  of  sixteen  distinct  types  of  processing 
maybe  specified  leading  to  the  display  of  maps  of  any  of  the  following  quan¬ 
tities:  (1)  exposure  rate  "normalized"  to  H  +  1  hour;  (2)  exposure  rate  at 
time  H  +  T1  hours;  (3)  integrated  exposure,  K  +  Tl  to  infinity,  accounting 
for  time  of  arrival;  (4)  Integrated  exposure,  H  +  Tl  to  H  +  T2,  accounting 
for  time  of  arrival;  (5)  fallout,  mass  (per  unit  area);  (6)  fallout  mass  (per 
unit  area)  deposited  between  times  H  +  Tl  and  H  +  T2;  (7)  integrated  expo¬ 
sure,  H  +  Tl  to  H  +  T2,  assuming  all  particles  have  arrived  by  H  +  Tl 
hours;  (8)  same  as  7  but  integrated  to  infinity;  (9)  concentration  of  an  indi¬ 
vidual  mass  chain  (curies/m2);  (10)  time  of  arrival;  (11)  time  of  cessation; 

(12)  smallest  particle  deposited;  (13)  largest  particle  deposited;  (14)  mass 
deposited  by  particles  in  the  size  range  Si  to  S2;  (15)  H  +  1  hour  "normal¬ 
ized"  exposure  rate  resulting  from  particles  in  the  size  range  SI  to  S2;  s* 

and  (16)  the  number  of  cloud  (model)  subdivisions  affecting  each  map  grid 
point.  The  user  is  free  to  specify  any  limiting  coordinates  and  scale  fac¬ 
tors  for  the  map  display  that  will  be  produced  and  can  also  cause  the  result¬ 
ing  map,  or  maps,  to  be  recorded  on  magnetic  tape  for  further  processing. 
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INTRODUCTION 


i  hi. s-  volume  is  intended  to  fulfill  two  needs:  (1)  to  provide  information  to  the 
perron  who  is  interested  only  in  understanding  the  Output  Processor  or  in  using  it 
nK  !i  nnd  (2)  to  provide  a  more  detailed  explanation; of  the  programs  and  their 
functions  to  the  researcher  or  programmer  who  would j make  uunJificatious  or  addi¬ 
tion:.,  The  sections  entitled  "Program  Description,  "  ^"Illustration  of  Output  Pro¬ 
cessor  Use,  "  and  "User  Information"  are  Intended  ro  fulfill  the  first  need;  the 

\ 

section  "Program  Details,  "  the  second  need. 

PROGRAM  DESCRIPTION! 

i 

i 

The  Basic  Operation  of  the  Output  Processor  j 

'fiie  Output  Processor  of  the  DOD  T.and  Fallout  Prediction  System  is  a  very 
flexible,  highly  modular  computer  program  for  use  in  the  interpretation  of  data 

i 

representing  grounded  subdivisions  of  the  radioactive  cloud.  In  simplest  terms, 

i 

it  is  the  task  of  the  Output  Processor  to  accept  descriptions  of  grounded  cloud 
subdivisions,  make  requests  for  particle  activities  or  mass  chain  concentrations 
from  the  Particle  Activity  Moduie,  interpret  them  into1  a  two-dimensional  memory 

array  or  map  image,  and  then  the  resulting  array  in  a  form  suitable  for  view- 

l 

ing  as  a  map.  The  processes  originally  required  of  the  Output  Processor  were 

t- 

x 

the  computation  of  (1)  exposure  rate  "normalized"  to  H  +  1  hour,  (2)  exposure 

V 

rate  at  a  specified  time,  (3)  exposure  accumulated  between  two  specified  times, 

(4)  particle  mass  deposited  per  unit  area,  and  (5)  concentration  of  a  user  specified 
mass  chain.  In  certain  of  these  processes  the  time  of  particle  arrival  on  the  fall- 
out  field  was  also  to  be  accounted  for.  \ 

1  I  ' 


*  i 

Exposure  rate  patterns  "normalized"  to  H  +  1  hour  are  intended  to  show  the 

exposure  rates  that  would  exist  on  the  fallout  field  one  hour  after  detonation  if  all 
radioactive  particles  that  ever  land  on  the  fallout  field  were  located  in  their  places 
of  deposit  at  that  time.  Obviously  this  differs  from  the;  exposure  rate  pattern  pre¬ 
dicted  to  exist  one  hour  after  detonation,  since  in  actudl  exposure  rate  prediction 
we  must  account  for  the  times  at  which  particles  arrive  on  or  near  the  ground. 


i 

i 

i 

|  1 

f 
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The  following  statements  of  requirement  and  intent  describe  some  of  our  initial 
motivations  and  justifications  for  the  approach  which  was  followed  in  the  construc¬ 
tion  of  the  Output  Processor: 

1,  Groat  flexibility  in  program  uee  should  be  allowed  both  in  terms  of  the 
nature  of  computations  and  tasks  and  in  terms  of  the  degree  of  precision 
in  both  modeling  and  display. 

2.  The  Output  Processor  should  be  c: vAhle  of  handling  a  large  set  of  grounded 
particle  data.  TV  size  of  this  set  might  vastly  exceed  available  memory 
space.  Thus  an  open-ended  philosophy  was  adopted  for  the  treatment  of 
data  on  grounded  particles. 

2.  The  position  and  scale  factor  of  the  map  should  be  under  the  direct  control 
of  the  researcher.  This  gives  the  user  or  researcher  the  ability  to  produce 
maps  of  any  desired  scale  factor  for  superposition  on  other  maps  and 
enables  him  to  achieve  either  a  microscopic  or  a  macroscopic  view  of  the 
predicted  fallout  field. 

4.  The  Output  Processor  should  be  capable  of  handling  output  maps  containing 
a  larger  number  of  map  grid  polnta  than  can  be  represented  in  the  computer 
memory  at  one  time.  TK  led  us  to  an  open-endedness  in  output  map  size. 

5.  fa  computing  exposure  rates  at  arbitrarily  specified  times,  it  is  deemed 

of  great  importance  to  avoid  reliance  on  a  single  exponential  decay  constant 
(such  as  -  1 . 2)  which  is  truly  applicable  only  to  a  mixture  of  unfractionated 
fission  products  —  not  in  general  to  isolated  samples  of  fallout  such  as  those 
that  appear  locally  In  fallcuc  fields.  Therefore,  the  Output  Processor 
Module  should  be  built  to  work  in  close  liaison  with  the  Particle  Activity 
Module  so  that  particle  activities  can  be  computed  directly  from  the  primary 
mass  chain  data  for  (deposited)  particles  at  the  particular  time  or  times 
specified  in  each  output  request.  Furthermore,  by  means  of  this  approach 
it  should  he  possible  to  compute  and  display  concentrations  of  any  (user) 
specified  mass  chain. 


6,  With  regard  to  display  of  the  fallout  map  data  produced  by  the  Output  Pro¬ 
cessor  we  are  faced  with  somewhat  conflicting  requirements:  (1)  we  desire 
a  numerical  display  of  the  data  rather  than  some  sort  of  purely  pictorial 

or  graphical  display  because  of  the  intended  research  application  of  the 
system,  whereas  (2)  an  automated  pictorial  or  graphical  display  relieves 
the  researcher  of  the  time  consuming  and  tedious  task  of  transferring  nu¬ 
merical  data  to  a  grid  and  hand  contouring  isoexposure  lines.  The  display 
actually  provided  is  a  compromise  between  these  two  extremes.  A  numer¬ 
ical  display  is  provided;  however,  it  is  in  a  format  that  allows 
strips  of  the  printed  computer  output  to  be  attached  side-by-side  so  that 
the  entire  fallout  prediction  area  is  included  on  the  assembled  paper.  Thus, 
the  printed  numbers  represent  exposure  rate  (or  some  other  output)  pre¬ 
dicted  for  each  of  the  points  of  a  regular  grid  that  can  be  arranged  to  be 
spatially  undistorted,  and  the  resulting  map  can  be  easily  contoured  di¬ 
rectly  on  the  printer  output  paper.  The  major  disadvantage  of  this  type  of 
display  is  the  large  size  of  the  maps  produced. 

7.  In  general,  it  was  desired  that  the  Output  Processor  be  simple  to  use  and 
be  reasonab'y  foolproof  and  automatic  with  respect  to  its  internal  oper¬ 
ations.  Since  the  sizes  of  input  and  output  data  sets  r/ere  assumed  to  be 
widely  varying,  this  led  us  to  a  certain  amount  of  essentially  "dimension 
free"  programming  with  the  objective  of  making  it  unnecessary  in  most 
situations  for  the  user  to  explicitly  modify  memory  allocations  (dimension 
statements)  and  recompile  programs  in  order  to  change  the  program's 
scale  of  operation. 

Flexibility  of  the  Output  Processor 

For  a  research  system  —  one  which  is  capable  of  aiding  the  researcher  in  his 
many  and  varied  tasks  —  no  single  approach  to  flexibility  is  sufficient  and,  conse¬ 
quently,  we  ha,re  designed  the  Output  Processor  with  three  modes  of  flexibility  In 
mind:  (1.)  program  modularity,  (2)  parameter  controlled  options;  and  (3),  as  a  mid¬ 
ground  between  these  two,  code  insertion  points. 
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First,  functional  subroutines  have  been  designed  to  operate  wherever  functions 
could  be  clearly  seen.  For  example,  within  the  Output  Processor  there  are  a  sub¬ 
routine  (CALC)  with  the  primary  function  of  interpreting  ground  cloud  subdivisions 
into  a  two-dimensional  array,  and  a  separate  subroutine  (MAP)  with  the  function  of 
composing  and  displaying  the  computational  results.  These  functional  subroutines 
may  be  relatively  easily  replaced  by  other  subroutines  having  similar  purposes. 

As  an  example  of  the  second  mode  of  flexibility,  in  CALC  all  currently  required 
computational  tasks  have  been  accounted  for  and  these  computational  alternatives 
are  treated  as  parameter-controlled  options  within  the  program.  In  order  to  select 
one  of  the  available  computational  options,  the  user  need  only  punch  on  an  input 
card  the  appropriate  numerical  value  for  an  input  parameter,  in  this  case  the  param¬ 
eter  NREQ.  Furthermore,  provision  has  been  made  in  the  CALC  option  selection 
procedure  (as  well  as  in  those  of  other  programs)  for  the  future  inclusion  of  other 
computational  option  codes  with  little  or  no  modification  to  the  control  programming 
within  the  subroutine.  Such  insertion  points  (noted  in  the  flow  charts  and  card  list¬ 
ings)  represent  the  third  mode  of  program  flexibility. 

Inputs  to  the  Output  Processor 

The  primary  input  to  the  Output  Processor  is  the  magnetic  tape  of  grounded 
particle  descriptions  which  is  produced  by  the  Transport  Program.  For  each  in¬ 
cluded  central  particle  (representing  a  cloud  subdivision)  this  tape  contains  the  two 
horizontal  coordinates  of  its  impact  point,  its  time  of  impact,  the  central  particle 
oize,  and  a  mass  per  h  'rizontal  unit  area  covered  by  the  cloud  subdivision.  In 
addition,  this  tape  contains  a  tabulation  of  particle  properties  (mass  and  sur- 
face/volome  ratio)  as  a  function  of  particle  size  range. 

In  addition  to  this  primary  input  the  user  must  communicate  to  the  program 
his  wishes  regarding  the  kind  of  output  computation  and  its  form  of  presentation. 

He  must  also  provide  run  identification  information  and  information  on  certain  im¬ 
portant  computer  features.  The  run  identifier  is  an  arbitrary  72  character  state¬ 
ment  which  the  user  can  set  to  identify  and  associate  outputs  and  inputs.  The 
other  inputs  are  needed  to  allow  the  program  to  adapt  to  some  degree  to  different 
computer  environments. 
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The  following  algebraic  sentence  summarizes  the  operation  of  the  Output 
Processor: 

Particle  impact  data 

+  available  tapes  and  printer  characteristics 
+  map  characteristics 
+  computation  option  specification 
+  display  option  specification 

=>  desired  presentation . 

By  available  tapes  we  mean  simply  the  identifiers  of  magnetic  tape  units  that  are 
available  for  temporary  use  by  the  Output  Processor.  The  printer  characteristics 
are  the  character  spacing  constants  for  the  off-line  printer.  The  map  character¬ 
istics  are,  at  this  time,  the  description  of  the  geographical  limits  and  data  point 
density  of  the  map  which  is  to  be  produced  as  output.  By  computation  option  we 
mean  the  choice  of  one  of  the  many  alternative  output  quantities  to  be  computed  and 
displayed.  By  display  option  we  mean  the  choice  of  a  particular  printed  map  for¬ 
mat;  perhaps  in  the  future  it  could  also  specify  a  format  for  another  kind  of  output 
device . 

Available  Options  for  Computation  and  Display 

The  following  is  a  listing  and  brief  discussion  of  the  major  options  for  compu¬ 
tation  and  display  which  exist  in  the  Output  Processor.  An  exhaustive  list  of  all 
currently  available  options  is  prodded  in  the  section  entitled  "User  Information.  " 

1.  Printed  descriptions  of  Impacted  particles 

Under  this  option  the  content  of  the  grounded  particles  tape  may  be  printed 
in  a  form  analogous  to  that  in  which  it  exists  on  the  transport  tape  (IPOUT). 
This  option  is  valuable  in  checking  the  execution  of  experimental  transport 
codes,  and  it  is  also  useful  in  providing  a  hard  and  readable  copy  of  the 
result  of  transport  production  runs. 
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Computational  options 

The  descriptions  below  apply  to  each  point  of  the  map  grid. 

a.  Count  of  grounded  cloud  subdivisions.  This  optional  computation  was 
of  primary  value  in  debugging  the  Output  Processor  but  may  also  be 
of  considerable  value  to  the  researcher  in  assessing  the  statistical 
validity  of  a  computed  map  quantity  at  any  particular  point  on  the  map. 

b.  Exposure  rate  "normalized"  to  time  H  +  1  hour.  This  is  the  recog¬ 
nized  standard  mathematical  construct  for  the  comparison  of  fallout 
patterns.  It  should  be  noted  that  differences  may  exist  between 
DELFIC  H  +  l  hour  normalizations  and  those  resulting  directly,  or 
indirectly,  from  backward  extrapolations  of  field  data  —  in  backward 
extrapolations  one  decay  constant  is  usually  used,  whereas  DELFIC 
provides  a  more  rigorous  modeling  of  radioactive  decay. 

c.  Exposure  rate  at  t'me  H  +  Tl.  This  is  actually  the  exposure  rate  at 
the  specified  time  taking  into  account  the  impact  times  of  all  cloud 
subdivisions. 

d.  Exposure  accumulated  from  time  H  +  Tl  to  infinity.  This  is  the  ex¬ 
posure  as  integrated  from  time  H  +  Tl  or  particle  Impact  time,  which¬ 
ever  is  later,  for  each  Impacted  particle. 

e.  Exposure  accumulated  from  time  H  +  Tl  to  time  H  +  T2.  This  is  the 
e::posure  as  integrated  from  time  H  +  Tl  or  particle  impact  time 
(whichever  is  later)  to  time  H  +  T2.  A  faster  alternative  treatment 
of  accumulated  exposure  not  accounting  for  particle  impact  time  is 
also  provided. 

f.  Total  mass  deposited.  This  is  the  mass  of  fallout,  both  radioactive 
and  Inert,  deposited  on  the  map  grid  points  during  the  entire  fallout 
period. 

g.  Total  mass  deposited  between  times  Tl  and  T?.  This  is  the  total 
mass,  both  active  and  inert,  deposited  during  the  specified  interval. 

h.  Activity  produced  by  a  user  specified  mass  chain  (curies/m^). 


The  undlstorted  map  option 

A  number  of  different  options  exist  with  regard  to  the  scaling  of  output 
maps.  As  stated  previously  (p.  3),  it  is  possible  for  the  Output  Pro  ¬ 
cessor  to  provide  a  numerical  presentation  of  the  fallout  data  on  a  spa¬ 
tially  undistorted  grid.  In  achieving  this  the  user  is  assisted  by  the  pro¬ 
gram.  However,  he  must  supply  the  printer  characteristics  (charac¬ 
ters/inch  both  horizontally  and  vertically)  and  nust  specify  that  the 
"undistorted  map  option"  is  to  be  exercised  Then  the  program  deter¬ 
mines  appropriate  grid  spacings  to  accommodate  the  printer  character¬ 
istics.  On  the  other  hand,  the  tiber  may  specify  the  gri>  intervals  and 
in  so  doing  he  can  obtain  any  rectangular  spatial  distortion  he  desires 
Also,  he  may  allow  the  program  to  make  small  adjustments  to  the 
specified  grid  intervals  to  achieve  a  faster  running  program  if  an  "n- 
distorted  map  is  not  a  requirement.  In  any  case,  an  overall  scale 
factor  for  the  map  must  be  specified. 

Numerical  display  options 

Two  options  exist  at  this  time  for  printing  the  numerical  values  of  the 
fallout  data  over  the  grid  points.  These  options,  which  can  bo  character¬ 
ized  as  the  two-line  E  format  and  the  two-line  F  11.3  format,  are 
explained  and  illustrated,  as  follows,  for  a  single  data  point: 

a.  The  two-line  E  format, 

NNNNNN 

iV.VW  , 

which  is  to  be  interpreted  as 
NNNNNN 

iV.VW  x  10 

b.  The  two-line  F  11.3  format, 

NNNNNN 

iV.VW  , 

which  is  to  be  interpreted  as 


i  NNNNNNV.  VW  . 


In  both  of  these  display  options  the  decimal  point  indicates  the  map  location 
of  the  grid  point. 

In  addition  to  the  two  options  for  printing  numerical  values,  a  third  option 
exists  which  allows  a  numerical  map  image  to  be  recorded  on  a  magnetic 
tape  referred  to  as  the  multiple  burst  tape.  Matched  pairs  of  multiple 
burst  tapes  may  be  processed  by  a  separate  program  (see  Volume  VII  of 
this  documentation)  to  form  printed  or  tape  recorded  maps  of  the  point-hy- 
polnt  sums,  differences,  products,  or  quotients  resulting  from  map 
superposition. 

Sequences  of  Processing  Requests 

The  Output  Processor  has  been  arranged  to  accept  as  input  a  sequence  of  re¬ 
quests  for  processing.  This  was  deemed  appropriate  because  of  the  lai„j  number 
of  different  quantities  which  mignt  be  of  interest  to  the  researcher  now,  and  also 
after  further  development  of  the  program,  and  because  of  the  usual  turn-around 
time  delays  which  plague  the  users  of  computer  centers  operating  with  a  batch  pro¬ 
cessing  system.  Rather  than  handling  requests  on  a  one -per -run  basis,  unlimited 
sequences  of  requests  are  accepted.  Thus,  the  program  is  open-ended  with  respect 
to  requests  on  a  given  grounded-particles  data  set. 

The  control  programming  of  the  Output  Processor  is  designed  to  allow  any 
number  of  maps  to  be  prepared  for  each  of  any  number  of  different  map  limits. 
Suppose  the  user  desires  sets  of  maps  to  be  prepared  separately  for  each  of  two 
different  Bets  of  map  limits.  For  example,  the  user  may  desire  large-scale  maps 
of  essentially  the  entire  local  fallout  field  for  (1)  exposure  rate  normalized  to  H  +  1 
hour,  (2)  total  accumulated  exposure,  and  (3)  activity  from  mass  chain  95,  He  may 
also  desire  the  same  map  options  *•  a  physically  larger  map  covering  a  geograph¬ 
ically  smaller  area  closer-in  to  ground  zero,  To  accomplish  this  Le  can  specify 
the  map  limits  and  scale  fact  f„r  the  large-scale  map  and  follow  it  by  the  needed 
map  option  request  cards,  luese  data  would  be  followed  in  turn  by  the  other  map 
limit  specifications  and  another  series  of  map  option  request  cards. 

Output  Processing  Independent  of  Other  DELFIC  Modules 

In  its  primary  roie  the  Output  Processor  acts  as  the  terminal  portion  of  the 
main  body  of  the  DELFIC  system;  yet,  since  it  consists  of  control  programs  and 
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a  set  of  subroutines,  it  can  also  operate  independently  of  the  other  programs  of  the 
DELFIC  system  except  for  the  Particle  Activity  Module.  (See  the  section  on  Inputs 
to  the  Output  Processor,  "  p,  4 . )  This  feature  can  be  used  to  advantage  if  the  ucer 
saves  the  magnetic  tape  results  of  the  transport  program's  execution.  In  this  way 
the  user  need  not  specify  all  desired  output  at  the  time  of  the  transport  execution  but 
can  make  subsequent  runs  of  the  Output  Processor  as  specific  questions  arise  during 
the  course  of  his  research.  The  U>e  and  card  inputs  to  the  Output  Processor  are 
the  same,  no  matter  which  way  the  program  is  used. 

General  Logic  of  the  Output  Processor 

In  this  section  we  present  a  general  description  of  the  logic  of  the  Output  Pro¬ 
cessor  Including  general  or  organizational  flow  charts.  A  more  detailed  description, 
which  includes  highly  detailed  flow  charts  of  all  subroutines  and  a  complete  discus¬ 
sion  of  tha  more  invoked  s^broutin^?.  is  given  in  the  Program  Details  bullion. 
Particular  emphasis  is  given  to  those  programs  which  we  feel  individual  users  may 
desire  to  modify  for  the  sake  of  adding  new  capabilities  to  the  DELFIC  system.  Ad¬ 
ditional  details  are  also  provided  in  the  Appendix, which  defines  the  various  arrays 
used  by  the  Output  Processor  programs,  and  by  the  glossary  of  primary  program 
variables  which  is  included  in  the  listing  of  the  first  control  program,  subroutine 
LINKS. 

The  Output  Processor  consists  of  two  control  programs  and  eleven  subroutines 
(a  brief  statement  of  the  purposes  of  these  programs  is  set  forth  in  Table  1).  In  ad¬ 
dition,  the  Particle  Activity  Module  subroutines  PAM1  and  PAM2  (see  Volume  V  of 
this  documentation)  are  called  within  the  Output  Processor.  It  should  be  noted  that 
subroutine  PAMl  requires  certain  card  and/or  tape  inputs  to  be  available  during  the 
execution  of  the  Output  Processor. 

Flow  chart  FC-l  gives  a  simplified  picture  of  the  logic  of  the  control  programs, 
LTNX8  and  LINK9,  and  is  a  suitable  point  of  departure  for  the  reader  who  wishes  to 
understand  the  Output  Processor  in  depth.  I r.  the  topmost  diamond  box  of  the  flow 
chart  we  note  that  one  of  the  optional  features  of  the  Output  Processor  allows  the 
user  to  have  the  content  of  the  grounded  particles  tape  printed,  and  have  the  pro¬ 
cessing  terminated  upoa  completion  of  this  or  continued  as  specified.  Following 
downward  one  can  see  the  hierarchical  nan  re  of  the  map  limits  specification  loop 
which  begins  with  the  reading  of  the  coordinate  limits  of  a  map,  and  the  map  request 
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TABLE  1 

OUTPUT  PROCESSOR  PROGRAM  SYNOPSIS 


Program  Name 

Purpose 

LINK  8 

Initialization  and  liaison  with  subroutine  PAMI 
(Particle  Activity  Module)  and  LINKS. 

LINKS 

Interpret  grounded  particles  into  the  output  map 
array  and  call  PAM2  (Particle  Activity  Module) 
for  particle  activities  or  mass  chain  concentrations. 

CALC 

Interpret  grounded  particles  into  the  map  array. 

COUNT 

Select  the  largest  sorted  data  set  for  dumping  onto 
memory  tape. 

CRDP 

Subordinate  control  routine  which  calls  SHIFT  to 
clear  out  most  of  the  particles  array  after  a  pass 
of  the  data  tape  has  been  completed 

DIFTJZl 

Expand  cloud  subdivision  areas  to  account  for 
diffusion. 

LET9GO 

Control  routine  for  the  situation  in  which  sorting 
onto  tape  is  required. 

MAP 

Display  the  Output  (Print  the  Map). 

PROC 

Subordinate  control  routine  which  eliminates  un¬ 
needed  particles,  sorts  and  counts  other  particles, 
and  calls  CALC  to  interpret  those  falling  within 
the  current  map  area. 

RUN1 

Control  routine  for  the  situation  in  which  no  data 
sorting  onto  tape  is  required. 

SHIFT 

Collect  a  selected  set  of  particle  descriptions 
and  write  them  onto  memory  tape. 

SLIDE 

Slide  the  content  of  the  right  buffer  zone  over  to 
the  left  zone  in  preparation  for  processing  the 
next  map  zone. 

ZERO 

Collect  blank  lines  at  the  top  of  the  particles 
arrays  in  preparation  for  reading  in  more  par¬ 
ticle  descriptions. 

loop  which  begins  with  the  reading  of  a  computation  request  (see  the  Sequences  of 
Processing  Requests  section  on  p.  8).  It  should  be  noted  that,  in  order  to  bring 
about  a  final  program  exit,  the  program  must  encounter  a  blank  card  (image)  to  ter¬ 
minate  the  map  request  loop,  and  then  another  blank  card  to  terminate  the  map  limit 
specification  loop.  Most  of  the  program's  complexities  are  found  within  the  two 
lowest  boxes  of  FC-1  and,  consequently,  the  remainder  of  this  discussion  will  be 
devoted  to  them. 


The  user  of  the  Output  Processor  must  specify  the  area  that  he  wishes  to  be 
mapped  by  indicating  its  limiting  coordinates.  This  area  is  rectangular  with  its  sides 
aligned  north-south  and  east-west.  The  positive  Y  direction  points  north  and  the  posi¬ 
tive  X  direction  points  easi.  Grounded  cloud  subdivisions  are  represented  on  the  data 
tape  by  the  impact  coordinates  of  one  central  particle  for  each  subdivision.  The 
shape  of  each  cloud  subdi  vision  is  a  square  and  the  length  of  the  side  of  this  square 
at  subdivision  definition  tiuie  is  communicated  to  the  Output  Processor  by  the  trans¬ 
port  program  via  the  parameter  BZ  on  the  grounded  particles  tape.  Since  central 
particles  which  fall  within  the  distance,  BZ/2,  from  the  edge  of  the  specified  map 
limits  will  affect  the  map  area,  we  must  consider  all  particles  falling  within  a  rec¬ 
tangle  which  includes  a  margin  of  v/idth  BZ/2  on  each  aide  of  the  map.  Figure  1  indi¬ 
cates  both  the  map  specified  by  the  user  (heavy  rectangle)  and  the  area  of  computa¬ 
tional  late  rest  (the  light  rectangle  surrounding  the  heavy  one).  Note  that  the  symbol¬ 
ogy  of  Figure  1  corresponds  to  that  of  th  FORTRAN  programs. 


XO,  YF 


XF,  YF 


1.  The  Map  Area  Without  Subdivision 


The  overall  map  area  hi  subdivided  into  a  grid  according  to  the  grid  intervals 
specified  by  the  user  (as  modified  by  the  program  if  the  undistorted  map  option  is 
specified  (see  p.  7)).  If  the  grid  intervals  specified  by  the  user  are  such  that  the 
data  for  the  entire  map  can  be  contained  completely  in  the  core  storage  map  array, 
no  further  subdivision  of  the  map  area  need  be  made.  This  Is  the  situation  repre¬ 
sented  by  Figure  1.  On  the  other  hand,  if  either  smaller  grid  intervals  or  a  larger 
map  (or  both)  are  specified,  the  number  of  gild  points  In  the  complete  map  may  be 
too  large  to  fit  ir  the  map  array  at  one  time.  In  this  situation  subdivision  of  the 
kind  shown  in  Figure  2  may  be  required.  This  type  of  subdividing  is  done  automat¬ 
ically  by  the  program  without  any  guidance  from  the  user.  The  snaded  strips  in 
Figure  2  represent  buffer  zones  between  map  areas  which  must  be  used  since  cloud 
subdivisions  falling  near  the  boundary  between  map  zones  could  otherwise  affect 
more  than  one  z ore.  *  The  region  marked  "Zone  1 "  will  be  treated  as  a  first  map 
and  any  particles  falliD*,  into  (or  near)  this  zone  will  be  immediately  interpreted  into 
the  map  array.  Particles  falling  in  the  other  marked  zones  (Zone  2,  Zone  3,  etc.) 
will  be  eventually  written  onto  zone  memory  tapes  for  later  interpretation  into  the 
map.  However,  since  particles  falling  into  the  buffer  /.ones  affect  more  than  one 
map  area,  they  must  be  written  onto  more  than  one  tape.  The  tape  memory  assign¬ 
ments  of  areas  are  indicated  in  Fig-ire  2  by  the  overlapping  arrows  marked  Tape  1," 
"Tape  2,  *  and  "Tape  3."  In  the  situation  represented  in  this  figure  there  will  be, 
after  the  first  map  area  is  printed,  effectively  a  one-to-one  correspondence  between 
subsequent  map  zones  and  tape  memory  data  sources.  A  slight  exception  to  this 
statement  is  brought  about  by  the  need  1.0  keep  the  'ighi-hand  buffer  zone  map  area 
in  the  map  array  *o  account  for  the  effect  of  Zone  S  particles  on  Zone  N  +  1. 

Figure  2  also  illustrates  the  meaning  of  tne  program  vaviab'es  X1.X2.X3.  and 
X4  for  the  prc^essing  of  the  second  map  zone.  XI  X2,  X3,  and  X4  denote  respectively 
the  X  coordinates  of  the  left  side  of  the  left  buffer  zone,  the  right  side  of  the  left  def¬ 
er  zone,  the  left  side  of  the  right  buffer  zone,  and  the  right  side  of  the  right  buffer 
zone.  These  variables  are  set  by  the  program  prior  to  the  start  of  processing  on  s 
new  map  zone  DELTAX,  which  ir  also  illustrated,  is  a  constant  set  the  program 
on  the  basis  of  the  required  map  dimensions  and  the  size  of  liK-'  v  at  latle  map  array. 

Buffer  zones  and  map  zones  should  not  be  confused  with  the  printer  strips  ,nio 
which  a  map  zone  is  actually  divided  (or  printing  and  subsequent  assembly.  The 
maximum  width  oi  printer  strips  is  fixed  by  the  number  of  characters  that  the  printer 
can  print  on  on<*  ltne. 
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If  there  were,  for  example,  only  two  tapes  available  for  sorting  instead  of  the 
three  illustrated  in  Figure  2,  it  would  be  necessary  to  use  a  tape  as  an  overflow 
memory.  In  this  situation  the  particles  falling  into  Zones  3  and  4  would  be  written 
onto  the  overflow  tape,  and  after  the  interpretation  of  the  data  on  Tape  1,  the  pro¬ 
gram  would  return  to  re-sort  the  data  on  the  overflow  tape.  The  use  of  this  overflow 
tape  memory  effectively  makes  the  program  entirely  open-ended  in  regard  to  the 
size  of  the  maps  it  can  produce.  At  the  same  time  the  use  of  tape  memory,  when 
(but  only  when)  it  is  required,  should  tend  to  make  the  program  acceptably  efficient 
for  smaller  tasks. 

It  is  important  to  note  that  the  "zones"  depicted  in  Figure  2  should  not  be  con¬ 
fused  with  the  strips  of  printer  paper  that  are  attached  side-by-side  for  direct  con¬ 
touring.  (See  the  discussion  on  p.  16.)  The  number  of  pa^er  strips  produced  is 
determined  by  the  interrelation  of  the  number  of  data  points,  the  grid  interval  in  the 
X  direction,  and  the  width  of  the  printer  line.  In  general  there  will  be  more  than 
one  strip  for  each  zone. 


XC,YF  XF.  YF 

XMIN.YMAX  XI  X2  X3  X4  XMAX,  YMAxl 


Buffer  Zones 


Figure  2.  Map  Area  With  Subdivision 
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AN  ILLUSTRATION  OF  OUTPUT  PROCESSOR  USE 

For  the  sake  of  illustration,  let  us  assume  that  you  as  a  researcher  are  prepar¬ 
ing  to  make  use  of  the  Output  Processor.  How  should  you  start  and  what  should  your 
strategy  be? 

It  is  suggested  that  the  first  things  to  be  done  in  an  unfamiliar  prediction  situa¬ 
tion  are  to  have  the  Output  Processor  print  a  hard  copy  of  the  grounded  particles 
tape,  and  in  addition,  to  call  for  a  large-scale  map  to  provide  an  overview  of  the 
fallout  field.  Both  of  these  tasks  can  be  carried  out  in  the  execution  run  that  does 
the  transporting.  The  grounded  particles  tape  produced  by  the  Transport  Module 
and  used  by  the  Output  Processor  in  the  initial  run  may  be  saved  and  reused  in  sub¬ 
sequent  runs  of  the  Output  Processor  and  Particle  Activity  Modules  to  produce  any 
additional  maps  desired.  In  order  to  specify  the  first  set  of  tasks  the  user  must  set 
about  a  dozen  control  parameter  values.  One  parameter  (IC(18)  >  0,  p.49)  causes 
the  grounded  particles  tape  to  be  printed.  Another  (IC(17)  =  0,  p.49)  causes  proces¬ 
sing  to  continue  after  that  printing  has  been  finished.  Four  others  (XMIN,  XMAX, 
YMIN,  YMAX,  p.49)  give  the  coordinate  limits  of  the  desired  map.  Two  others 
(DGX  and  DGY)  give  the  map  point  grid  intervals  in  the  X  and  Y  directions.  One 
other  specifies  that  the  map  should  be,  for  example,  a  map  of  exposure  rate  "nor¬ 
malized"  to  H  +  1  hour  (NREQ  =  2,  p.51). 

In  the  second  and  subsequent  execution  runs  after  the  printed  list  of  grounded  par¬ 
ticles  and  the  large-tjcale  mr  p  have  become  available,  the  user  may  request  that  any 
number  of  more  detailed  maps  be  printed.  These  larger  and  more  precise  maps  can, 
of  course,  portray  any  of  the  possible  output  quantities  of  the  DELFIC  system. 

Figure  3  illustrates  how  a  user  may  arrange  his  requested  map  areas  and  scale 
factors  to  expose  prediction  details  of  interest  to  him.  Area  coverages  and  map  dimen¬ 
sions  shown  are  merely  illustrations  ar.d  in  no  sense  are  meant  to  imply  any  restric¬ 
tions  in  the  use  of  die  Output  Processor  since  their  characteristics  are  completely 
under  the  control  of  the  user.  Map  1  is  a  large-scale  overview  that  indicates  the  shape 
and  location  oi  the  fallout  field,  but  is  necessarily  crude  because  it  is  small  in  actual 
size  (2x2  ft)  and  thus  contains  a  small  number  of  points.  Maps  2,  3,  and  4  represent 
much  more  precisely  the  predicted  fallout  field  within  10  miles  of  ground  zero.  They 
are,  for  example,  6x8  ft  in  size  and  contain  nine  times  as  many  data  points  as  the 
overview  (Map  1),  but  represent  an  area  of  the  fallout  field  less  than  one  tenth  as  large 
as  the  overview.  They  may,  for  instance,  portray  predictions  of  mass  deposition, 
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MAP  6,  4x2  ft 


MAP  5,  4x2  ft 

MAPS  2,  3,  4 
Close-In 
Details,  6x6  ft 

MAP  1 

Overview.  2x2  ft 

exposure  rate  at  H  +  2  hours,  and  the  concentration  of  mass  chain  89.  Maps  5  end 
6  continue  the  representation  of  mass  deposition  but  do  so  with  less  precision  than 
that  of  Map  2  since  they  are  about  4  ft  long  and  cover  larger  fallout  areas. 

The  task  of  arranging  the  card  inputs  to  the  Output  Processor  is  not  a  complex 
oue,  but  it  does  require  that  the  user  make  certain  decisions  about  what  he  wishes  to 
have  portrayed  and  how  he  wants  it  to  be  portrayed.  He  should  know  what  coordinates 
were  used  to  identify  the  location  of  ground  zero  within  the  transport  module.  He 
should  know  that  all  distances  (coordinates)  are  measured  in  meters  from  the  same 
coordinate  origin.  He  should  have  at  least  a  rough  idea  of  the  direction  of  the  winds. 
Beyond  that,  he  need  only  know  which  of  the  display  options  he  wants  to  have  portrayed. 
A  card-by-card  and  parameter-by-parameter  explanation  of  the  input  card  deck  for 
the  Output  Processor  is  given  in  the  User  Information  section. 

The  map  produced  by  the  Output  Processor  will  consist  of  a  sequence  of  numbered 
"strips"  of  printer  paper  that  can  be  detached  at  the  boundaries  between  successive 
strips  and  assembled  side-by-side  into  a  single  map  of  the  overall  area  covered  by 
the  specified  map  limits.  When  so  assembled  (the  strips  are  numbered  in  sequence 
from  left  to  right)  and  hung  on  the  wall  for  viewing,  the  data  point  with  minimum  X 
and  minimum  Y  coordinates  will  be  found  in  the  lower  left  hand  corner  of  the  map(i.e. , 
the  lower  left  hand  corner  of  strip  number  one).  The  coordinates  of  this  point  will  be 
(XMIN  +  DGX,  YMIN  +  DGY).  This  point  need  not  be  either  the  origin  of  coordinates 
or  ground  zero. 
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PROGRAM  DETAILS 


Control  Program  LINK 8  (FC  -2) 

The  purpose  of  LINK8  is  to  (1)  initialize  the  Output  Processor  system,  (2)  print 
out  an  impacted  paiticle  list,  and  (3)  call  PAM1  of  the  Particle  Activity  Module. 

(See  Volume  V  of  this  documentation.) 

Upon  entrance,  LINK8  first  sets  a  number  of  program  constants  which  denote 
such  things  as  the  logical  identifiers  of  certain  essential  tapes  and  the  maximum 
sizes  of  certain  arrays.  When  installing  the  Output  Processor  at  a  new  computation 
center  or  on  a  new  computer  system,  it  is  essential  that  the  system  ta^e  identifiers 
be  checked  to  see  that  they  conform  to  established  requirements.  Note,  however, 
that  changing  the  tape  number  assignments  at  the  beginning  of  LINK 8  will  suffice  for 
all  programs  of  the  Output  Processor  since  all  tape  references  have  been  made  via 
established  tape  names,  e.g. ,  ISIN,  ISOUT,  IPOUT. 

In  a  similar  way,  whenever  it  becomes  possible  to  run  the  Output  Processor 
on  a  computer  that  has  a  high  speed  memory  larger  than  32,  768  words,  adjustments 
should  be  made  to  the  sizes  of  certain  arrays  within  the  program.  If  a  larger  memory 
is  available,  the  programs  will,  of  course,  function  correctly  without  adjustment, 
but  improved  program  efficiency  can  be  easily  attained  if  adjustments  are  made.  If 
the  program  must  be  run  in  a  smaller  segment  of  memory  than  it  now  uses,  only  a 
few  simple  adjustments  must  be  made  !n  order  to  scale  down  the  program,  but  pro¬ 
gram  efficiency  will  necessarily  suffer.  The  program  array  OMAP  and  the  compan¬ 
ion  array  size  variable  NMAP  may,  for  instance,  be  adjusted  downward  (or  upward) 
together  to  change  the  size  of  the  usable  map  array.  Changing  the  OMAP  dimension 
statement,  and  the  statement  which  sets  the  value  of  NMAP,  is  all  that  one  needs  to 
do  since  all  usage  of  the  map  array  is  based  on  the  parameter  NMAP.  Of  course, 
the  new  size  of  the  OMAP  array  as  specified  in  the  dimension  statement  and  the  value 
of  NMAP  should  coincide,  and  all  Output  Processor  programs  should  be  recompiled 
(with  adjustments  made  io  all  dimension  card  decks). 

After  parameter  setting,  LINK8  attempts  to  check  the  identification  of  the 
grounded  particles  tape,  and  if  successful,  it  proceeds  to  read  and  record  all  pre¬ 
vious  identifier  records  that  are  on  that  tape.  Next,  particle  size  frequency 
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distribution  data  are  read  into  memory  from  the  grounded  particles  tape.  This  is 
followed  by  a  list  of  identifiers  of  tapes  available  for  use  in  sorting  which  is  read 
from  the  system  input  tape.  A  partial  check  is  then  made  to  remove  from  this  list 
any  inadvertent  references  to  system  tapes  or  the  grounded  particles  tape.  A  du¬ 
plicate  copy  of  the  corrected  list  is  stored  in  the  array  IIT(J)  for  eventual  recon¬ 
struction  of  the  entrance  condition.  Next,  the  overall  control  array  IC(J)  is  read 
and  a  branch  is  made  on  the  value  of  IC(18)  which  indicates  whether  or  not  a  hard 
copy  of  the  grounded  particles  tape  is  desired  at  this  time.  After  producing  a  hard 
copy  of  the  grounded  particles  tape,  a  test  of  IC(17)  indicates  whether  or  not  further 
processing  is  desired,  and  if  it  is,  the  program  execution  proceeds  to  statement 
number  5111. 

At  statement  5111  printer  character  spacing  data  (IH  and  IV)  are  read  and  a 
title  page  is  written  onto  the  system  output  tape.  Then  control  is  passed  to  sub¬ 
routine  PAM1  of  the  Particle  Activity  Module  where  occurs:  (1)  sorting  and  editing 
of  nuclear  transitions  and  decay  subchain  data,  (2)  computation  of  Freiling  FR  fac¬ 
tors,  and  (3)  computation  of  part  1  of  the  induced  activity  contribution.  Upon  return 
from  PAM1,  control  is  transferred  to  LINK9. 

Control  Program  LINK9  (FC-3) 

Under  the  control  of  program  LINK9  of  the  Output  Processor,  output  maps  are 
prepared  and  printed  and  subroutine  PAM2  of  the  Particle  Activity  Module  is  called 
as  needed  to  compute  particle  activities  or  specific  mass  chain  concentrations  (see 
Volume  V).  At  statement  119  a  sum  Of  map  ordinates  is  printed  if  one  was  Com¬ 
puted  earlier.  At  statement  1191  a  specification  of  the  limiting  coordinates,  grid 
intervals,  and  surface  roughness  factor  of  a  desired  map  are  read.  If  the  surface 
roughness  factor  is  input  as  zero,  it  is  set  to  1.0.  The  sum  of  the  grid  intervals 
is  tested  for  zero  as  a  termination  condition.  This  is  the  correct  exit  condition 
and  leads  to  an  on-line  printer  comment  for  the  notification  of  the  computer 
operators. 

If  acceptable  grid  intervals  are  specified,  the  local  control  array  JC(J)  is  read, 
parameters  are  set  for  the  processing  request  loop,  and  a  test  is  made  on  the  printer 
description  parameters  IV  and  IH. 


Statement  1211  is  a  return  point  where  a  special  terminal  record  is  put  onto 
the  multiple  burst  tape  if  one  is  in  use.  At  1209  the  sum  oi  map  coordinates 
is  printed  if  it  has  not  been  printed  previously.  At  statement  1219  we  enter  the 
processing  request  loop  and  read  a  request  in  the  form  NREQ,  Tl,  T2,  MASCHN. 
After  input  checking  and  loop  initialization,  a  test  is  made  on  NREQ  for  the  exit 
condition  (NHEQ  *»  0)  which  leads  to  a  return  to  the  processing  loop  at  statement 
119.  If  a  valid  request  number  is  encountered,  tape  IP  OUT  is  repositioned  so 
that  the  next  read  statement  can  bring  in  grounded  particle  data;  a  request  title  page 
f'  ,nueu,  iWii  .  ‘~ance  parameters  are  set  for  both  the  Output  Processor  and  sub¬ 
routine  PAM2.  If  JC(18)  and  JC(16)  so  Indicate,  an  adjustment  of  grid  intervals  is 
made  so  that  an  undistorted  map  will  be  produced,  but  if  an  undistorted  map  is  not 
required  and  small  grid  interval  adjustments  are  permitted,  a  transfer  to  1301  is 
made  so  that  small  adjustments  can  be  made  to  the  grid  Intervals  to  yield  a  more 
rapid  program  execution. 

At  statement  140  final  grid  Intervals  have  been  arrived  at  and  the  width  of  the 
buffer  zones  is  set  equal  to  the  width  of  an  integral  number  of  grid  intervals.  Com¬ 
putations  are  then  made  of  the  numbers  of  grid  points  covered  by  the  map  in  its  two 
principal  directions.  NZ,  the  number  of  memory  map  zones,  is  also  computed  here. 
NZ  is  one  less  than  the  total  number  of  zones  into  which  die  map  must  be  divided  so 
that  it  can  be  produced  within  the  available  map  array.  Next,  NOX,  the  number  of 
map  grid  intervals  between  buffer  zones,  and  DELTAX,  the  distance  between  buffer 
zones,  are  computed  for  later  use.  Than,  if  grin  interval  adjustments  were  made, 
a  record  is  made  as  part  of  the  program’s  printed  output. 

At  statement  1405  more  initialization  is  performed.  At  300,  parameters  MIN, 
JIN,  and  LAST  are  set  on  the  basis  of  die  aw.iber  of  available  tape  units.  If  NZ  is 
zero,  only  one  map  area  must  be  computed,  i.e. ,  all  of  the  required  map  will  fit 
within  the  map  array  at  one  time.  Note  that  &  Single  pass  through  RUNi  and 
MAP  is  required.  Return  is  made  to  1211  which  leads  to  the  reading  of  the  next 
request.  If  NZ  is  positive  but  less  than  the  number  of  tapes  available,  a  single 
sorting  pass  followed  by  a  sequence  of  Single  area  interpretations  is  called  for. 

Also,  subroutine  LETSGO  is  called  to  process  the  first  map  area  and  sort  the 
remaining  particle  data  onto  available  tapes.  Then  a  loop  which  c;  is  SLIDE,  RUNI, 


and  MAP  la  executed  to  interpret  the  sorted  data  and  complete  the  printed  map. 

This  time  return  is  also  made  to  1211  for  the  uext  request.  If  NZ  is  greater  than  the 

number  of  available  tapes,  more  than  one  sort  pass  is  required.  Again 
* 

LETSGO  interprets  the  first  map  area,  sorts  data  for  a  number  of  other  adjoining 
a  re  as.  and  writes  die  remaining  data  onto  a  separate  overflow  tape  for  subsequent 
sorting.  Return  in  this  case  is  made  to  statement  300  so  that  a  second  sort  of  the 
overflow  tape  may  be  carried  out.  Eventually,  a  return  for  the  next  request  will  be 
made  via  one  of  the  previously  mentioned  transfers  to  statement  1211. 

Subroutine  CALC  (FC-4) 

The  purpose  of  this  subroutine  is  to  interpret  a  cloud  subdivision  into  the  map 
array  OMAP.  This  Interpretation  consists  of  first  the  selection  of  the  appropriate 
computation  code  which  is  shown  on  the  flow  chart  as  the  branching  operation  based 
on  the  value  of  the  parameter  NREQ. 

After  the  branch,  a  map  ordinate  increment  is  computed  and  stored  in  the  vari¬ 
able  F.  Next,  starting  at  statement  number  100,  CALC  computes  appropriate  stor¬ 
age  indices  for  the  control  of  the  program  loops  which  actually  carry  out  the  map 
array  Incrementing  operation. 

Last  die  increment  F  is  added  to,  or  compared  with,  the  concert  of  the  selected 
map  array  points  by  means  of  the  two  nested  FORTRAN  do  loops  starting  just  beyond 
statement  number  19. 

It  should  be  noted  that  code  insertion  points  have  been  provided  within  the  com¬ 
putation  selection  branching  operation  so  that  new  computational  options  may  be 
easily  added  to  the  program.  To  insert  the  first  new  code  one  need  only  make  sure 
that  the  code  insert  begins  with  statement  number  109  and  ends  With  the  statement 
'Go  to  100.  "  The  user  must  put  the  insert  in  the  place  of  the  card  *109  CONTINUE" 
and  then  recompile  subroutine  CALC.  Of  course,  the  insert  must  not  contain  any 
statement  numbers  previously  used  within  CAL C. 

Some  explanation  of  the  procedure  used  to  compute  the  storage  control  indices 
is  called  for  since  this  procedure  is  somewhat  involved  and  very  central  to  the 
Output  Processor.  As  noted  previously,  the  basic  shape  of  the  cloud  subdivision 
in  plan  view  is  a  square.  The  length  of  the  side  of  this  square  is  communicated  to 
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the  Output  Processor  and  to  subroutine  CALC  in  the  parameter  BZ,  The  procedure 
under  discussion  begins  at  statement  number  100  with  the  calculation  of  the  coordi¬ 
nates  of  the  corners  of  the  grounded  cloud  subdivision  square.  These  coordinates 
are  stored  in  the  variables  WXL,  WXR,  WYB,  WYT.  Any  map  grid  point  which 
lies  within  this  square  will  have  the  quantity  F  added  to  its  value.  Next,  the  proce¬ 
dure  computes  NOL,  die  X  index  of  the  grid  point  which  lies  just  to  the  rignt  of  the 
square’s  left  boundary,  and  JX,  the  X  index  of  the  point  which  lies  just  to  the  right 
of  the  square's  right  boundary.  NWX,  the  difference  between  these  two  indices. 
Indicates  how  many  grid  points  are  to  be  incremented  in  the  X  direction.  The  result 
of  this  subtraction  may  be  zero,  in  which  case  no  grid  ooint  will  be  incremented  by 
CALC.  A  similar  treatment  is  given  to  the  Y  dimension,  and  the  final  set  of  results 
is  stored  in  the  variables  NOL  and  NWX  for  the  X  direction,  and  NOB  and  NWV  in 
the  Y  direction. 

At  statement  number  19,  loop  control  variables  K,  MM,  and  NN  are  computed 
in  order  to  ridge  the  logical  gap  between  the  previously  computed  two-dimensional 
indices  such  as  NOL  and  the  one-dimensional  array  in  which  the  map  ordinates  are 
really  stored. 

Certain  complexities  sre  added  to  the  foiegoing  by  the  grid  interval  adjustment 
options  which  are  available  and  by  the  need  to  correctly  handle  boundary-effects 
situations  in  which  a  cloud  subdivision  falls  partially  off  the  map  area.  The  details 
of  these  treatments  are  provided  oy  the  program  listings. 


Subroutine  COUNT  (FC-5) 

This  subroutine  selects  the  currently  largest  sorted  and  identified  set  of  par¬ 
ticle  descriptions  for  dumping  onto  an  external  (tape  )  memoir  device.  In  making 
this  selection  use  is  made  of  the  particle  counter  array,  NP(  ),  in  which  NP(J)  for 
even  J  contains  the  count  of  the  number  of  cloud  subdivision  central  particles  that 
have  fallen  into  the  J/2th  map  zone  (exclusive  of  adjoining  buffer  zones).  The 
NP(J-l)  and  NP(J  ♦  1)  for  even  J  contain  the  counts  of  centra!  particles  falling  into 
the  buffer  zones  to  the  left  and  right  (west  and  east)  of  the  Jth  interbuffer  map  zone. 
In  general,  MAX(NP(J))  for  even  J  represents  the  maximum  number  of  particle 
descriptions  that  can  be  sct"illj  removed  from  the  particle  arrays  by  a  single 
dumping  operation. 


Since  cloud  subdivisions  whose  central  particles  fall  into  a  buffer  zone  affect 
the  two  map  zones  on  either  side,  it  is  necessary  to  dump  buffer  zone  particles 
along  with  the  particles  from  the  surrounding  map  zones.  However,  buffer  zone  par¬ 
ticle  descriptions  must  also  remain  in  memory,  since  they  also  affect  some  other 
map  zone.  (The  subroutine  shift  reidentifies  these  duplicate  particle  descriptions 
by  changing  their  classifications  as  recorded  in  the  KTR(  )  array  in  order  to  avoid 
double  counting  later.)  In  general,  the  number  of  particle  descriptions  which  Will 
actually  be  dumped  onto  memory  tape  exceeds  MAX(NP(J)  for  even  J)  by  NP(J  - 1) 

+  NP(J+  1).  A  slight  exception,  however,  must  be  made  for  the  first  map  zone 
since  the  left-reaching  effects  of  the  particles  in  its  left  buffer  zone  have  already 
been  taken  into  account  (in  the  in-core  map  array)  and,  therefore,  the  particles 
falling  into  its  left  buffer  zone  need  not  be  duplicated. 

After  selecting  MAX  [  NP(J)  for  even  J,  NP(1)  +  NP(2)] ,  subroutine  COUNT 
adds  the  appropriate  buffer  zone  counters  to  compute  the  number  of  particle  de¬ 
scriptions  to  be  dumped  and  puts  the  result  in  parameter  NS.  Next  it  sets  the 
parameter  KTAPE  equal  to  the  identification  number  of  the  tape  on  which  the  forth¬ 
coming  dump  operation  (subroutine  SHIFT)  should  write  and  then  returns. 

Subroutine  CROP  (FC-6) 

The  purpose  of  subroutine  CRDP  (core  dump)  is  to  empty  the  particle  arrays 
of  all  particle  descriptions  except  those  that  can  be  processed  into  the  next  load  of 
the  map  array,  OMAP(  ).  CRDP  is  actually  a  control  program,  since  it  uses  sub¬ 
routine  SHIFT  to  do  the  writing  of  sets  of  particle  descriptions  onto  appropriate 
map  zone  memory  tapes.  (SHIFT  does  not  actually  clear  (store  zeros  in)  all  words 
of  each  affected  particle  description  in  the  particle  arrays,  but  merely  sets  the  ap¬ 
propriate  KTR(J)  entries  to  zero  to  indicate  that  the  Jth  lines  in  the  particle  arrays 
are  available  for  reuse.)  CRDP  also  makes  the  adjustments  required  to  keep  the 
particle  counter  array  NP(  )  current.  It  should  be  noted  that  CRDP  does  not  dump 
the  descriptions  of  particles  falling  into  the  first  sort  zone  (J  =  2)  onto  tape.  These 
particles  are  left  in  the  memory  array  for  immediate  processing  into  the  next  printed 
map  zone.  Note  also  that  CRDP  writes  a  final  zero  on  each  sort  tape  as  an  indica¬ 
tion  of  the  end  of  the  data  on  the  tape . 
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Subroutine  LETSGO  (FC-7) 


This  subroutine,  like  RUN1,  is  a  specialized  executive  program,  but  it  may 
call  subroutines  COUNT,  PROC,  SHIFT,  and  ZERO.  It  is  used  onlv  when  the  data 
on  the  current  input  tape  must  be  sorted  because  the  required  output  map  is  too 
large  to  fit  completely  into  high  speed  memory  at  one  time.  LETSGO  first  sets  NE, 
the  counter  of  empty  lines  in  the  particle  arrays,  at  the  full  size  of  the  arrays,  and 
then  clears  the  zone  counters  NP(I)  and  the  particle  zone  indicators  KTR(I).  Next, 
it  reads  a  particle  block  count  NIJ  and  tests  NU  for  the  termination  condition 
(NU=0),  In  the  event  of  termination  CRDP  is  called  to  see  that  the  content  of  the 
particle  array  is  written  out  onto  the  appropriate  tape  (if  required).  If  NU  is  posi¬ 
tive,  other  tests  are  performed  and,  if  necessary,  COUNT  and  SHIFT  are  called  to 
make  room  for  the  incoming  data  block.  Thereafter,  except  for  the  first  pass, 

ZERO  is  called  to  group  the  number  of  empty  particle  array  lines,  required 
for  the  incoming  particle  block. 

The  particle  block  Ir.  then  read  into  the  computer  and  subroutine  PROC  is  called 
to  process  the  particle  (or  cloud  subdivision)  data.  This  processing  consists  of 
determining  which  numbered  map  zone  or  buffer  zone  the  central  particle  falls  into 
and  recording  that  zone's  identification  number  in  the  central  particle's  zone  indi¬ 
cator  parameter  KTR(  ).  At  this  time  PROC  also  uses  CALC  to  process  into  the 
map  array  OMAP(  )  all  those  central  particles  that  affect  the  map  area  currently 
being  evaluated. 

Subroutine  MAP  (FC-8) 

This  subroutine  writes  complete  fallout  maps  on  the  system  output  tape  ISOUT 
for  batch-printing.  It  writes  a  map  title,  a  description  of  what  quantity  the  map 
portrays,  and  an  indication  of  the  map's  style  of  presentation  (format).  It  divides 
the  output  map  into  printer  strips  on  the  basis  of  the  printer  width  parameter  INC. 

It  prints  a  strip  count  (MAPRUN)  at  the  top  of  each  strip  for  identification  purposes. 
Since  an  individual  map  may  consist  of  more  than  one  map  array  full  of  data,  it  is 
necessary  that  MAP  operate  correctly,  even  if  an  individual  map  must  be  produced 
by  a  sequence  of  calls  to  MAP.  This  feature  is  facilitated  by  the  parameter 
MAPRUN,  which  is  zero  at  the  time  of  the  initial  call  of  subroutine  MAP  and  is 
positive  thereafter. 


Following  the  detailed  flow  chart  of  subroutine  MAP,  we  see  at  the  beginning  a 
transfer  on  the  basis  of  MAPRUN  to  a  first, -pass  code  if  MAPRUN  equals  zero.  In 
this  first-pass  code  parameter  initializations  are  performed,  a  map  title  is  written, 
the  display  option  control  parameter,  JC(1),  Is  checked  for  an  acceptable  value,  and 
then  a  branch  transfer  is  made  to  a  code  that  writes  the  presentation  style  title  and 
makes  control  transfer  assignments  within  the  map  writing  loops. 

Between  the  statement  numbers  102  and  170  a  two-part  title  describing  the  quan¬ 
tity  presented  in  the  map  is  written.  Between  statement  numbers  170  and  initial¬ 
izations  are  made  for  the  three  nested  map  writing  loops.  At  the  time  when  P^  is 
first  reached  M  contains  the  number  of  printer  strips  that  are  to  be  produced,  and 
LEFT  has  the  number  of  columns  that  should  appear  on  the  last  printer  strip. 

At  P  ,  which  is  the  return  point  for  the  outer  map  writing  loop  (printer  strip 
loop), MAPRUN,  the  couuter  of  printer  strips, is  incremented  and  strip  title  is  written. 
Also,  KL,the  lower  index  for  retrieval  from  the  one-dimensional  map  array  OMAP, 
is  set  at  its  initial  value.  Note  that  in  the  iteration  KL  progresses  from  its  largest 
value  to  its  smallest  value  to  invert  the  map  which  is  stored  numerically  inverted 
in  the  map  array. 

At  P2,  the  return  point  for  the  middle  map  writing  loop  (printer  line  loop),KH, 
the  upper  index  for  retrieval  from  the  map  array,  is  set  and  KDC.an  index  for  the 
printer  line  integer  array  JMAP.is  initialized. 

At  P_,  the  return  point  for  the  last  map  writing  loop  (data  point  loop),  KDC  is 

u 

incremented  and  a  transfer  is  made  to  the  desired  presentation  code  c"  the  basis  of 
previous  assignment.  The  two  furnished  printer  display  codes  take  their  inputs 
from  the  map  array  and  place  their  results  back  into  the  map  array  and  into  the 
integer  printer  line  array  JMAP.  All  map  producing  codes  return  to  statement  num¬ 
ber  300. 

Below  300  the  printer  lines  are  written  onto  the  output  tape,  suitable  indexing 
operations  are  performed,  and  return  is  made  to  deal  with  either  the  next  line  in 
the  current  strip  or  the  first  line  (and  title)  on  the  next  strip,  or  a  final  return  is 
made  to  the  calling  program.  Note  that  if  entrance  is  made  to  MAP  with  MAPRUN 
set  positive  as  a  consequence  of  a  previous  entrance,  the  overall  titles  will  not  be 
printed  again  and  strip  counting  will  be  resumed  where  it  had  been  left  off. 


33 


MAPBUN 


BV BUQDtMT  f1#1\ 
pa ns 


TOUT  PAM 
(MC  ■  U 

KKL>  NBZXI  ♦  1 

NX  -  xxkaf  - 1  •  ntn 

Ltrr  -  no (ax  -  m  /  dock. 

TM04.  fl  .DQX 
XCOOfU>XMD*DOX 
VD*OWC 
xrmrv  v»c»DGX 


'non 

MAP 

iTTO/ 


This  subroutine  is  called  by  LETSGO  to  sort  (classify'  panicle  descriptions, 
and  it  calls  CALC  to  actually  interpret  into  the  map  array  those  particles  which 
affect  the  current  map  image.  If  the  value  of  ICTR  is  zero  upon  entrance,  a  trans¬ 
fer  is  made  to  a  code  which  discards  part,  le  descriptions  that  fall  outside  the  user's 
area  of  interest  and  consolidates  thote  remaining  at  the  top  of  the  particle  array. 

If  ICTR  is  not  zero,  this  code  is  bypassed  and  a  transfer  is  made  immediately  to  a 
code  which  assigns  values  to  the  KTR  array  for  tho**e  particle  descriptions  which 
cannot  be  immediately  interpreted  into  the  map  array,  and  CALC  is  called  to  inter¬ 
pret  those  remaining. 


The  value  assigned  to  the  KTR  for  each  particle  description  indicates  into  which 
map  zone  or  buffer  zone  the  particle  has  fallen.  This  value  is  computed  by  the  loop 
beginning  at  the  point  labeled  o  where  the  east-west  distance  between  the  particle 
position  and  XO,  the  eat  er.i  boundary  of  the  extended  area  of  interest  (see  Figure 
2)  is  computed  and  stored  in  R.  Thereafter,  the  program  determines  the  classifica¬ 
tion  parameter  KTR  for  each  particle  by  alternately  subtracting  the  width  of  a  buffer 
zone  (BZ2),  and  the  width  of  a  map  zone  exclusive  of  buffer  zones  (DELTAX),  from 
R.  and  performing  a  test  against  zero  after  each  subtraction.  While  this  loop  of 
subtractions  is  being  carried  out,  the  parameter  J  is  used  to  record  the  number  of 
the  map  zone  counter  array  element  pertinent  to  the  particle.  Note,  nowever,  that 
central  particles  falling  into  the  area  of  the  in-core  map  array  can  be  processed 
immediately  into  the  map  array  by  a  call  to  CAI.-C  and,  therefore,  need  not  be  stored 
permanently  in  the  particle  arrays.  Thus,  their  classification?  .eed  not  be  stored 
in  their  KTR(J)  which  remain  zero  It  should  also  be  noted  that  the  indexing  of  die 
counter  array  NP(  )  is  offset  by  two  from  the  classification  index  ,j .  This  leads  to 
a  usage  of  NT(J)  as  the  counter  for  the  number  of  central  particles  failing  into  the 
first  buffer  zone  to  the  right  of  the  first  printed  map  area,  and  KP(2)  as  the  counter 
for  the  first  map  zone  to  the  right  of  the  first  printed  map.  and  so  forth. 


Subroutine  Rl'Nl  (FC-10) 


This  subroutine  is  a  specialized  executive  program  which  calls  only  subroutine 
CALC.  It  is  used  in  the  situation  where  only  one  pass  of  the  input  data  tape  is  re- 
quired  in  order  to  fully  account  for  the  data  on  the  tape.  Th.s  situation  arises  eitner 
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when  all  of  the  required  map  will  fit  within  the  map  array  at  one  time  or  when  a 
previously  sorted  input  data  tape  is  being  processed.  In  the  first  case  each  input 
coordinate  point  must  be  checked  to  see  if  its  associated  cloud  subdivision  affects 
the  map.  If  not,  the  impact  point  is  not  interpreted  into  the  map.  In  the  second 
case,  when  a  prior  sorting  operation  guarantees  that  only  necessary  particle  de¬ 
scriptions  are  on  the  input  tape,  this  detailed  checking  procedure  is  bypassed. 

Subroutine  SHIFT  (FC-11) 

The  basic  function  of  this  subroutine  is  to  collect  into  the  top  of  the  particle 
arrays  NS,  particle  descriptions  of  the  selected  class  (those  having  their  KTR 
variable  equal  to  JOUT)  and  then  write  them  onto  the  appropriate  tape  (KTAPE). 
Adjustments  are  made  to  the  count  of  the  number  of  members  in  the  selected  class, 

to  NE,  to  the  count  of  available  lines  in  the  particle  arrays,  and  to  the  KTR  vari- 

\ 

ables  of  those  particle  descriptions  which  are  read  out  onto  tape  either  for  the  first 
or  second  time,  in  accordance  with  the  need  to  duplicate  descriptions  of  particles 
that  fall  into  buffer  zones. 

The  actual  sort  procedure  is  one  which  usually  moves  only  the  minimum  amount 
of  information  and,  at  worst,  a  temporary  description  storage  may  be  used  once 
(for  one  description)  during  each  execution  of  SHIFT.  In  this  procedure  JT  is  an 
index  that  proceeds  from  the  top  of  the  arrays  (JT  =  1)  toward  the  bottom,  and  JB 
proceeds  from  the  bottom  toward  the  top.  The  program  first  examines  the  top  line 
to  determine  if  it  is  a  member  of  the  class  to  be  dumped  (KTR(JT)  =  JOUT).  If  it 
is,  it  is  left  in  place  with  its  KTR  set  to  zero,  and  then  the  top  mdex  JT  is  incre¬ 
mented  so  that  the  next  line  can  be  considered,  if  the  top  line  is  empty  (KTR(JT) 

=  0),  the  program  goes  to  the  bottom  of  the  array,  using  index  JB  to  try  to  find  a 
particle  which  is  to  be  dumped  and  thus  can  be  moved  into  the  empty  line,  JT.  If, 
on  the  other  hand,  the  top  line  contains  a  particle  description  that  is  not  to  be  dumped, 
the  program  goes  lu  the  bottom  of  the  array  to  seek  an  available  line  to  which  the 
top  particle  can  be  moved.  In  this  general  way  the  program  proceeds  b,  skipping 
particles  to  be  dumped  when  they  are  found  in  the  top,  moving  particles  to  be  dumped 
to  the  top  when  they  are  found  in  the  bottom,  and  moving  particles  that  are  not  to  be 
dumped  from  the  top  to  the  bottom.  The  sorting  stops  when  NS  particles  of  the  class 
to  ow  d-’vapeu  •.'•ve  been  collected  into  contiguous  cells  in  the  top  of  the  particle  arrays. 
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Next,  the  counter  NS  followed  by  NS  particle  descriptions  are  written  onto  the 
memory  tape  identified  by  KTAPE.  NE,  the  counter  of  empty  (available)  lines  in 
the  particle  arrays,  is  Increased  by  the  number  of  lines  that  have  just  been  made 
available  and  appropriate  adjustments  are  made  to  the  affected  members  of  the 
NP(  )  class  counter  array  .  Then  the  program  returns. 

Subroutine  SLIDE  (FC-12) 

It  is  the  purpose  of  subroutine  SLIDE  to  move  the  (incomplete)  map  ordinates 
which  e  ?*t  In  the  right  (eastern)  buffer  zone  to  the  left  buffer  zone  and  then  to  blank 
out  the  mldsectidn  and  right  buffer  zone  parts  of  the  map  array.  The  common 
argument  variables  NBZX2,  NOX,  NXMAP,  and  NY MAP  are  used  to  communicate 
the  layout  of  the  map  array's  buffer  zones  and  midsection  to  subroutine  SLIDE. 

Subroutine  ZERO  (FC-13) 

This  subroutine  merely  scans  down  from  the  top  of  the  particles  array  and 
collects  NIJ  empty  particle  descriptions  in  the  top  NIJ  array  positions.  Empty  par¬ 
ticle  descriptions  are  denoted  by  a  zero  in  the  associated  position  (same  index)  of 
the  KTR  array. 

Subroutine  DIFUZ1  (no  flow  chart) 

This  subroutine  accounts  for  the  effects  of  atmospheric  diffusion  through  the 
use  of  an  elementary  model  that  adjusts  the  sizes  of  cloud  subdivisions  on  the  basis 
of  only  the  amount  of  time  that  the  subdivision  spends  in  flight.  See  Volume  IV  for 
a  discussion  of  the  basis  of  this  program. 
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FC-12.  Subroutine  SLIDE 


USER  INFORMATION 


Card  Inputs 

/Like  each  of  the  major  subdivisions  of  DELFIC,  the  Output  Processor  requires 
a  deck  of  input  cards  to  identify  and  control  its  operation.  In  brief,  these  inputs 
consist  of,  first,  a  single  set  of  identification  and  overall  control  data  which  will 
be  referred  to  as  "initi  d.zation  data"  and,  second,  a  series  of  "local"  data  sets 
which  indicate  the  geog  phic  limits  of  the  map  to  be  produced  by  the  program. 

Within  each  "local"  data  set  there  may  be  any  number  of  individual  map  requests 
of  the  form: 

NREO  T1  T2  MASCHN  , 

where  NREO  is  a  code  integer  denoting  the  kind  of  computation  to  be  performed,  T1 
and  T2  are  time  arguments  to  be  used  in  deposited  mass,  exposure  rate,  and  ac¬ 
cumulated  exposure  computations,  and  MASCHN  is  a  mass  chain  number  that  is 
specified  only  whei  the  output  is  to  be  in  terms  of  that  single  mass  chain.  Each 
map  request  re3’ilts  in  the  preparation  and  printing  of  a  separate  map  showing  the 
output  indicated  by  the  request  code  NREQ.  The  end  of  a  "local"  data  set  is  indi¬ 
cated  to  the  Output  Processor  by  a  map  request  card  having  a  zero  in  the  NREy 
code  field  (a  blank  card  will  suffice). 

Table  2  gives  details  of  the  card  formats,  program  variable  names,  and  the 
meanings  and  uses  of  variables  for  initialization  data  and  local  'lata  sets.  Not'  that 
a  final  blank  card  is  required  to  indicate  the  e  id  of  the  data  deck  and  cause  the 
program  to  terminate  correctly, also,  that  certain  eard  inputs  to  the  particle  activity 
module  (subroutine  FA  Ml)  may  be  called  for  following  the  fourth  input  cara  for  the 
Output  Processor. 

Card  1.  Output  Processor  Run  Identification.  This  card  can  be  used  to  uniquely 
identify  the  current  run  of  the  Output  Processor.  The  content  of  ti  is  card  is  made 
part  of  the  hard  copy  output  produced  by  the  run. 

Card  2.  Available  Tape  Identifications.  When  required,  extensive  use  is  made 
of  secondary  tape  memory  30  that  tha  Output  Processor  can  be  reasonably  efficient 
in  producing  large  output  maps  or  tabulations.  This  card  shouid  contain  the  "logical" 


TABLE  2 


DETAILS  OF  THE  OUTPUT  PROCESSOR  DATA  DECK 


Name  of 

Data  Set 

Card  No. 

Content 

i 

Program 
Variable  | 

Format 

r  i 

72  character  identifier  for  the 
Output  Processor  run. 

OPID(J) 

(12A6) 

Initialization 
Data  \ 

2 

List  of  logical  tape  numbers  of 
those  tape  units  available  for 
use  in  sorting 

IOT(J) 

ff8I4) 

3 

Overall  control  variables 

IC(J) 

(1814) 

Particle 

Activity 

Data  Set  ' 

(see  Voi.V) 

L 4 

r  --- 

L  ... 

Printer  description 

IH,  IV 

(214) 

. 

r 

5 

Map  parameters  limiting  co¬ 
ordinates  and  grid  intervals 

XMAX, 

XMIN, 

YMAX, 

YMIN, 

DGX  DGY, 
GRUFF 

(7F10.3) 

First  J 

6 

L^cal  control  variables 

JC(J) 

(1814) 

"Local"  1 

Data  Set 

7 

First  processing  request  on 
current  map 

NREQ.Tl, 
T2,  MAS  CRN 

(14,  2F10. 3, 14) 

8 

Second  request 

V 

Request  termination  card  (blank) 

blank 

r 

Next  map  specification 

See  Card  No.  5 

Next  j 

"Local"  i 

Data  Set 

j 

i 

i 

Next  local  control  variables 

See  Card  No,  6 

Next  deck  if  processing  requests 
“ee  Card  No.  7 

i 

c 

Request  termination  card(hiank) 

Final 

Data  Card 

Data  deck  termination  card 
(blank) 

blank 
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identification  numbers  of  ill  tapes  which  are  available  for  use  by  the  Output  Pro¬ 
cessor.  The  program  checks  the  input  values  to  exclude  FORTRAN  system  tapes 
and  the  grounded  particles  tape,  and  use  of  the  remaining  ava;1-ible  tapes  is  made 
only  when  required  for  sorting  particle  data.  As  man>  as  18  tape  numbers  may  be 
listed  but  at  least  one  is  required  by  the  program. 

-Card  3.  Overall  Control  Variables.  Provision  has  been  made  fo**  the  specifi¬ 
cation  of  18  unique  overall  control  variables  whose  values  are  stored  in  the  program 
in  the  array  [  IC (J >,  J  =  1, 18] .  At  present  only  two  of  these  variables  have  been 
given  functions  within  the  program,  and  the  others  remain  for  use  in  control  and 
interprogram  communication.  The  functioning  variables  are  as  follows: 

IC (1 7)  Controls  t.e  entrance  to  the  Output  Processor.  IC(17)  >  0  causes 
the  program  to  stop  without  entering  the  Output  Processor 
proper.  This  setting  is  used  if  only  a  printing  of  the  grounded 
particles  tape  is  desired. 

IC (1  7)  =0  causes  a  normal  entrance  to  the  main  body  of  the  Output 
Processor  regardless  of  whether  the  grounded  particles  tape  has 
been  printed. 

IC (1 8)  Controls  the  option  to  print  the  content  of  the  grounded  particles 
tape.  IC (1 3)  >  0  causes  the  grounded  particles  tape  to  be  printed. 
IC(1S)  ~  0  bypasses  the  printing  of  the  grounded  particle  tape. 

r  ard  4.  Printer  Description.  To  simplify  the  production  of  spatially  undistorted 
maps,  the  Output  Processor  needs  constants  which  describe  the  character  spacing  of 
the  off-line  printer  to  be  used.  These  constants  IH  and  IV  give  respectively  the 
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coordinates.  To  allow  flexibility,  the  scaled  spacing  between  grid  points  on  the  out¬ 
put  map  has  been  arranged  to  be  set  by  the  user.  The  variables  DGX  and  DGY  indi¬ 
cate  the  intergrid-point  distances  in  the  east-west  and  north-south  directions  res¬ 
pectively.  It  should  be  noted  that  on  the  printed  map  the  actual  physical  spacing  of 
the  data  points  is  fixed  by  the  printer's  character  and  line  spacing.  Map  printing 
formats  have  been  arranged  to  achieve  the  greatest  reasonable  data  point  density  on 
the  printed  Huge,  and  on  IBM  printers  this  amounts  to  thiee  lines  per  grid  interval 
in  the  vertical  direction,  and  six  characters  per  grid  interval  in  the  horizontal 
direction.  If  the  user  wishes  to  have  a  map  of  some  particular  scale  produced  by 
the  Output  Processor,  he  must  set  the  parameters  DGX  and  DGY  to  account  for  both 
the  character  spacing  of  the  printer  as  well  as  the  interdata-point  character  counts 
used  by  the  program  (3  lines  per  interval  in  the  vertical  direction  and  six  characters 
per  interval  in  the  horizontal  direction).  Obviously,  zero  values  should  never  be 
assigned  to  DGX  or  DGY. 

An  option  exists  within  the  Output  Processor  to  ~ause  it  to  adjust  the  grid  inter¬ 
vals  put  in  by  the  user  so  as  to  yield  an  undistorted  map— a  map  on  which  the  same 
scale  factor  applies  in  all  directions.  If  the  user  has  specified,  via  parameter 
JC(1G),  the  automatic  undistorted  map  option,  the  program  makes  use  of  either 
DGX  or  DGY  as  the  scale  factor  basis,  depending  upon  which  of  these  two  param¬ 
eters  will  yield  the  largest  undistorted  map  (smallest  scale  factor).  The  last  item 
on  this  card  is  a  ground  roughness  factor  (GRUFF)  by  which  the  program  multiplies 
ail  computed  exposures  and  exposure  rates  before  display. 

Card  6.  Local  Control  Variables.  Provision  has  been  made  for  the  sped.  .ion 
of  18  unique  local  control  variables  whose  values  are  stored  in  the  program  in  the 
array  [JC(J),  J  =  l,  18],  At  present,  only  four  of  these  variables  have  been  given 
functions  within  the  program,  and  the  others  remain  for  use  in  control  and  inter¬ 
program  communication  at  the  local  level.  The  functioning  variables  are  as  follows: 

JC(1)  Output  format  control  variable  ,JC(1)  =  1  results  in  the  printing  of  the 
output  map  in  a  two-line  E  format  which  has  the  power  of  ten  printed 
on  one  line  and  the  associated  multiplier  printed  immediately  below 
it  (see  p.  7  ). 
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JC(i)  -  2  results  in  the  printing  of  a  two-lino  Fll.3  format  which 
has  the  six  highest  order  characters  printed  on  the  first  line  and 
the  five  lowest  order  characters  on  the  second  line  (see  p.  7  ). 

JC(i)  =  3  causes  the  Output  Processor  to  write  a  map  image  onto 
the  multiple  burst  tape  (the  unit  identified  in  parameter  MBTAPE 
as  logical  11).  This  tape  is  written  in  a  format  acceptable  to  the 
separate  multiple  burst  tape  processing  program  MULTIB  (see 
Volume  VII  of  this  documentation).  When  using  the  multiple  burst 
option,  care  should  be  taken  to  see  that  the  tape  unit  identified  by 
parameter  MBTAPE  is  not  also  specified  as  being  available  for  use 
during  the  sorting  operations  of  the  Output  Processor  (Card  2). 

JC(15)  Diffusion  control  parameter  JC(15)  >  0  brings  about  the  use  of  the 
diffusion  subroutine  DIFUZ1.  JC(15)  =  0  bypasses  the  diffusion 
model. 

JC(16)  Automatic  undistorted  map  parameter  JC(16)  =  0  results  in  the 

automatic  adjustment  of  the  grid  interval  DGX  or  DGY  to  yield  an 
undistorted  output  map. 

JC  (16)  ^  0  results  in  no  adjustment  to  the  grid  intervals. 

JC(18)  Grid  interval  adjustment  control  parameter 

JC (18)  =0  indicates  the  user's  permission  for  the  program  to  make 
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efficiency.  This  adjustment  may  result  in  4ncreased  map  resolution 
but  cannot  result  in  decreased  resolution. 

JC(18)  >  0  indicates  the  user's  wish  to  have  no  adjustment  made  to 
the  grid  intervals.  JC(18)  >  0  overrides  JC(16)  =  0,  i.  e. ,  for  an 
aut>matic  undistorted  map,  JC(18)  and  JC(16)  must  both  equal  zerc. 

Cards  7,  8  ,  ■ ,  Processing  Requests.  Table  3  presents  the  meanings  of  the 
computation  codes  NhEQ  and  arguments  T1  and  T2  for  currently  available  computa¬ 
tion  options.  MASCHN  is  the  mass  chain  number  if  the  output  is  to  be  for  a  single 
mass  chain.  Otherwise,  Us  field  may  be  left  blank. 


TABLE  3 

AVAILABLE  COMPUTATION  CODES 


Computation  Code 

NREQ 

Computation  Type 

Description 

0 

Termination  of  Hie  set  of  requests 

1 

Count  of  wafers  covering  each  output  point 

2 

Exposure  rate  normalized  to  time  H  +  1  hour 

3 

Exposure  rate  at  time  H  +  T1  hours 

4 

Integrated  exposure,  H  +  T1  to  »  accounting 
for  time  of  arrival 

3 

Integrated  exposure,  H  +  T1  to  H  +  T2  account¬ 
ing  for  time  of  arrival 

6 

Total  mass  deposited 

i 

Total  mass  deposited  from  time  H  +  T1  to  H  +  T2 

8 

Integrated  exposure,  H  +  T1  to  H  +  T2  assuming 
all  particles  have  arrived  by  H  +  T1  hours 

9 

Same  as  8  integrated  to  infinity 

10 

Concentration  of  an  individual  mass  chain 
(curies/m2) 

11 

Time  of  arrival  of  first  fallout  particle 

12 

Time  of  deposit  of  last  fallout  particle 

13 

Smallest  particle  size  deposited 

14 

1 

Largest  particle  size  deposited 

15 

Mass  from  particles  in  size  range  T1  to  T2 

16 

H  +  1  hour  "normalized”  exposure  rate  result¬ 
ing  from  particles  in  size  range  T1  to  T2  microns 

Binary  Input 

The  Output  Processor  Module  takes  binary  input  from  the  grounded  particles 
tape  (IPOUT)  that  is  produced  by  the  Transport  Module.  The  structure  and  contents 
of  this  tape  are  described  in  detail  in  Table  4. 
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TABLE  4 

THE  GROUNDED  PARTICLES  TAPE,  IPOUT 
(Binary  input  to  the  Output  Processor  Module) 


Logical 
Record  No. 

Record  Content 

Variable  Names 

... 

1 

Identification  word  (IPOUT) 

IPOUT 

2 

Fission  yield,  mass  of  soil  lifted, 
solidification  temperature,  time  of 
solidification,  spare,  time  at  which 
transport  was  terminated,  width  of 
cloud  subdivisions  at  time  of  defini¬ 
tion.  density  of  fallout  particles,  X, 

Y,  and  time  coordinates  of  ground 
zero. 

FW,  SSAM,  SLDTMP, 
TMSD,  SIGMA,  TW, 

HOB,  NCL.T LIMIT, 
BZ,ROPART,XGZ, 
YGZ.TGZ 

3 

Run  identifiers  for  Initial  Conditions, 
Cloud  Rise,  Cloud  Rise -Transport 
Interface,  Transport,  anu  Wind  Field 

(DETID(J),  J=I.  12), 
(CRID(J),  J=l,  12), 
(PSEID(J),  J=l,  12), 
(TID(J),  J=l,  12), 

(WID(J),  J=l,  12) 

4 

Number  of  particle  size  ranges 

NPS 

5 

i 

1 

Central  particle  size,  associated 
mass,  maximum  particle  size,  and 
surface -to-volume  ratio  for  each 
size  range 

PS(T),FMASS(J), 
PACT(J),  SV(J), 

J=l,  NPS 

6 

Topography  identifier 

TOPID(J),  J=l,  12 

7 

Number  of  particle  descriptions  in 
the  following  data  block 

N 

8 

X  coordinate,  Y  coordinate,  time, 
particle  size,  and  mass  associated 
with  6s„n  of  N  particles 

NP(J),  YP(J).  TP(J), 

1  PS(.D  FUTASUU 

J=l,  N 

9 

Same  as  record  7 

10 

Same  as  record  8 

• 

Pairs  of  records  like  7  and  8  are 
repeated  until  all  grounded  particles 
are  recorded 

Last  record 

The  end  of  the  grounded  particles 
data  set  is  indicated  by  a  particle 
count  of  zero 
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♦  #.W*e**K##*»*-IH:**#K**«'t*******HHH,»*>,'*****l**»***********#»*******>***> 


CZUN 

3ZLN 

CiJUN 

CZUN 

C2UN 

>U.v 


4  j 


:zun 


]  t- 

lc 

l? 

1  c 
1  ^ 
20 
21 
22 

2  3 
2  4 
25 
20 
27 

2  o 

29 

30 

31 

32 

33 

3  A 
33 

36 

37 
33 
39 
uo 

41 

42 

4  3 
44 
4- 


CZUN 

4  6 

MAX*NP< 1 ) +NP ( 2 ) 

CZUN 

47 

JPUT -2 

CZUN 

48 

IF  (JIM.LT.4)  G?  TP  6 

CZUN 

49 

00  4  J  =  4  .JIN, 2 

CZUN 

50 

I F ( MAX-NP I J ) )  3.4,4 

CZUN 

51 

T 

N‘AX  =  NP  I  J  ) 

CZUN 

52 

J0UT=J 

CZUN 

53 

4 

continue 

CZUN 

54 

1 p ( MAX-NP ( 1 ) -NP 12)15*6*7 

CZUN 

55 

CALL  ERROR  l  PRZGRf’  »-f>  *  I S2UT  ) 

CZUN 

56 

6 

NS  =  V.AX  +  NP  0*1 

CZUN 

57 

GO  T?  o 

CZUN 

58 

7 

NS*NP( JZUT)+NP( J0UT+1 ) 

CZUN 

59 

Avai!ab!o  Codv 


e  JJ  =  J7UT/2 
KTAPf *!?T< JJ) 
RETURN 
END 


CZUN  60 
COUN  61 
CZUN  62 
CZUN  63 


62 


i! 
i! 
II 
i! 


0 


I 

I 

E 

I 

R 


(i 


l! 

I! 


tTge-Tr  rpr.ov  L  1ST  .•'ECX  «M94/2 
SUBROUTINE  CRD0 

C  P.  PLUSSrR  TECHNICAL  ?PfRATIi-NS  RESEARCH  SK  CR^P 

C 

'  2fcF™  67 

C***AETER  ALL  the  °ARTICLr  PARAVETE-’S  >.AVE  DEEN  RL^O  Kv-..  TA»E 
r«*#tP?'JT*  SL'oRO’JTlNr  cR^p  is  c  A  L  L r  ^  AND  '.-.'RITES  O'.T  THOSE 
C***Pa.KT  I  Cl  “  PARAMETERS  THAT  ARE  LEFT  1  N  C-/RC  Jti  ThL  APPR-PR l«Tr 
C#*#TApES.  *  LEAVES  THOSE  PARTICLE"  THAT  FALL  INTO  THE  ZONl 
C**#V/hJCh  I-'  3  EE  TREA  TF"  NEXT  IN  LOSE  AtO  SETS  NIJ«-1  Ah  a 

C«**SIGNAL  T.  ELBSECL  E“.TLv  CALLED  5  L n RO -T  I f>E S  TO  TREAT  THESE 
C***PaRTICLFS  BEFORE  READING  Nf/.  ONEh  FROi-  TAPE. 

✓ 

r  *#****»***#»  {*«4»**fl<Hi*»*#»**i>**4**»«i*’t*1 


e.RL'P 
CROP 
CROP 
CKO  h 
CROP 
CR^P 
CROP 

^  V  H 

C  3rP 
0  R P 
CRT13 
'R-;P 
E  :<  [  '  P 

CROP 


r  ***•*■»«■*»**»*«■***■********  r.  -•*Hf»*»**»»****1 

C 

^  *****#»****•#.-*  !}*#**#**»»#**#****’<'**'»'»»*>' ‘H.##*******#**#*.***..**** 

r  *#*##*•***«■***►»*****■»**•(,;;:“  >  {5t;t»«J»>»f***il***»»******‘ 

C 

<* J 1 N-2 

C**«  IS  LAST  ZONE  SORTED...  IF  YES,  wRITE  OUT  ALL  40NES  * 

C*#* i r  \eT»  TREAT  LAST  0?NE  DIrFrR"NT lY , 

I  *"  t  LAST  )  EC, ED, 51 

e:  k»k+2 

51  IFtK.LT.21  GO  TO  52 
D?  1  J  *  2  ♦  K  *  2 
jeuT*j 
jJ*J/2 

KTAPEM0T1 JJ) 

1  F ( J-2 131*2*31 


i  * 
14 


CO*- 

"ON  / S E T 1  / 

oROP 

L  J 

1 

D  I A  V 

»r'rTlD(l?i  .IRIS" 

.  I  rXEC  . 

IS  IN  .  I? OUT 

.  CROP 

1  6 

? 

*  SPA  3  *  PiSA-** 

t  T  L  » 

TMP1  .  TYP2 

,  cr:'P 

17 

T  *>\* 

,  u  ,  vpr 

•  * 

H2JPST  .  CLCHr 

•  CROP 

L  5 

U 

7  I  D  t  4' 

R,.  ,  1-tliTR 

.  o  PARI  . 

y  Z  T  A  P  Z  •  Fee  ■“* 

.  CROP 

19 

CP  J"  A 

.  S P A R  5  *  SPA <6 

.  SPAR7  . 

SPARS  *  SPAR9 

CROP 

20 

CO1 

••ON  /SET’/ 

CROP 

21 

A  7 

,<=Z2 

»3ZZ 

» BZ22 

CROP 

22 

» D  E  L  T  A  X 

ir-GX 

.  DGY 

OIFC1-N 

CROP 

23 

•> 

>  D ! !’  A  ~  J 

.rMAFIEC.') 

,  F MASS  (  200 

*  ICt 18  > 

CROP 

24 

L 

♦  i  C  0  K 

,  I  C  'i  i 

*  I  r4. 

*  1  0  T  i  1  3  1 

CROP 

25 

.IP 

.  T  (■•  OUT 

.  ITT ( ; 3  ) 

.IV 

CROP 

26 

c. 

•  jctir 

.JIN 

♦  J^UT 

*  J  P  J  'J  T 

CROP 

27 

? 

.XT'5  (  eP") 

.^TAPf 

.LAST 

, MAP RUN 

22 

r 

.'■•A  >®AY 

.•'IN 

.'■'XREO 

CROP 

29 

r 

.  N  A 

*  Nl  Z  X 

♦  iNZ<ZX2 

CROP 

30 

.NT  7  Y 

♦  NCL 

» NE 

♦  NF 

CROP 

31 

.  N 1  j 

.NMAP 

.  NY  AX 

.  NL'X 

CROP 

32 

.NP!  n  ! 

*  N  R  E  C  • 

*  NS 

.NT  APE  S 

CRIP 

33 

4 

-NT  A  P  ■  T 

.NTASK 

t  NX YAP 

♦  NY, MAP 

CROP 

.‘4 

,  Y  1 N 

.PS ( EDC i 

.PSIZEI2CC 

)  .PAST, 200) 

CROP 

3  5 

4 

,  REPORT 

» ? V ( 20C  1 

*  T (500) 

♦  T 1 

CROP 

36 

•) 

.  T? 

« T  L I V I  T 

♦XI50C) 

.XF 

CROP 

37 

c 

»>'C 

.XV  AX 

.  XVI N 

♦  XNV.AP 

CROP 

h8 

Q 

-XI 

.X2 

♦  X2 

♦  X4 

CROP 

j? 

1 

,v(R«n 

.  VF 

.YO 

.  Y’lAX 

CROP 

40 

COMH-'N  /EFT4/ 

d"ad  1 4 cc° i 

CROP 

41 

CROP 

;  itfrftvk^^pp 

CROP 
CROP 
CROP 
CROP 
CROP 
CROP 
CROP 
CROP 
CROP 
CROP 
CROP 
CROP 


AD 
44 
4  5 

46 

47 

48 

4  9 

5  C 

51 

52 

53 

54 

55 

56 

57 

6  6 
59 


63 


31 

r<S*N?!JHN?(  J+l) 

CROP 

60 

IFlNS)  2*2-3 

CROP 

61 

3 

CALL  5HlcT 

CRDF 

62 

2 

L»0 

CROP 

63 

WR»Tf  (HTAPE)L 

CROP 

64 

1 

RFVIND  <Ta°E 

CRDP 

6  '• 

c*#* 

IS  LAST  Z2NF  WRITTEN  3UT... 

CRDP 

66 

52 

I F ! LAST )  8.10.8 

CRDP 

67 

C*** 

can  we  dusp  all  particles  New... 

CRDP 

68 

a 

IF ( NO ( J I N )~NMAX )  5. 5.6 

CRDP 

69 

6 

N  ,%,-f-sX 

CRDP 

70 

NPf  JIN)  *"iP  (JIN)  -NS 

CRDP 

71 

c-e  TO  7 

CRDP 

■?*? 

R 

NS^NP (JIN) 

CRDP 

73 

7 

J2UT  *J I N 

CRDP 

74 

KTA°E*I?T ("INi 

CROP 

75 

*F'NS)9.9.71 

CRDP 

76 

71 

CALL  SHIFT 

CRDP 

77 

C***ARE  AN*  PA°.TICLc5  LEFt... 

CRDP 

78 

IF(NPIjrN))  9*Vie 

CRDP 

79 

9 

L*0 

CRDP 

80 

WRITE  (KTAPE)L 

CRDP 

81 

REWIND  KTaPE 

CRDP 

82 

K 

NIJ«-1 

CRDP 

S3 

RETURN 

CRDP 

84 

END 

CRDP 

85 

«4 


SIRFTC  0IFUX1  L  1ST  *0?CK*M94/2  DIFU 

SUBROUTINE  DIFUZ1 (RZ2)  DIFU 

C  T.W.  SCHWENKE  TECHNICAL  OPERAT 1 55NS  RESEARCH  DIFU 

C  11  3CT  66  DIFU 

C  DIFU 

C  THIS  SUBROUTINE  cXPANOS  THE  CLOUD  SUBDIVISIONS  AS  A  MEANS  OF  DICU 

C  APPROXIMATING  THE  EFFECTS  OF  ATMOSPHERIC  DIFFUSION.  DIFu 

c  D 1  F  u 

£  ##*«##•*###**##»#*■#*###***##•»##«■  #**.****#***  *-»***#***-»****lf**  w***r.  j  ru 

c  D 1  Fu 


COMMON  /  S  E  T 1  / 


DIAM  *DET [ D ( 1 0 ) ♦  I R I SE 
S"  ♦  SPAR  ♦  SSA v 

T?M  .  U  ♦  VPR 

TIDt&D) ♦  Rf'IN  .  IDISTR 
SPAR4  *  SUORAD  *  RAD^’AX 


I  EXEC 
TME 

W 

E  PARI 
XGZ 


1SIN 

TVRl 

HB.JRST 

MBTAPE 

YGZ 


isout 

TMP2 

-CLDHE 

fsum 

TGZ 


##*####•»##*#«##**#*'»*****'***■*<'* *«*■#■*******#**#**#*•»*•* in.**# j  rU 

DIFU 

COMMON  /SET3/  D !  FU 


1 

ez 

.TZ2 

.  BZZ 

♦  3Z2  2 

DIFU 

20 

. DELTA/. 

»DGX 

.  DGY 

♦DIFCON 

DIFU 

21 

*2 

»D I FADJ 

. CMAS ( ECC  ) 

» FMASS ( 2CC ) 

.JC1 12  ) 

DIF. 

22 

4 

.ICON 

.  I  CTR 

.  IH 

.IDT (18) 

DIFu 

23 

e 

•  ID 

. | POUT 

»  1 T  T  ( 1  2  ) 

» 1  •/ 

DIFU 

24 

6 

.  JC ( IF) 

.JIN 

.  J2UT 

. JPOUT 

DIFU 

23 

7 

♦  FTRliO'1! 

.KTAPE 

.LAST 

.MAP RUM 

D I  rU 

25 

P 

.VARRAY 

.min 

.MXPEC 

DIF. 

27 

4 

♦  N 

»N  A 

» NEZX 

.N32X2 

DIFU 

2a 

1 

.  NBZ  Y 

•  NCL 

.NE 

»NF 

DIF. 

2  V 

? 

.NIJ 

»NMAP 

♦  MM  AX 

.NOX 

DIFU 

.0 

.  NPt  21  ) 

.  N  R  E  Q 

.NS 

.MTAPES 

DIFU 

31 

4 

.ntapet 

.NTASK 

t.NXMAP 

.NYMAP 

DIFU 

32 

p 

.  YV  T  M 

.PS  (SCO 

.P:IZ£(200) 

♦PACT ( 2CC ) 

DIFU 

33 

6 

.ROPART 

.SVI2CP) 

.TI5CCI 

♦  T1 

DIFU 

34 

7 

*T? 

.TLJMIT 

. X  1 500 ) 

♦  XF 

DIFU 

j  t. 

r 

.XO 

.XMAX 

.  X  M I  M 

♦ XNMAP 

DIFU 

36 

o 

.XI 

»X2 

» X3 

♦  X4 

DIFU 

37 

A 

i 

.Y(?ro i 

♦  ”r 

,Y0 

♦  YMAX 

DIFU 

32 

'•#•*•*#*####*•11**##*  #  # ! 


COMMON  /SETA/  OMAP(ACCC)  DIFU 

DIFu 

###»##*#####*»  *#»#**#**#*■**#*******#?*■»*****»*#***■*•******» ’■■*if*-»»if’fifif»Kp)j[ru 

DIFU 

THE  VARIABLE  CIFCON  HOLDS  AT  THIS  PD!,\T  THE  PRODUCT  2.C*D1FCUN  DIF„ 
WHERE  D I F CON  ORIGINALLY  WAS  PUT  IN  Au  AN  ATMOSPHERIC  DlFOoION  i.1F« 

CONSTANT.  DIFU 

“0*2  has  the  square  OF  BZ/2.C  .THE  KI3IH  OF  THE  CLOUD  SUBDIVISION. DIFu 

DIFU 

l  F ORMAT ( ? 3H  DIF^USlV"  GROWTH  O'7  CLOUD  SUBDIVISIONS  IS  EXCESSIVE,  .3X0 |ru 
1,3HX»  FI  1 .3  *3HY=  F11.3)  DIFU 

DIFU 

##########*  ###*###  ##«*#####*#####**###**#*#*########*#»  *##**#*#**  *##»*~,[p-,_, 

DIFU 

DATA  NG3?C  /O/  DIFU 

DIFU 

RZ2  =  SQRT(''IFCDN*(T(  IP)-TGZI+FZ20!  DICU 


65 


ITfRZ?  .LF.SZ?) 

RETURN 

DIFU 

60 

RZ2*C'Z? 

DIFU 

6) 

NC.ee  *»NC,erc+i 

DIFU 

62 

IFiNGSeF  .CiT.n 

•AND. NG2CF  .LT.5C)  WRITE! 1S2JT.1 )X( IP) iYIIP) 

DIFU 

63 

RCTURN 

DIFU 

64 

END 

DIFU 

65 

66 


S I B PTC  LETSC  LIS"  >p'ECK»M94/2 

SU3R0UT !  NC  LETSG? 

C  26  ~Z  P  67 

c  P.  EL’JSSER  TECHNICAL  ZPERmTU'NS 


RESEARCH  SR  LETSG2 


C 

c##* 

c*## 

c### 


THIS  SU3R0UT  INE  REAPS  THE  PARTICLE  PARAMETERS  FROM  TAPE  IPLST 
AND  CALLS  THE  APPROPRIATE  SUER2JTIMES  TO  PROCESS  THtM. 
SUBROUTINES  CALLED  Pv  LFTSGC, ... 

C###PR0C .ZERO. COUNT. SHI  FT. CROP. 

NUMBER  n c  E’-PTV  SPACES. 

CG#*  WARRAY=  OU’1:  NSI ON  SF  PARTICLE  ARRAY 
c***  NP(I)*  MUMPER  LF  PARTICLES  FALLING  IN  ITH  Z0NE 
C*««  KTR(I)=  ZZNc  INDENT I F I  CAT  1 2N  NUM3FR  of  ITH  PARTICLE 
C*tf#  MIJ=  NO.  OF  PARTICLES  IN  NECT  DATA  BLOCK • 


#*#*#***«< 


COMMON 


1 

2 

•x 

4 


C0MM3N 


/SET  1/ 
DIAM 
SB 
T2M 

T I  h ( 40 
SPAR4 
/SFT  3/ 


.0ETIDI12) .’RISE  » 


SPAR 

U 

RMIN 

spar; 


SSAM 
VPR 

id; str 
SPAP6 


R2 

.EElTAX 
.01  PAD J 

•  ICON 
.IP 

.JC(IS) 
tlCTR  (  son ) 

.MAPRAY 
.N 

.N3ZY 
.NIJ 
.NPI21  ) 

»N1 APE T 

•  Yf'IN 
•RCPART 
.T2 
.XC 
.XI 

. Y ( 5  DC  ) 
C2VM2N  / S E T 4 / 


»PZ  2 

.OGX 

. FMAS ( SCO ) 

♦  1 C  TR 
.IP2UT 
.JIN 
»KT  APE 
.MIN 
.NA 
.NCL 
.NMAP 
.nrfo 
.NTASK 
.PSI5C0) 

. 5 V ( 200 ) 

♦  T  L  I M I T 
.  XMAX 
.X2 

»  VF 

CY AP ( 4C00 ) 


OATA  PR2GRM/6HLFTSG0/ 

*##*«*****#*»**»  i 
************************  5 { 

NE “MAPRAY 
K“J I N+l 
03  2  I»1 «< 

?  NP( I ) *0 

03  3  I-l.MARRAY 

3  <TR( Il-O 
ASSIGN  KO  T?  NSN 

4  (WAD  (I  POUT)  NIJ 


IEXEC 

TMS 

W 

SPAR1 

SPAR7 


I  SIN 
TMP 1 
HBJRST 
MBTAPE 
SPARS 


I  SOUT 
TMP2 
SCLDHU 
F  SUM 
SPAR9 


•  B2Z 
.  DSY 

»F MASS ( 20C ) 
.  IH 

.  ITT ( 18  ) 

.  JOuT 

.  LAST 

.MXREO 

»N5ZX 

.  NE 

»NMAX 

.NS 

.NXMAP 
.PSIZEI200) 
.  T  (500) 
.X'SCO) 

»  Xfl  I N 

*  X  3 
.Y2 


,BZ22 

♦  D  ’  FC'JN 

.  I C 1 1 8  ) 

.  1 3T ( 1 8 ) 

.  IV 

» JPO JT 
.MAPRUN 

.N32X2 
.Nr 
»NZX 
.NTAPES 
.NY MAP 
.PACTI 200) 
»Tl 
.  XF 

*  XNMAP 
.  X4 

.  YM.AX 


LETS 

LETS 

LETS 

LETS 

LETS 

LETS 

LETS 

LETS 

LETS 

LETS 

LETS 

LETS 

LETS 

LETS 

LETS 


LETS 
LETS 
LETS 
LETS 
LETS 
LETS 
LETS 
LETS 
LETS 
LETS 
LETS 
LETS 
LETS 
LETi 
LETS 
LETS 
LETS 
LETS 
LETS 
LETS 
LETS 
LETS 
LETS 
LETS 
LETS 
LETS 
LETS 
LETS 
LETS 

*<t  #*###  »*l_r  Ti 
LETS 
LETS 
LETS 

ETS 
LETS 
LETS 
LETS 
LETC 
LETS 
LETS 
LETS 
LETS 
LETi 


0 

1 

2 

3 

4 

5 

6 
7 

3 
9 

10 

11 

12 

13 

14 
lc 
16 

17 

18 

19 

20 
21 
22 

23 

24 
26 
26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

4  c 

46 

47 

48 
40 
‘0 
51 
62 
v3 

54 

55 

56 

57 

58 
39 


67 


C***  ARP  v'p  O0NF...  P«YES.  9“N0. 

LETS 

60 

IF(N’J)  5.8.9 

LETS 

61 

5  IRR0R--5 

LETS 

62 

7734  CALL  ERROR  ( PROGKM  .  i  RROR .  I S3UT  ) 

LETS 

63 

C»**DUMP  AND  RETURN 

LETS 

64 

8  CALL  CRDP 

LETS 

65 

GO  Te  18 

LETS 

66 

9  I F ( N1 J-MARRAY )  11.11.10 

LETS 

67 

1C  IRR0R»-1 0 

LETS 

69 

G 0  TO  7734 

LETS 

69 

C###  js  THCRr  ENOUGH  RO0M  TO  READ  NTX 

T  DATA  BLOCK... 

LETS 

VC 

C**»  12=  YFS.  PROCEED.  SHIFTING  ZFR0S 

IF  NECESSARY. 

LETS 

71 

11  I F ( NI J-NE  )  12.12.15 

LETS 

72 

12  GO  TO  NSN • ( 100 .20C ) 

LETS 

73 

IOC  ASSIGN  200  TO  NSN 

LETS 

74 

C*##  13*  READ  PARTICLE  PARAMETERS  DEC 

RE. ME NT  EMPTY  jPACES.  PROCEjo 

LETS 

75 

C***PaRTICLES  READ  IN  AND  READ  NUMBER 

OF  PARTICLES  IN  NEXT 

LETS 

76 

C**#  nATA  BLOCK. 

LETS 

77 

1 ?  RFAO  ( I POUT 1IX(I).YII).T(I).PC(I) »FMA3 (I).I  =  1»NIJ) 

LETS 

73 

NE=NF-NIJ 

LETS 

79 

CALL  °ROC 

LETS 

80 

CC  T?  4 

LETS 

81 

2^0  CALL  ZERO 

LETS 

82 

GO  TO  13 

LETS 

83 

15  CALL  COUNT 

LETS 

84 

call  shift 

LETS 

85 

GO  TO  11 

LETS 

86 

If  RETURN 

LETS 

87 

END 

LETj 

88 

68 


sipftc  lnag  LiSTxoro. xwa/2 

.'.LfiRaLT  1  NE  LINK  fa 
v-  T  • /« •  oCriX L  >xn. L  Tii^i  iN  i  C/-u_ 


Ti./.4j  •<  »_oi_  AriCn  wuTP.jT  AKi/Cc.oo.uR 


LNKo 
LNRo 
i  lsao 


c 

!;•  FFo  5  7 

L  r-,  o 

3 

c 

F  I  r<  T  HALF  OF  T  ri£ 

.■„TP..T  PKOCEoodR 

LNK.O 

♦4 

c 

1 1 1 1  3  PR*/  uK^''!  1  .>i 

J  T  ImL  i  v  m  » *  L*  k<i*!  1  o  f  i  iiuo 

F  Thl  oL  T  Ps-f  r 

Li*k6 

1/ 

c 

•  ThL 

IT  ifti-Lj  JHt:  r  i  *<^T  ^mf<T  o  F 

t He  part  iC'-L 

II V 1 T  Y 

UixXv  0 

U 

c 

“L/OtLL  (  P  A>‘,  i  )  r 

c 

C  t  L  ;)m  T  m  x^wiilu  JY  f  rl€_ 

zjcC  ZNis  Pa r  T  o f 

THt 

Lilso 

7 

c 

i-’mKTICi-E  activity 

t-.i/ywlL  ixi»ICn  ••  ILL  U  ^  c^ULL 

\j  THL 

LNxe 

o 

c 

FaE-CCTIJN  Zr  L*1 

'Six.  V  • 

i  ■ 

U.NKU 

7 

c 

. 

f 

L.Wno 

lU 

(- 

■Vi  v 

>■  It  i-  x  ir  H  -H  ik  K  -Jr  V  >  V  *■  ->  A-  it  it  ir  #  #  #  4MHf 
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1N0T  BE  GUARANTEED.  TnJj  RwN  *  AS  CoNTJNjED  «.ITh  oSIg  INTERVAwo  um.i 

7x 

23H  ADJUSTED  FOP  MAXIMUM  EFFICIENCY.) 

Li.kS 

7c 

32 

FORMAT! 14. 2F1C.3  *14) 

LNA> 

Tj 

33 

FORMAT  (2>t-.0vNACCLPTAtiLL  Rl^ULST  ...1,4) 

LAN  7 

7h 

34 

FORMAT  (////  / /  I  6X  *1  5hRE*-GE  oT  Ng..3ER  l4///i  jX  *5rlTYPE  14. 

.OXShTi  *  Flu. is/ 

1  j 

1 C . 1  * i CX *6HT2  a  FI G. 1  * lOX * VHMAoCHN  *  14) 

LNk  7 

7o 

35 

FORMAT ( 15  ) 

Li\K7 

77 

C 

LNK7 

7  o 

C  #*##**•»  *«#**«  a*  »*>#*****»****»*•  »**  *  *****  «».*•*** »*«***» »*»*#* 

N  ir  Ar  ■«  rfifc*  *  *  k  ^  4  (j 

(•  ##«**»  ^  vtt -1  *«<.««««  *«  *»  «»**#*»*«  .******** 

it  *k  A*  h  kir*  k  *  H  w ,  >  N  fj 

C 

LNk  7 

G  A 

data  prbgR'!  /6n  links/ 

LNKS 

ft  2 

LOGICAL  JD.KD2S.IG0 

L  A  K  v 

0-> 

c 

LNK  5 

04 

I  GO  =  •  TRL  £ . 

LN/.7 

NLL  =  0 

LNk  v 

OJ 

FswM*C  .  0 

Li,k7 

o  7 

11  V 

IF  <  FSGV « L-i.O.<3  )  GO  TO  .151 

LNk  7 

00 

WRITE! ISOuT. 21  FEU 

LNK.  7 

o  ; 

FStnO.G 

LN.K7 

>0 

c 

L'>K  7 

>1 

Cl  1 5] 

READ  LIMITS  ON  ARE/,  OF  INTEkE-T 

LNk  7 

115] 

READ  (  1 0  I A  •  5  )  Al-'AX  *Xi I'**  Yi»  AX  *  Y;.;jN  *DGA  *DGT  *GRGF  F 

LNK.  5 

IF  1  GRUFF ,"0.0.0  CRUFF  *  1.0 

LNK.7 

>4 

1 F  ( P6X+DGY  1120*120  *121 

LNK  7 

12C 

WRITE  ( ISOuT  *17) 

LNK  7 

REWIND  IPOLT 

L  NK  > 

PRINT  17 

Li.kS 

>0 

STOP 

Li\K7 

>  > 

LNk  7 

1  uO 

£ 

RUG  OTHER  SPECIFIC  INPuT 

LNK  7 

1UA 

121 

READ  (  I S I N*  1  5  )  ( JC  1 J )  *0  =  1  *1  o  ) 

LNk  V 

ntask=ntask+i 

LNK  3 

iw/i 

MQ=C 

LNK  7 

iw4 

C 

L  NK  7 

iv  J 

C 

CHECK  PRINTER  DESCRIPTION^ 

LNK7 

ICO 

NI  =  I H* I V 

LNk  V 

1J7 

I  F  (  N I  )  6.01  .601*122 

LNK  9 

Ig'O 

601 

JC ( 16 ) =1 

LNk  7 

WRITE  (  I S JUT .31) 

LNK  7 

4iO 

C 

LNk  7 

in 

C  122  WHITE  A  LOCAL  rtLADING 

LNk  7 

122 

WRITE  (  IoOUT  *23  )NTA3K 

L..K7 

i  X  0 

a-RITF  (  IuC«T  *24  )  AMIN  » AM AX  » YMI N  *  YMAa  *dgx  *dgy 

LNk,  7 

U4 

write  ( IS0UT  .25 ) t JC! J)  *J  =  1 *18  ) 

LNK7 

i  A  0 

WRITE  US0UT.4)  GRUFF 

LNK7 

j.  lb 

GO  TO  1205 

LNK7 

1 1 7 

r 

y. 

I.NK  5 

i  Lo 

1211  in JCm.E-J.3)  *RITF  ((■'.bTAPEINUw.NoL 

LNK* 

1 1  > 

76 


n  n  n  nnn 


c 

1209  IF  ( FSUM.EQ.Q.O )  6Z  T0  1219 
WRITE  (IS0UT.2)  FSOM 

C  1219  READ  A  REQUEST  F0R  PROCESSING 
1219  REAO  ( IS IN*32 )NREQ*T 1 *T2»MASCMN 

I F I MA5CHN»EU»0»AND*NREU.NE . 1 0 ) GZ  T0  1210 
IFiXASCHN.GT.71.AND.MASCHN.LT.162i  G«  TU  1210 
CALL  £RR0R(PR0GRM*12O9.IS0UT) 

MASCHN-95 

1210  00  935  I»1 *NMAP 
935  0MAP(I!«O.O 

00  936  1*1 tKARRAV 

X(I)*O.C 

936  vm*oso 
1C0N*0 
MAPRUN*0 
NRG*NRC+1 
FSLM  *  C.O 
NTAPES*NTAPET 

02  937  1*1 *NTAP£T 

937  I 0T ( ! ) *  1  TT (  1 ) 

IS  NRE0  AN  ACCEPTABLE  REQUEST 
NC  T0  1215 

IF(NRE5iiZl2,1212.12*3 
121;  I F (  NREQ-iVXREU )1214*1214»l2i5 
1215  WRITE  ( loOLT  *33)NRE(J 
G2  T0  1209 


LNKV  i20 
LNKV  121 
LNKV  122 
LNKV  123 
LNK9  124 
LNK9  125 
LNK9  126 
LNK9  127 
LNKV  1 2o 
LNKV  129 
LNKV  1 jO 
LNKV  IjI 
LNK.9  U2 
LNK5  1j3 
LNkV  u 4 
LNKV  1j-j 
LNkV  i3o 
LNkV  1-7 
LNKV  130 
LNK9  139 
LNK9  140 
LNKV  141 
LNKV  142 
LNK9  143 
LNKV  144 
LNKV  143 
LNk7  1*6 
LNkv  1*7 


LNkV  l*to 


MX^Ew  I  u  MAX I MuM  NLNcj£<t  k*E  CAoCJLmT  1  *11  ALL2/*oO  F*R  I*.  LmLC* 

121?  A"  MORE  RE-UEoTo*  PRE^ARc.  T <u  R-T JRN  To  Rc«G  L*CAL  i/ATa. 

1212  CONTINUE 
GO  TZ  119 
1214  IPJUT-JPOUT 
RtrtlND  IPSoT 
REAOI  IP0UTITST 
REAO  (IP2UT)(T*T.J*1.16) 

«EAD (  I PoLT  J  iDETlD ( j) *J=1 *12) . (DcTiDio) *o*i*A2 ) *<Dfc.TlO( J)  *J  =  i *12) * 
1 (OETIOI J) *J=1 *12  i  .(OETliVi J) *J  =  1 *i*) 

RLaD<  IP«iLT).»P6 

RE  AO 1 1  POUT  )  (  PS  1 2E  (  ')  *F.-.aSw(  J)  *PAC  T  ( J  )  *oVI  - )  *-.=  i*NPs} 

REaDI IP2-T) (OtTlOIJ) *J*1  .12) 

I F ( NRG-1 11221 *.222*122* 

1221  WRITE  t IS30T  .3  ) 

122?  WRITE  (  I  jZw  T  *  34 )  NRw' *  NRLi*  *  T 1  *  T2  * -lAiCHN 
RUr S«N  =  SSAM/GRLFF 
JG2  =  1 
I TAB=NP3 

FiSNUM*rW»l  .45E13 


LNKV  149 
LNKV  130 
LNkV  lJi 
LNKV  .  132 
LNKV  ivj 
LNKV  1-* 
LNk7  L;j 
LNKV  lvo 
‘-'.KV  1v7 
LNkV  l-o 

LNK7  I-* 
LNkv  ioO 
LNKV  16i 
LNKV  lo2 
LNKV  lo3 
LNkV  164 
LNK>  ib3 
LNKV  loo 
LNKV  1 o7 
LNKv  loo 


1 F  (  NREO •  u  •  1 3  •  OR* NRE <-  >u<  i o  J  o-  T-  122- 
Tl*TIKSEC(Ti*3*0) 

T2*TJ,VSEC(  T2.3.0) 

1223  TlXt*Tl-T0- 
TtNTER=TIXE 
TEXIT*T2-TG2 
JO*. TRUE. 
kD0S  =  . FALSE . 


LNkv  i o  7 
LNkV  170 
LNKV  i  1 1 
cNKV  172 
LNKV  1 7 j 
LNKv  i 7* 
LNKV  ! 7 v 
LNK V  i 7o 


GZ  T0  (  8- * 8 1  *7 V  •  78  *8 2*-6Q.*6v  * 8 2  *  76  *  63  » 30  *60  *30  *80  *80  *81  * 60  *20  *60  *dJLi\K'v  177 
1  )  *NRE»1  LNKV  1  7o 

81  T 1ME*3600 .0  LNKV  i7V 


77 


60  T0  79 

LNKV 

leO 

•2 

KD»6a.TRUE. 

LNKV 

loi 

G0  T0  78 

LNKV 

io2 

63 

JG0-2 

LNKV 

103 

6 1 SNUMaF 1  $NUM*1 .  E*4 

LNKV 

104 

RUFSAK-SSAM 

LNKV 

1 05 

G0  T0  79 

LNKV 

10b 

78 

JD». FALSE. 

LNKV 

107 

F 1 $NUM«  f j  SNUM*  3600 • 

LNKV 

100 

79 

CALL  PAM2 

INK* 

lev 

80 

continue 

lnkv 

1  vO 

C 

LNK» 

1  »* 

C 

IN1T1ALUE  F0N  PN0CES6ING 

LNKv 

1  v2 

ICTRaO 

LNk» 

id 

622aB2#0.3 

LNKV 

1V4 

c)Z22«BZ2*BZ2 

LNKV 

1/3 

C 

LNkV 

1  til 

C 

IS  USE  0F  DlFFubUN  MBUEl  REwLEaTtlv.  YE-  T-  121<1 

LNKV 

iv7 

IF ( JC ( 1 3 1 11217.1217**216 

LNKV 

1  VO 

c 

LNKV 

1  V  V 

0216 

ADJUST  BUFFER  LONE  «1utH  T-  AL**.'  F-n  F-R  1  HC-rili.G 

UlFF-aWL  uHu..ThLi«KV 

2lO 

C 

OF  CL3UD  SObDIVISif No. 

LNK» 

20  * 

1216 

ftZ2"BZ2*DIFADJ 

LNKV 

202 

C 

LNsv 

2u  j 

0217 

16  A  GRID  INTERVAL  A  J-wSP.lNT  ALA.-il  TTtD  Yt6  T« 

1  *30  lNk» 

204 

1217 

lFUCUai-1113u.131.13U 

ui.kV 

2u3 

L 

LivKv 

2ub 

C 

Nu  ADJUSTi-iEi.T  PERMITTED  -PTI..N 

LltKT 

2-7 

131 

NuZX»-l 

LNKV 

2uo 

NbZY«-l 

LNnv 

2uv 

NaZX2«BZ2/06X+*.C 

LNK» 

2*0 

1311 

TST-NoZAZ 

LNkv 

211 

8Z2»T$T*DGX 

LNKV 

2*<t 

C 

TH;  bUFrER  2J.NE  la  N*.  AN  INTEGRAL  .wMoEr  -F  GRID 

INTER VaLo  ..lot  L..s» 

2*3 

61)  T2  140 

LNKV 

2*4 

C 

LNK» 

2*3 

C  1 3t 

16  AN  UNDISPUTED  MAR  DESIRED  YES  T-  1  iOc 

LNKV 

2*6 

1 30 

I F ( JC ( 1 6 ) )la01.1So2.13ul 

LNKV 

<*7 

1  3C? 

DloY»lV 

LNKV 

2*e 

DISXMH 

LNAy 

iii 

KD«2.0»DISY/D1SX 

LNKV 

it  0 

D13X»RD*DGY 

LNKV 

in 

I F l DISX-DGX 113:3.1 304 .1304 

LNKV 

2  22 

1303 

DGX-DISX 

LNKV 

223 

1304 

DUY«DUX/RD 

LNKV 

224 

GO  T0  131 

LNkv 

223 

13C1 

IF(  JCIlS)  .uT.O  1  u6  T  <i  131 

LNKt 

220 

t 

L  NK  > 

a  7 

C 

EFFICIENCY  ADJUSTMENT 

L«K7 

iiO 

No2X>bZ/DGx 

L>.nv 

it  V 

TaT»Na2X 

lnkv 

£3M 

T6T»TST*DGX 

LNKV 

£31 

lF(oZ-TSTH33. 1341. 1361 

LiVKV 

2j2 

1361 

IFINBZXU33. 13o»134 

LNKV 

233 

136 

NoZX— 1 

LNkv 

2a4 

NbZX2«l 

LNKV 

233 

GO  T0  137 

LNK. 

2ab 

133 

13RSR»-1 33 

LNKV 

2a7 

Gt  T0  333 

LNKV 

230 

C 

LNKV 

2)7 

78 


C  114 

ADJUST  DGX  T9  MAKE  62  AN  INTEGRAL  MULTIPLE  UP  IT. 

LNX9 

240 

1)4 

NBZX-NBZX+l 

LNXV 

241 

TST-NBZX 

LNXV 

242 

DQX>BZ/TST 

LNK9 

243 

1)41 

NBZX2‘(NBZX*l)/2 

LNK.9 

244 

C 

LNA9 

243 

C 

New  PBR  THE  Y  D I MENS  I  BN 

LNK9 

246 

1)7 

N8ZY»BZ/DGY 

LNXV 

247 

tst-nbzy 

LNXV 

246 

TST-TSTWDGY 

LNXV 

249 

IPlBZ-TST )  133.1)11.1371 

LNX9 

250 

1)71 

IF (NBZY) 135 *13 8.139 

LNX9 

251 

1)5 

IHRUR— 135 

LNxV 

2>2 

G0  T9  333 

LNxV 

253 

1)6 

NBZY—1 

LNxV 

254 

00  TB  1311 

,NA9 

255 

139 

N3ZY-N8ZY+1 

LNXV 

256 

TST-NBZY 

LNXV 

25  7 

06Y-BZ/TST 

'.NX  9 

25o 

GB  T0  1311 

LNXV 

2sV 

C 

1  A>  *  i  • 

ioO 

14C 

XB»XMIN-3Z2 

LNXy 

26a 

Y0»YMI N-BZ2 

LNXV 

Out 

XF»XMAX+bX2 

LNX> 

205 

YF«YMAX*tlZ2 

LNXy 

oo*. 

C 

LNXV 

2o3 

C 

PREPARE  U  RRUCEoo  .juTP^T 

LNXV 

2ob 

u.XV 

2o7 

NYMAP  «  JYMAX  -  YMINI/OGY 

LNXV 

Zoo 

NXMAP«(XF-X0)/DGX 

LNXV 

2o!; 

NST  ■  NMAP/NYMAP 

LNXV 

270 

IF(NXMaP-NST)  1503. 150j. 1402 

LNX9 

27  i 

140? 

NXMAP»N5T 

LNXy 

27  2 

XMAP«NXMAP 

LNXy 

27j 

ZZ*  t  XF-XZ  )  /  (  XMAP*DGX  ) 

LNX9 

274 

NZ«ZZ 

LNXV 

275 

To! »NZ 

LNK.9 

276 

IF(ZZ-TST)  150v. 1501.1401 

LNXV 

277 

1500 

1RNOR«-130>- 

LNXy 

27o 

GO  T?  333 

LNXV 

27  y 

1501 

NZ«NZ-1 

LNXV 

2o0 

GO  T?  1401 

LNXV 

2oi 

1503 

NZ«*. 

LNXV 

C 

LNXV 

2oj 

C****»**«  EM3  „F  CALCULATION 

LNXV 

2*4 

14C1 

NJX«NXMAP-NBZX2-NBZ  X2 

LNXV 

l  F ( N2X ) 1403 .14^3.1404 

LNXV 

2.6  6 

14C3 

IRR3R  — 1403 

LNXV 

20  7 

333 

CALL  ERROR  l  RROGRi'.  .  I  .IRoR .  I  ^OOl  ) 

LNxV 

2oo 

00  TS  1211 

LNXV 

2oV 

1404 

ex«Nzx 

LNXV 

2yO 

OELTAX«eX*DGX 

LNXV 

2y  1 

C 

LNXV 

hi 

1502 

WRITE-1  I50LT .27I0GX.DGY 

LNXV 

2»3 

C 

LNaV 

2y4 

X1-X0 

LNXV 

2V5 

X2-X1*BZ2 

LNX9 

2V6 

X3*X2*DELTAX 

LNXV 

2y  7 

X4«X3*8Z2 

LNXV 

2  Vb 

300 

IFJNZ-NTaPEM  200.200.201 

LNXV 

2Vy 

79 


200 

MIN-NZ 

LNX9 

300 

LAST-0 

LN<9 

301 

04  Tf  202 

LNK9 

302 

201 

M1N-NTAPES 

LNK9 

303 

LAST-1 

LNK9 

304 

202 

J1N-2*MIN 

LNK9 

303 

IP  I N2 ) 203 1 204  >20} 

LNK9 

304 

20} 

1RR0R— 203 

LNK9 

307 

60  T0  333 

LNK9 

3J0 

204 

CALL  RUN1 

LNK9 

309 

CALL  MAP 

LNK9 

3i0 

60  T0  1211 

LNK9 

311 

20} 

I F  <  LAST  1 204*207 • 20V 

LNXV 

3x2 

204 

IRR0R— 204 

LNK9 

3  *3 

60  T0  333 

LNX9 

314 

207 

CALL  LETSG0 

LNX9 

3i3 

ICTR-l 

LNX9 

3  46 

D0  208  INDEX-1 tMIN 

LNK9 

317 

CALL  MAP 

LNK.9 

3io 

IP0UT-I0T (INDEX) 

LN<9 

319 

X1-X3 

LNXS 

3X0 

X2-X4 

LNK9 

3X1 

X3«X4*DELTAX 

LNN9 

3X2 

X4-X3+BZ2 

LNK.9 

3X3 

CALL  SLIDE 

LNK» 

3X4 

CALL  RUN1 

LNX9 

3X3 

2oe 

REWIND  1P3uT 

LNX9 

3X6 

CALL  MAP 

LNA9 

3X7 

60  T0  1211 

LNK9 

3X6 

2C9 

CALL  LETSG0 

LNK» 

3X9 

K.IN2-IP2UT 

LNk.9 

330 

ICTR-l 

LNKV 

331 

AIN-MIN-1 

LNkV 

3  3X 

D0  210  INDEX-i.XIN 

LNK9 

333 

CALL  MAP 

LNK.* 

334 

lP0oT«I  <4T  ( I NOEX ) 

LNK9 

333 

Xi»X3 

LNK9 

330 

X2-XA 

lnk> 

337 

X3-XA*DcLT**X 

LNK.9 

330 

X4-X3-*dZ2 

LiNk* 

339 

CALL  SLIDE 

LNK9 

340 

CALL  RUN1 

LNK.V 

341 

21  v 

RE-IND  IP3UT 

LNK.9 

34* 

CALL  MAP 

LNK9 

343 

IP0LT -18T (MIN) 

LNK9 

3*»4 

IFI ICaN)2111 .2112 Ill 

LNK9 

3**3 

2112 

NTAPES-NTAPLS-1 

LNKV 

3*»0 

IC0N-1 

LNK.9 

3*7 

oa  T0  2113 

lNkv 

3**o 

2111 

IdT (MINI -AIN2 

LNk.9 

3*9 

2113 

REWIND  iPeuT 

LNK9 

330 

NZ-NZ-MIN 

LNKV 

361 

X1-X3 

LNX9 

33X 

X2-X4 

LNX9 

333 

X3-X2*DELTAX 

LNX.9 

334 

XA-X3*eZ2 

LNKV 

33S 

CALL  SLIDE 

LNK9 

336 

60  T0  300 

LNX9 

337 

ENO 

USA  9 

3>6 

80 


6IBPTC  MAPX  LIST .DECK .M94/2 

MAPX 

0 

SUBROUTINE  MAP 

MAPX 

1 

C 

26  FEB  67 

MAPX 

2 

C 

T.W.SCKWENXC 

TECHNICAL  OPERATIONS  RESEARCH  SR  MAP 

MAPX 

3 

C 

MAPX 

4 

C  « 

•••MAPX 

3 

c 

MAPX 

6 

COMMON  /SlU/ 

MAPX 

7 

1  DIAM 

•DCT 10(123 . IRISE 

•  1EXEC  •  I  SI N 

•  ISOUT 

.  MAPX 

0 

?  so 

.  SPAR  •  5 SAM 

•  TME  .  TMP1 

.  TMP2 

•  MAPX 

9 

3  7XM 

.  y  •  VPR 

•  W  •  HBURST  •  SCLDHB 

*  MAPX 

10 

A  710(40) 

»  RMIN  .  IDISTR 

»  SPAR1  .  H6TAPE  .  FSUM 

»  MAPX 

11 

9  SPAR4 

»  SPARS  •  SPAR6 

*  XGZ  •  YGZ 

.  TGZ 

MAPX 

12 

COMMON  /SETS/ 

MAPX 

13 

1  82 

•  BZ2 

•  BZZ 

•  BZ22 

MAPX 

14 

2  tOELTAX 

•pox 

•  DGY 

•DIFCON 

MAPX 

15 

9  .01 FADJ 

•fmasisooi 

. FMASS(200) 

.  IC(  18 ) 

MAPX 

16 

4  *!C«N 

•  ICTR 

»IH 

. 10T( 18 ) 

MAPX 

17 

5  *IP 

.IPOUT 

•  ITT! 18  ) 

»IV 

MAPX 

18 

6  » JC ( 18  ) 

•  JIN 

»JOUT 

* JPOUT 

MAPX 

19 

7  •KTRI500) 

•ktape 

.LAST 

.MAPRUN 

MAPX 

20 

8  .MARRAY 

.MIN 

•MXREO 

MAPX 

21 

9  »N 

•  NA 

•  NBZX 

.NBZX2 

MAPX 

22 

1  *NBZY 

•NCL 

•  NE 

.NF 

MAPX 

23 

2  .NIJ 

•nmap 

.  NMAX 

.NOX 

MAPX 

24 

3  »NP! 21 ) 

•nreq 

•  NS 

.ntapes 

MAPX 

23 

4  .ntapet 

•ntasx 

•NXMAP 

.NYMAP 

MAPX 

26 

3  *YMJN 

•PS ( 500 ) 

•PSIZEI200) 

•PACT! 200) 

MAPX 

2? 

6  .ROPART 

•SVI2C0) 

•T 1300) 

.T1 

MAPX 

26 

7  *T2 

.TLIMIT 

•  X 1300 » 

.XF 

MAPX 

29 

e  *xe 

•  XMAX 

•  XM1N 

•XNMAP 

MAPX 

30 

9  .XI 

•  X2 

»X3 

.  X4 

MAPX 

31 

1  «Y (SCO  > 

•  YF 

.YO 

*YMAX 

MAPX 

32 

COMMON  /SF74/ 

Of 'API  40T*C  1 

MAPX 

33 

DIMENSION  JMAP!20) 

MAPX 

34 

L 

MAPX 

13 

r  ** 

***M*PX 

16 

c 

MAPX 

17 

i 

FORMAT (IH1.5HSTRIPI3 i 

MAPX 

38 

2 

FORMAT ( / 1 X  *  1 9 1 6 ) 

MAPX 

39 

3 

FORMAT  1 15X21HTW0-L INC  C  FORMAT  MAPi 

MAPX 

40 

4 

FORMAT (3X.19F6. 

3) 

MAPX 

41 

3 

FORMAT!  15X26MTW0-UNE  FH.3  FORMAT  MAP.) 

MAPX 

42 

6 

FORMAT! 16MOOISPLAY  METHOD  14.33m 

IS  NOT  AVAILABLE.  USED  METHOD 

1.IMAPX 

43 

7 

FORMAT  1 // 1 3X  »28mTHE  OUTPUT  PRESENTATION  IS  A) 

MAPX 

44 

8 

FORMAT l//15X»25HTHf  QUANTITY  PRESENTED  IS) 

MAPX 

45 

9 

FORMAT! 15X.27MA 

COUNT  0*  grounded 

WAFERS. ) 

MAPX 

46 

10 

FORMAT  1 13X.38HDOSE  RATE  NORMALIZE 

D  TO'  TIME  M*  1  HOUR.  ) 

MAPX 

.7 

11 

FORMAT  USX .2 DMDOSE  RATE  AT  TIME  H.F10.1.9M  SECONOS. ) 

MAPX 

46 

12 

FORMAT  115* .32HD8SE  ACCUMULATED  3E  TWEEN  TIME  haF10.1.22H  SECONDS 

ANMAPX 

49 

ID  INFINITY.) 

MAPX 

50 

13 

FORMAT! 13X.32HOOSE  ACCUMULATED  BE  TWEEN  TJME  M*Fl0.i»12H  AND  TIME  HMAPX 

31 

1 *F10. 1 »9H  SECONOS.) 

MAPX 

52 

14 

FORMAT! 15*. S4MT0TAL  MASS  OF  OEPOS 

ITEO  PARTICLES. 

) 

MAPX 

53 

13 

FORMAT U5X.44MT0TAL  PARTICLE  MASS 

deposited  between  times  F10.1 

. 5HMAPX 

54 

1  AND  FJ0.1.9M  SECONDS.) 

MAPX 

55 

16 

FORMAT  I /lX»r7.0 

F12.0.3X) .20X.5 

MAPX 

56 

17 

FORMAT! 13X. 41 MASSUM*S  ALL  PARTICLES  APE  GROUNDED 

BY  Tl.i 

MAPX 

57 

18 

rSNMAT ! 1 5X . 1 7MACT IVf TV  aT  TjNF  FI0.1.19M  OUE  TO 

mass  chain  t4) 

MAPX 

50 

19 

FORMAT !15X.26MMULTIPlE  BURST  BINARY  TAPE) 

MAT  X 

59 

81 


rw*»  r%  r% 


20 

FORM AT ( 1 3  X » ) 1 MGR0UN0  ZERO  IS  LOCATED  AT  X 

■  F10 

•  1 »8H  .  Y  -  F10. 

1MAPX 

60 

11 

MAPX 

61 

21 

FORMAT ( 1H1 .41X.36HY-C00RDINATE  SCALES  FOR  SIDES 

OF  MAP/ 1 HOI 

MAPX 

62 

22 

FORMAT ( //IX .F13.0.I2X. FI 3.0) 

MARX 

63 

23 

FORMAT (19X.48HT1MC  (SECONOS)  OF  ONSET  if  FALLOUT 

DEPOSITION.) 

MAPX 

64 

4 

FORMAT fl9X»90HT IKE  (SECONDS)  OF  CESSATION  OF  FALLOUT  DEPOSITION.) 

MAPX 

65 

29 

FORMAT  1 1 9X • 9  OHO  I  /.METER  (MICRONS)  OF  SMALLEST  DEPOSITED  PARTICLE.) 

MAPX 

66 

n 

FORMAT! 13X.49H0 ! AMETER  (MICRONS)  OF  LAROEST 

DEPOSITED  PARTICLE.) 

MAPX 

67 

21 

FORMAT!  19X.90HMASS  DEPOSITED  (KOM/MM*)  BY 

PARTICLES  IN  THE  SIZE  RMaRX 

60 

1ANGE  .F12.9.4H  TO  .F12.5.  9H  MICRONS.) 

MAPX 

69 

28 

FORMAT ( 19X .73HM41  HOUR  NORMALIZED  DOSE  RATE 

RESULTING  FROM  PARTICtMAPX 

70 

ITS  IN  THE  SIZF  RANG'  •F12.5»4*  TO  .F12.5.9H 

MICRONS.) 

MAPX 

71 

C 

MAPX 

72 

c  «• 

*mapx 

c  ••• 

•mapx 

74 

c 

MAPX 

75 

DATA  B’TLUM* INC.LREW/  6HMULTJB.i9.!V 

MAPX 

76 

c 

MAPX 

77 

IF(MAPRUN)  101.100.101 

MAPX 

78 

100 

TINC«9.?»D0X 

MAPX 

79 

XC00RD-XMIN496X 

MAPX 

30 

VINC-INC 

MAPX 

81 

XC INC"VINC*DGX 

MAPX 

82 

KKL»NBZX2*1 

MAPX 

83 

NX«NXKAP-N8ZX2-NBZX? 

MAPX 

84 

c 

LEFT  IS  USED  HERE  AS  A  Tf  'ORARY  STORAGE 

MAPX 

85 

LEFT  » ( XMAX-X2 ) /9GX 

MAPX 

86 

c 

PRINT  MAP  TITLE 

MAPX 

87 

WRITE  (ISPUT.T) 

MAPX 

88 

c 

SELECT  APPROPRIATE  DISPLAY  OP,  N  CODE 

MAPX 

89 

IF(JC (1)1147. 147. 131 

MAPX 

90 

131 

I F( JC ( 1 ) ~6 1132.132. 147 

MAPX 

91 

13C 

JCI 1 )»i 

MAPX 

92 

13? 

N1 » JC  >  1 ! 

MAPX 

93 

Cl  It  (141. 142. 143. 144. 141. 1461 »N1 

MAPX 

94 

141 

ASSIGN  150  TO  N2 

MAPX 

95 

WRITE  ! I S9UT  »3  J 

MAPX 

96 

Cl  TO  1C2 

MAPX 

97 

142 

ASSIGN  151  TO  N2 

MAPX 

98 

WRITE  (ISOL‘T.3  5 

MAPX 

99 

Gf  T?  102 

MAPX 

100 

143 

WRITE! ISOUT. 19) 

MAPX 

101 

ASSIGN  301  TO  N2 

MAPX 

102 

IFILREW.NE.01  GO  T?  1431 

MAPX 

IC3 

LRfW«’. 

MAPX 

104 

REWIND  W8TAPP 

MAPX 

105 

14«l 

WRITE  (W9T.APF  >SITLCv 

MAPX 

106 

WRITE (MSTaPE ) XKIN.XMAX.YMIN.YMAX »OG*»rGY 

MAPX 

1C7 

GO  TO  102 

MAPX 

108 

c 

MAPX 

109 

£•••••••••••••••••••••«••  |h$r-<Tl^3  ►'a  1  NT  5  •' 

*«apx 

110 

1  *4 

CeNTINUE 

MAPX 

111 

14* 

CONTINUE 

MAPX 

112 

146 

CONTINUE 

MAPX 

113 

j-e^r  |N$f ST  1# N  POINTS  *< 

•MAPX 

114 

C 

MAPX 

115 

147 

WRITE  CISeuT.6iNl 

MAPX 

116 

GO  TO  l?p 

MAPX 

117 

101 

KKl«l 

MAPX 

118 

NX«NXMAP-N6ZX2 

MAPX 

119 

c 

CERT  1$  USED  MERE  AS  A  TEMPORARY  STORAGE 

MAPX  120 

ICET-IXMAX-XIJ/DGX 

MAPX  Hi 

TO  1702 

MAPX  122 

C  102  PRINT  ORDINATE  DESCRIPTION 

MAPX  123 

C 

MAPX  124 

102 

WRITE  (IS0UT.8) 

MAPX  123 

00  TO  (161.162.163tl64.165tl66.167.169.169.171.172.173 

.174.175. 176MAP*  126 

1.177.179.179.170.1701 .NREO 

MAPX  127 

161 

WRITE  (1S0UT.9) 

MAPX  126 

GO  TO  170 

MAPX  129 

162 

WAITE  USOUT.lOl 

MAPX  no 

GO  TO  170 

MAPX  131 

163 

WRITE  ( ISOUT .11 ) T 1 

MAPX  132 

GO  TO  170 

MAPX  133 

164 

WRITE  (ISOUT. 121T1 

MAPX  134 

GO  TO  170 

MAPX  13S 

169 

WRITE  (ISOUT. 131T1.T2 

MAPX  136 

GO  TO  170 

MAPX  137 

166 

WRITf  (ISOUT. 14) 

MAPX  138 

GO  TO  170 

MAPX  139 

167 

WRITE  (I50UT.191T1.T2 

MAPX  140 

GO  TO  170 

MAPX  141 

168 

WRITE  (ISOUT. 131T1.T2 

MAPX  142 

WRlTEdSOUT.17) 

MAPX  143 

GO  TO  170 

MAPX  144 

169 

WRITE  (ISOUT. 121T1 

MAPX  145 

WRITE  (ISOUT.17) 

MAPX  146 

GO  TO  1 70 

MAPX  147 

17J 

WRITE  ( IS0UT.181T1 .MASCHN 

MAPX  148 

WRITE ( ISOUT .17) 

MAPX  149 

GO  TO  170 

MAPX  1 5C 

17? 

WRITr  nseL‘T.23) 

MAPX  151 

GO  TO  170 

MAPX  152 

173 

WRITE  ( I  Sfi.'T  >24  ) 

MAPX  153 

GO  TO  1 7C 

MAPX  154 

174 

WRITE  (ISOUT, 25) 

MAPX  155 

Ge  TO  170 

MAPX  156 

175 

WRITE  I  I SCUT  .26 ) 

MAPX  157 

GO  Te  170 

MAPX  158 

176 

WRITE  (ISOUT ,/l)  T1.T2 

MAPX  1d9 

GO  TJ  170 

MAPX  160 

177 

WR I Tr  (ISOUT. ?8)  T1.T2 

MAPX  161 

GO  T?  170 

MAPX  162 

r 

V 

MAPX  163 

£••••« 

178 

CONTINUE 

MAPX  165 

179 

CONTINUE 

MAPX  166 

(ttMl 

c 

MAPX  163 

17C 

WRIT"  I ISOUT »?:i  »G2.TG^ 

MAPX  169 

I E (  JC  ( 1  ) -EC?.  3  i  03  TO  1  :C2 

MAPX  1 70 

c 

MAPX  171 

c 

PRINT  A  PAIR  OF  PASTE -3.N  Y  SCAi.ES  HERE 

MAPX  172 

WRITE  (ISOUT, 21) 

MAPX  173 

(fjYVAPl 

MAPX  174 

Of  1701  J«1 .NYVAP 

MAPX  175 

WRITr  I ISOUT ,22)v» •»» 

MAPX  176 

!7ri 

YY»YY-rOY 

MAPX  177 

170? 

irnffT-Nii  ic2i ,io:2.:c22 

MAPX  178 

1071 

NX»Ef FT 

MAPX  179 

S3 


v  n  r>  n  r%  r> 


102?  MH-NXMINC1 

C  LEFT  IS  USED  HERE  AS  THE  NUMBER  fF 
C  PRINTER  STRIP 

LEFT»NX-HM*I INC! 

STRIPS 

Of  110  ISTRIP-.I.M 

MAPRUN«MAPRUN+1 

IF  (JCm.EQ.3)  G#  Tf  102) 

XC2-XC«PR04TINC 

XC3*XC2*T| NC 

WRITE  I ISBUT  #1 1MAPRUN 

write  useuT.i6ixce»R0*xc2*xc) 

102)  XL-XXL*INVVAl»-n»NXMAP 


MAP*  180 
MAPX  181 

PRINT  COLUMNS  IN  THE  LAST  MAPX  182 

MAPX  183 
MAPX  18* 
MAPX  183 
MAPX  186 
MAPX  187 
MAPX  188 
MAPX  189 
MAPX  190 
MAPX  191 
MAPX  192 
MAPX  19? 
MAPX  194 


IF! I STRJP-MI103.1 04.103 

MAPX 

195 

10* 

AINC-LCFT-1 

MAPX 

196 

VLEFT-LEFT 

HA*>X 

197 

XCIN*VlEFT*DGX 

MAPX 

198 

G«  Tf  1C31 

MAPX 

159 

10) 

KINC«INC-1 

MAPX 

200 

XC I N»XC INC 

MAPX 

2C1 

10)1 

CfNTINUE 

MAPX 

202 

XL  INK  •  K.INC  +  1 

MAPX 

203 

I F  (  JC 1 1  )  .00.3  )  WRITE  (M8TAPC  1NYI1AP  .XL  INX 

MAPX 

29* 

MAPX 

2'5 

RfWS 

MAPX 

206 

OZ  2CC  J-l.NYMAP 

MAPX 

207 

KH«KL*K INC 
K9C-C 

92  201  X-KL.XH 
2A1  FSUM»FSL'y*2MAP«M 

numbers  within  Rews 

30?  K»KL  tXH 
X2C"X0C*1 

TRANST  ER*TC  ceoE  FfR  SELECTED  PRESENTAT I JN 
GZ  Te  N2. 1159.151  .301  ) 
c 

C  150  CZ9E  E?R  PewfB  ZF  Tr«(  DISPLAY 
1 5C  IFlfVAPIXI 1125. 1C6. 127 
105  ASSIGN  1?1  la  N) 

PVAP ( < I *-ZVAP(K  ) 

GC  T?  1 r  9 


MAPX  203 
MAPX  209 
MAPX  219 
MAPX  211 
MAPX  212 
map*  213 
MAPX  214 
MAPX  213 
MAPX  216 
MAPX  217 
MAPX  213 
MAPX  219 
MAPX  220 
MAPX  221 
MAPX  272 
MAPX  223 


107  ASSIGN  ??0  TZ  S3 
109  H  *  ALZG1 0<  CV*.P  (XI! 


MAPX  22 4 
MAPX  225 


Hi »AM?9(M. 1 »C ) 

JMAPIK9C 1 «H-M1 

ir ijMAPtxoc i .ro.ot jv \pixoc ! »o 
fyAP(A)«:o.o»*Mi 
IF tZMAP!* 115.12 5. 1091 
1091  CMAPtXJ.OVAPtK  '1C.0 
JMAP I  ADC  1  JM*P  ( «  "C  )  ♦  l 
GZ  T3  115 


MAPX  226 
MAPX  227 
MAPX  228 
MAPX  229 
MAPX  230 
MAPX  231 
MAPX  ?>2 
M*.  'X  2  33 


1*>6  J^APIKDC  >  »o 


MAPX  23* 


GZ  Tf  )0P 

115  C8  T 8  N).t J0C.121  1 
C  121  RESET  SIGN  8F  maP  C.8R9INATE 
121  "HA:  •«. )»  — 2MA?t£  I 
oe  «e  300 


MAPX  235 
MAPI  236 
MA°X  237 
MAPX  23d 
MAP*  239 


84 


n 


f  i 
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MAPX 

2*0 

191 

CfDE  FIR  TWe-U IKIE  FI  1.3  DISPLAY 

HAPX 

2*1 

151 

JMAPfKDC) ■BMAPIK)/10.0 

MAPX 

2*2 

ZKAP-JMAPIKDO 

MAPX 

2*3 

0MAP <  X 1 -0MAP t  X  J  -  < ZMAP»1 0.0) 

MAPX 

2** 

300 

CONTINUE 

MAPX 

2*5 

WRITE  US8UT.2)  <JMAPIX).X«1  .»  , 

MAPX 

2*6 

WRITE  I  IS0UT.A) (8MAPIX) .X*Xl .\M> 

MAPX 

2*7 

GO  T0  200 

MAPX 

2*8 

301 

WRITE  (MBTAPE)  (0MAPIX)  »X»XL«XK) 

MAPX 

2*9 

200 

XL*XL-NXMAP 

MAPX 

250 

it  ucm.eo.’.:  G0  to  no 

MAPX 

251 

WRITE  1 IS8UT.16JXC00RD.XC2.XC3 

HAPX 

252 

XC00RD-XC00RD+XCIN 

MAPX 

253 

110 

XXL»XXL*INC 

MAPX 

25* 

RETURN 

MAPX 

255 

ErlD 

MAPX 
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0 

SUBROUTINE  PROC 

PRC 

1 

'C  26  FEB  67 

PRC 

2 

C  P.  FLUSSER  TECHNICAL  OPERATIONS 

RESEARCH  SR 

PROC 

PRC 

3 

C***THI5  SUBROUTINE  COMPUTES  A  NUMBER t 

KTR<!>»  WHICH 

determines 

PRC 

4 

C***INT0  WHICH  Z0N£  OR  BUFFER  ZONE  THE  I  Th  PARTICLE 

HAS  LANDED, 

PRC 

5 

C»**  IT  SETS  <TRin»0  IF  THE  1  TH  PARTICLE  LANDED  IN 

the  first  zone. 

PRC 

6 

C***  AND  CALLS  SUBROUTINE  CALC  TO  COMPUTE  THE  CONTRIBUTION  TO  THE 

PRC 

7 

C***  FINAL  RESULT  OF 

this  partio-e.  it  also  computes 

A  NEW  VALUE  FOR 

PRC 

8 

C*»*  FOR  NE*  THE  NUMBEC  OF  EMPTY  SPACES  IN  THE  PARTICLE  ARRAY. 

PRC 

9 

C»**  FINALLY.  THE  NUMBER  OF  (ARTICLES 

that  have  fallen  in  each  zone 

PRC 

ID 

C*«»  IS  COMPUTED  <NP<m. 

PRC 

11 

C 

PRC 

12 

Q  %  *##**#•»######**#***#«  ft  ft*************************************  **#***#**Prt£ 

13 

C 

HRC 

14 

COMMON  /SET1/ 

PRC 

16 

1  DIAM 

•OETI D  C 1 2  J  *  I R I SE 

»  lEXEc  •  ISIN  »  tSOUT  , 

PRC 

16 

2  SD 

.  SPaR  .  SSAM 

.  TME  »  TMP1  .  TMP2  , 

PRC 

17 

3  T2M 

.  U  »  VPR 

.  W  .  HBURST  ,  SCLDHB  » 

PRC 

18 

4  TID140) 

.  RMIN  .  IDISTR 

.  SPAR1  .  M8TAPE  *  FSUM  , 

PRC 

19 

5  SPAR4 

,  SPAR5  »  SPAR6 

.  SPAR7  .  SPARS  ,  SPAR9 

PRC 

20 

COMMON  /SET3/ 

PRC 

21 

1  BZ 

•  BZ2 

,BZZ 

»BZ22 

PRC 

22 

2  .DELTAX 

»DGX 

,DGY 

.DiFCON 

PRC 

23 

3  .DIFADJ 

.FMASI500) 

.FMASS ( 200 ) 

►  i ct  ia » 

PRC 

24 

4  .ICON 

» ICTR 

»!H 

» 1 1ST  <  18  ) 

PRC 

26 

5  .IP 

•  I  POUT 

»  ITT ( 18 ) 

,IV 

PRC 

26 

6  .  JC (18) 

•  JIN 

.JOUT 

* JP3UT 

PRC 

27 

7  »KTR i 500  J 

•KTAPE 

.LAST 

.maprun 

PRC 

28 

8  .MARRAY 

•  MIN 

.MXREQ 

PRC 

29 

9  »N 

,NA 

.NBZX 

.NBZX2 

PRC 

30 

1  »NBZY 

»NCL 

,NE 

.NF 

PRC 

31 

2  .NIJ 

,NMAP 

»NMAX 

.NOX 

PRC 

32 

3  «NP( 21 ) 

»NREC 

•  NS 

.NTAPES 

PRC 

33 

4  .NTAPET 

•NTA5K 

.NXMAP 

.NYMAP 

PRC 

34 

5  .YMIN 

»PS ( 500 ) 

,PSIZE(200) 

•PACT ( 200 ) 

PRC 

36 

6  .ROPART 

•SV ( 200 ) 

«T (500) 

»T1 

PRC 

36 

7  .12 

.TLIMIT 

»X  I 500 ) 

.XF 

PRC 

37 

e  .xo 

•  XMAX 

»  XMIN 

.XNMAP 

PRC 

36 

9  .XI 

»X2 

»X3 

*X4 

PRC 

39 

1  »Y ( 500 ) 

.YF 

,Y0 

» YMAX 

PRC 

40 

COMMON  /SET4/ 

0MAP \ 4000 ) 

PRC 

41 

C 

PRC 

42 
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*PKC 

43 

c 

PRC 

44 

DATA  PR0GRM/6H 

PROC  / 

PRC 

45 

c 

PRC 

46 

•-  #V##  *  *****  *** 

##*** ##**#**#«###* 

.**#****»**»-»#**»* 

♦PRC 

47 

PRC 

4b 

C  NGONE  IS  THE  NUMBER  OF  PARTICLES 

DISCARDED 

PRC 

49 

IF(ICTR)  1.1,2 

PRC 

50 

C  INTO  OR  NEAR  ENOUGH  TO  THE  AREA  OF  INTEREST 

PRC 

51 

1  JCTR  «  0 

PRC 

52 

DO  107  I *1 »NI J 

PRC 

53 

IF(XU)-XO)  107 

•108.108 

PRC 

54 

108  IFtXF-XI m  107 

.109.109 

PRC 

55 

109  IF ( Y ( I > -Y0 I  J07 

.110,110 

PRC 

55 

11C  IFIYF-Yini  1C7 

,111,111 

PRC 

57 

111  JCTR-JCTR+l 

PRC 

58 

X! JCTRJ-XU  ) 

PRC 

59 

86 


rtJcT*)*mi 

PRC 

60 

TcjcTR)«rm 

PRC 

61 

T»SUCTW»»PSm 

*ftC 

62 

FMAS( JCTRI*PMAS4 1 1 

PRC 

63 

10? 

CONTINUE 

PRC 

64 

c 

WG0NF  IS  THE  NUMBER  0F  PARTICLES  DISCARDED 

PRC 

65 

WG0Nfc*N! J-JCTR 

PRC 

66 

1FING0NE1  20.2.52 

PRC 

67 

20 

1RR0R*-2C 

PRC 

68 

GO  TO  7734 

PRC 

69 

52 

JPTR=JCTR*1 

PRC 

70 

9IG«X0-1.0 

PRC 

71 

DO  104  I-JPTR.NIJ 

PRC 

72 

KTRU)-0 

PRC 

73 

104 

X(  I  )*BIG 

PRC 

74 

NE=NE+NG@N£ 

PRC 

75 

NIJ*JCTR 

PRC 

76 

IF(.'CTR!2l»15»2 

PRC 

77 

21 

IRR0R»-21 

PRC 

78 

7734 

CALL  errorsprogrm.irror.isouti 

PRC 

79 

c  *****  end  of  particle  (cloud  subdivision!  discarding  code  ********* 

PRC 

2 

DO  9  1*1. NIJ 

PRC 

81 

r*x(  n-xi 

PRC 

82 

K*JIN*2 

PRC 

83 

DO  8  J-2.K.2 

PRC 

34 

R*R-BZ2 

PRC 

85 

!F(R}  3.3.6 

PRC 

36 

3 

IF( J-21  4.4.5 

PRC 

67 

4 

KTR ( I j  *0 

PRC 

88 

1C 

N£*NE+1 

PRC 

89 

11 

IP*I 

PRC 

93 

CALL  CALC 

PRC 

91 

GO  TO  9 

PR<t 

92 

5 

KTR( I )*J-3 

PRC 

93 

NP(J-3)*NP(J-3)+l 

PRC 

94 

IF( J-4) 22 .25  *9 

PRC. 

95 

22 

IRR0R--22 

PRC 

96 

GO  TO  7734 

PRC 

97 

6 

R"R-DELTAX 

PRC 

98 

I F  ( R )  7.7.8 

PRC 

99 

7 

KTR(I»*J-2 

PRC 

100 

IFU-2)  10.10.12 

PRC 

101 

12 

NP( J-2 ) *NP ( J-2 ) +1 

PRC 

102 

GO  TO  9 

PRC 

103 

25 

KTR( I)*2 

PRC 

104 

GO  TO  11 

PRC 

105 

e 

CONTINUE 

PRC 

106 

KTR( I !«JIN 

PRC 

1C7 

NP( J1N)«NP( JINJ+1 

PRC 

108 

9 

CONTINUE 

PRC 

109 

15 

RETURN 

PRC 

110 

END 

PRC 

111 

SIBFTC  RNNl  L !ST»DECK.M94/2  RNNl  © 

SUBROUTINE  RUN1  RNNl  1 

C  26  FEB  67  RNNl  2 

C  P.  FLUSSER  TECHNICAL  OPERATION'  PcSEARCH  SR  RUN1  RNN1  4 

C  RNN1  4 

C»*#THIS  SU8R0UTINE  IS  CALLED  IF  AND  0NLY  IF  ALL  THE  PARTICLES  THAT  RNNl  5 

C##«ARE  0N  THE  TAPE  WHICH  WILL  BE  READ  NEXT  EITHER  FALL  INT0  RNNl  6 

C*#*  THE  AREA  CURRENTLY  BEING  CONSIDERED  OR  CAN  BE  DISCARDED  RNNl  7 

C»**  ALTOGETHER.  THIS  SUBROUTINE  CALLS  CALC  WHICH  THEN  COMPUTES  RNNl  8 


C**«  THE 

CURRENT  REQUEST. 

RNNl 

9 

C 

R*N! 

10 

#***##  #***♦**####*#*####“## 

11 

c 

R'  11 

12 

C0FM0N  /SET1/ 

RNNl 

13 

1 

DIAM 

.DETIDU2)  .IRISE 

»  IEXEC  . 

I  SIN  .  IS0UT 

• 

RNNl 

14 

2 

SD 

.  SPAR  .  SSAM 

»  TM£  . 

TMP1  ,  TMP2 

• 

RNNl 

15 

3 

T2M 

•  U  .  VPR 

.  W  . 

HBURST  .  SCLDH3 

• 

RNNl 

16 

4 

TIO I 40 I 

.  RMJN  ,  IDISTR  .  SPAR1  . 

MBTAPE  »  FSUM 

» 

RNNl 

17 

5 

SPAR4 

♦  SPAR5  .  SPAR6 

.  SPAR7  . 

SPAR8  .  SPAR9 

RNNl 

18 

C0MM0N  /SET3/ 

RNNl 

19 

1 

9Z 

.RZ2 

•  BZZ 

.  BZ22 

RNNl 

20 

2 

.DELTAX 

♦  DGX 

.DGY 

.DIFC0N 

RNNl 

21 

3 

.OIPADJ 

.FMaS ( 500 ) 

•FMASSI20C) 

.  ICC  IB ) 

RNNl 

22 

4 

,  ICON 

♦  ICTR 

.IH 

.I0T<18! 

RNNl 

23 

5 

.IP 

.  I  POI'T 

.ITTU85 

.IV 

RNNl 

24 

6 

»JC( 18) 

.JIN 

» J0UT 

. JPZOT 

RNNl 

25 

7 

.KTRI500) 

.KTAFE 

.LAST 

.MAPRUN 

RNNl 

26 

e 

.MARRAY 

.MIN 

.MXREQ 

?NN1 

27 

9 

•  N 

.NA 

»  NBZX 

“.NBZX2 

3NN1 

28 

l 

»NBZY 

.NCL 

♦  NE 

.NF 

RNNl 

29 

2 

.NIJ 

.NMAP 

*NMAX 

.NOX 

RNNl 

30 

3 

*NP( 21 ) 

.nreq 

.NS 

•NTAPES 

RNNl 

31 

4 

.ntapet 

.ntask 

.NXMAP 

.NYMAP 

RNNl 

32 

c 

.YMIN 

.PS ( 500 ) 

.PSIZEC200) 

.PACT (200) 

RNNl 

33 

6 

.RePART 

.SVI200) 

.T (500) 

*T1 

RNNl 

34 

7 

«T2 

.TLIMIT 

♦XI 500) 

.XF 

RNNl 

35 

e 

.X0 

.XMAX 

•  XMI N 

.XNMAP 

RNNl 

36 

9 

.XI 

»X2 

.  X3 

.X4 

RNNl 

37 

1 

.Y ( 500 ) 

.YF 

,Y0 

»YMAX 

RNNl 

38 

COMMON  /SET4/ 

0MAPI4OOC) 

RNNl 

39 

C 

RNNl 

40 

£  *###*##*##*#tf###*###*#*****#****###*##***#*******#****##****#*##*##***KNNl 

41 

C 

RNNl 

42 

DATA  PR0GRM/6H 

RUN1  / 

RNNl 

43 

C 

RNNl 

44 

£  ###*#*#*#*###*######*#*#*#####*■#**###*#**#*##**##**#**##***##*###**##*^1^ 

45 

e  *#####*##*###*#*#*##*##*####*###**##*##**####***##*#***#####*#**##*#**k.\ki 

46 

C 

RNNl 

47 

c***  are 

THERE  ANY  CURRENT  PARTICLES 

LEFT  IN  CORE. 

•  • 

RNNl 

48 

c***  HAVE 

THOSE  PARTICLES  WHICH  FALL 

0UTSIDE  THE  AREA  0F  INTEREST 

RNNl 

49 

c**»bfen 

DISCARDED... 

RNNl 

50 

10  I F  < ICTRJ2.3.4 

RNNl 

51 

2  IRR0R— 2 

RNNl 

52 

7734  CALL  ERR0R(PR0GRM»IRR3R,IS0UT) 

RNNl 

53 

3  ASSIGN  200  TO  N 

RNNl 

54 

G0 

TO  5 

RNNl 

55 

C*»«ARE  THERE  ANY  CURRENT  PARTICLES  LEFT  IN  CORE. • 

• 

RNNl 

56 

4  IF(NIJ)41 ,42.42 

RNNl 

57 

h1  ASSIGN  200  T0  N 

RNNl 

59 

N!  J 

-MARRAY 

RNNl 

59 

88 


60  T0  «<J 

RNN1 

60 

42 

ASSIGN  100  Tf  Cl 

RNN1 

61 

T* 

HEAD  l  rPfOTTWTJ 

RNN1 

62 

IFINIJI  6  >7 10 

RNN1 

63 

C«**ART  WC  DONE  •  •  • 

RNN1 

64 

6 

IRR0R— 6 

RNN1 

66 

j0"T§  773* 

RNN1 

66 

8 

READ  ( JP0UT ) tX t I 1 »Y | I ) *T( J ) *PSJ I ) »FMAS (I)iI"l»NIJ) 

RNN1 

67 

00  300  1*1  *PTIJ 

RNN1 

68 

C***PRELIMINARY  CHECK  0F  PARTICLES 

RNN1 

69 

G0  T0  K»f ICO *200) 

RNN1 

70 

200 

IF(X(I)-X1)300. 201.201 

RNN1 

71 

201 

IF(Xin-X4»202. 202*300 

RNN1 

72 

202 

IFirm-Y0>  300.203*203 

RNN1 

73 

203 

IFCYm-YF)  100*100.300 

RNN1 

74 

100 

IP*I 

RNN1 

76 

CALL  CALC 

RNN1 

76 

300 

CONTINUE 

RNN1 

77 

GO  T0  10 

RNN1 

78 

7 

RETURN 

RNN1 

79 

END 

RNN1 

80 

or*  noonoooooono  oooooooooooooooo  n  onnn  or* 


8IBFTC  SHIP  LIST. DECK .M94/2  SHIP 

SUBRIUTINE  SHIFT  SHIP 

26  FEB  67  SHIF 

P.  FLUSSER  TECHNICAL  OPERATIONS  RESEARCH  SR  SHIFT  SHIF 

THIS  SUBROUTINE  WRITES  0N  THE  APPR0PRIATE  TAPE  THE  PARTICLE  SHIF 

PARAMETERS  0F  TH05F  PARTICLES  WHICH  FALL  IN  THE  M0ST  DENSELY  SHIF 

POPULATED  ZONE  0F  THE  AREA  0F  INTEREST.  IT  ALS0  COMPUTES  A  NEW  SHIF 

VALUE  0F  NE.  THE  NUMBER  0P  EMPTY  SPACES  IN  THE  PARTICLE  ARRAY.  SHIF 

SHIF 

SHIP 

NS  •  NUMBER  0F  PARTICLES  T0  BE  WRITTEN  BUT  SHIF 

NE  ■  NUMBER  0F  EMPTY  SPEACES  CURRENTLY  AVAILABLE  SHIF 

KTR(I)  ■  INDEX  INDICATING  INT0  WHICH  ZONE  THE  I  TH  PARTICLE  HAS  SHIF 

LANDED  SHIF 

J0UT  «  INDEX  0F  Z0NE  T0  BE  WRITTEN  OUT  SHIF 

X.Y.T.ID.PSI.  «  PARTICLE  PARAMETERS  SHIF 

I0TIJ0UT/2)  «  TAPE  NUMBER  0F  TAPE  TC  BE  USED  IN  CURRENT  WRITE  SHIF 

KTAPE  *  Tape  number  OF  Tape  TO  be  USED  IN  CURRENT  WRITE  SHIF 

JIN  «  LARGEST  INDEX  APPEARING  IN  PARTICLE  CLASSIFICATION  SHIF 

NMAX  «  MAXIM  I M  NUMBER  0F  UNSSRTED  PARTICLES  T2  BE  WRITTEN  OUT  SHIF 

IN  0NE  DATA  BLOCK  SHIF 

JT  ■  TOP  COUNTER  SHIF 

JB  «  BOTTOM  COUNTER  SHIF 

SHIF 

*•****«»*»•»•»*•»••*«****•••**»»•••••»**•»••••***•»»•••***••••* **#aHjp 

SHIF 


IF  NTH.NBH  OR  NTPH  ARE  ZERO •  THERE  IS  A  HOLE  IN  THE  TOP.  BOTTOM  ORiHlF 


TEMPORARY  STORAGE  RESPECTIVELY.  IF  THESE  VARIABLES  ARE  1.  SHIF 

THERE  IS  NO  HOLE  THcRE  SHIF 

NHCTR  »  INDEX  HELPING  TRACK  OF  SPACE  AVAILABLE  FOR  INSERTION  OF  SHIF 
CONTENTS  OF  Tl MPORaRY  STORAGE  AT  THE  END  3F  EXECUTION.  SHIF 
NPIII  «  NUMBER  OF  PARTICLES  WITH  CLASSIFICATION  NUMBER  I  SHIF 

IF  LAST-O.  THE  LAST  ZONE  HAS  BEEN  SORTED.  IF  LAST-1  THIS  STILL  SHIF 
NEEDS  TO  BF  DONE.  SHIF 

MARRAY  «  DIMENSION  OF  PARTICLE  ARRAY  SHIF 

SHIF 


0 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 
13 
16 
17 

ie 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 
51 

32 

33 

34 

35 

36 
■»  *» 
38 


COMMON  /SET1/  SHIF  39 


1 

DIAM 

iDET ID  < 12 ) .IRISE 

•  iexec 

»  I SIN  .  ISOUT 

»  SHIF 

40 

2 

SD 

.  SPAR  .  SSAM 

.  TME 

.  TMP1  .  TMP2 

»  SHIF 

41 

3 

T2M 

•  U  .  VPR 

.  w 

»  HBURST  .  SCLDHB 

.  SHIF 

42 

4 

TID { 40 1 

.  RMIN  .  IDISTR 

.  SPAR1 

.  MBTAPE  •  F SUM 

»  oHlF 

43 

5 

SPAR4 

,  SPAR5  .  ' 0ARf 

.  SPAR7 

.  SPARS  .  SPAR 9 

SHIF 

44 

COMMON  /SET3Z 

SHIF 

45 

1 

BZ 

»BZ2 

*BZZ 

.B 122 

SHIF 

46 

2 

•deltax 

»OGX 

.DGY 

.DIFCON 

SHIF 

47 

3 

•RIFADJ 

»rMAS(50C) 

.FMASS1200)  » I C  £  1 8  ) 

SHIF 

48 

4 

•  ICON 

*  ICTR 

.IH 

.  IOT ( 18 ) 

SHIF 

49 

5 

•  IP 

. I  POUT 

*  ITT  1 18  1 

* !V 

SHIF 

50 

6 

.JCI18) 

.JIN 

.JOUT 

* JP0UT 

SHIF 

51 

7 

.KTRI50C) 

.KTAPE 

.LAST 

.maprun 

SHIF 

52 

8 

.MARRAY 

.MIN 

•MXREO 

SHIF 

53 

9 

.N 

.NA 

•  NBZX 

•NBZX2 

SHIF 

34 

1 

*NB2Y 

»NCL 

.NE 

»NF 

SHJF 

53 

2 

•  HU 

•  NMAP 

.NMAX 

»N0X 

SHIF 

56 

3 

.NPI21 1 

•nreq 

.NS 

.NTAPES 

SHIF 

57 

4 

.NTAPET 

.ntask 

.nxmap 

.NYMAP 

SHIF 

58 

5 

•  YMIN 

.PSI500I 

.PSIZEI200)  .PACT i 200 ) 

SHIF 

59 

90 


T 


fl 

0 

0 

0 

[1 

0 

0 

I 

I 

I 

I 

I 

B 

0 


(1 

y 


6 

7 

I 

9 

1 


.ROPART 

»T2 

~»Xf 

•XI 

•YI500) 


•5VC200r 
•TLIMIT 
•KM  AX  ' 
»X2 
•YF 


.mooi 

. XI500) 
•  AMIN 
*X3 
*Y0 


»T1 

•  XF 

.XNWAP 

•  X4 

•  YMAX 


COMMON  /SET4/  0MAPI4OOOI 


DATA  PR0GRM/6H  SHIFT/ 


ARE  WE  SORTING  THE  LAST  2CME 
IF! LAST  I  3.3*60 
60  ARE  WE  DUMPING  LAST  ZONE 
6C  IF(  JOl'T-JIN)  3.3.6 
IRR0R’-6 
G0  T0  7734 
SET  UP  COUNTERS 
3  JY  ■  1 

JB  «  MARRAY 
NTH»1 
N3H»1 
N1 ®H»C 

nhctr^c 

IS  THERE  A  HOLE  IN  7 HE  T0P 


3-N0.  6*Y£S 
5"YES.  3«N0 


)))  IF  YES.  EXAMINE  BOTTOM) 


1C  IFIKTHCwY) )  9.12*11 
9  IRROR-9 

7734  CALL  ERRORkPROGRM.lRRfR.lSOUT I 

12  NTH«G 
G0  T?  15 

C  IS  JT  TH  PARTICLE  T0  BE  DUMPED  ...  IF  YES.  EXAMINE  NEXT  T0 
C  PARTICLE.  IF  NOT.  EXAMINE  BOTTOM  PARTICLE) 

11  IF(  KTRI JT  J-JOUT )  14.13.14 

C**»  THIS  PARTICLE  DOFS  N0T  FALL  IN  BUFFER  Z0NE.  ZERO  0UT 
C*«#lDEHTIFICATI0N  NUMBER. 

13  XTP(JT)«0 
G0  777 

14  IF(XTR( JTI-1-J0UT)  15.16.15 
C*«»  THIS  PARTICLE  FALLS  IN  BUFFER  ZONE. 

C***  WRITTEN  2UT  AND  RLTAINET •  INCREASE 
C»»*  BY  ONE. 

16  KTR( JT) *KTRI JT).l 
JT*JT.1 

tFIJT-NS)  10.10.4C 

IS  THERE  A  HOLE  IN  THE  B0TT0M...  IF  YES.  SEE  IF 
PARTICLE  THAT  WANTS  T0  C0ME  D0WN  FROM  THE  T0P.  IF  N0T .  EXAMINE 
PART ICLE • 

IF(NrH)  17.18*17  _  . 

DOES  PARTICLE  WANT  TO  BE  DUMPED...  JF  YESi  MOVE  IT  TO  THF  TOP.  I 
NOT .  EXAMINE  NEXT  BOTTOM  PARTICLE 
IFIKTRIJB))  20.19.20 
NBH-0 

ifinthi  7g.8ea.70 
MOVE  TOP  TO  BOTTOM 
XUB)»X(JT> 

Yf J9)«YI JT) 

T ( JBI *T I JT ) 


777 


IT  NEEDS  TO  BE  BOTH 
IDENTIFICATION  NUMBER 


THERE  IS 
IF  NOT. 


15 


17 
19 

18 

78 


shit 

40 

SH1F 

61 

SMIf- 

« 

SHIF 

63 

SHIP 

64 

SHIF 

65 

SHIP 

66 

•SHIF 

67 

SHIF 

68 

SNIP 

69 

SHIF 

70 

•SHIF 

71 

•SHIF 

72 

SHIF 

73 

SHIF 

74 

SHIF 

75 

SHIF 

76 

SHIP 

77 

SHIP 

78 

SHIF 

79 

SHIF 

80 

SHIF 

81 

SHIF 

82 

SHIF 

83 

SHIF 

84 

SHIF 

85 

SHIF 

86 

SHIF 

87 

SHIF 

88 

SHIF 

89 

SHIF 

90 

SHIF 

91 

SHIF 

92 

SHIF 

93 

SHIF 

94 

SHIF 

95 

SHIF 

96 

SHIF 

97 

SHIF 

98 

SHIF 

99 

SHIF 

100 

SHIF 

101 

SHIF 

102 

SHIF 

103 

SHIF 

104 

SHIF 

105 

SHIF 

106 

SHIF 

107 

SHIF 

108 

SHIF 

109 

SHIF 

110 

SHIF 

111 

SHIF 

112 

SHIF 

113 

SHIF 

114 

SHIF 

115 

SHIF 

116 

SHIF 

117 

SHIP 

118 

SHIF 

119 

91 


'WJVfiPvrcr} . 

FMASI JBI*PHA$I JT) 

KnrrD6r=«TftrjTT  ’ 

NTN»0 

NBH-1 

C  ARE  WE  FINISHED  ..  IF  NOT *  CONTINUE*  IF  YES*  DUMP  PARTICLES 
C  COMPUTE  A  NEW  VALUE  FfR  NE 
C»M  ARE  we  done... 

118  JB«JB-1 

IFUB-NS)  *0.17*1? 

20  IF(KTRlJB)-JOUT)  2* *23*24 
C***  THIS  PARTICLE  DOES  NfT  FALL  IN  BUFFER  ZBNE. 

C»**  XfR8  BUT  IDENTIFICATION  NUMBER. 

23  KTR( JB)"0 
Gt  TB  26 

2*  IFIKTRJ JB)-1-J8U7)  888123*888 
C»»*  THIS  PARTICLE  FALLS  in  (UFFER  ZONE*..  IT  NEEDS  TB  BE  BOTH 
C***  WRITTEN  BUT  A NO  RETAINED.  INCREASE  ID.  NO.  BY  ONE. 

23  KTR(JBI«KTR(JB)*1 

26  IF(NTH)  27*87*27 

27  IFINTPHI  30.77.30 
30  IRROR*-30 

GB  TB  773* 

C  MOVE  TOP  TO  TEMPORARY  STORAGE 

77  XEMP-XUT) 

YEMP*Y ( JT ) 

TEMP-T(JT) 

PEMP*PS ( JT I 
FEMP«FMA$  (JT) 

XEMP-KTRIJT) 

NTPH*1 

C  MOVE  BOTTOM  TO  TOP 

87  X!JT)»XIJ8I 
YUTI-YIJBJ 
T ( JT  J  »T ( JB ) 

PS(JT)>PS(J6) 

FMASUTI-FMASIJBI 

KTR(JT).KTPCJB) 

C*«**  TO  AVOID  DUPLICATION  OF  PARTICLES*  ZERO  OUT  IDENTIFICATION 
C»»»  NUMBER. 

KTRIJBl-0 

N3H»0 

NTH-1 

C***  KEEP  TRACK  OF  EMPTV  SPACE  IN  BOTT.'M 
NHCTR-JB 
GO  TO  777 

AC  JF(NTPH)  67.5000*67 

c***  move  temporary  storage  to  bottom. 

67  X(NMCTR!-XEMP 
YINHCTRI-YEMP 
tinhctrj-temp 

PS  C  NHCTR ) »PEMP 
FMAS(NHCTR)*FEMP 
KTRJ NHCTR) -KEMP 
3000  WRITE  (KTAPE)NS 

WRITE  UTAPEl(X(I)*Y(n.T<n*PS(I)*FMASm*I*l*NS) 

C»»*  ADJUST  T.*  NUMBER  OF  EMPTIES  ANO  THE  PARTICLE  COUNT 
IFIJ0UT-2)  100.102*101 
100  IRR8R—100 

GO  TO  773* 


SHIP  iztr 

SHIF  121 
SHIT  127 
SHIF  123 
SHIF  12* 
SHIF  123 
SHIF  126 
SHIF  127 
SHIF  128 
SHIF  129 
SHIF  130 
SHIF  131 
SHIF  132 
SHIF  133 
SHIF  13* 
SHIF  135 
SHIF  136 
SHIF  137 
SHIF  138 
SHIF  139 
SHI*  1*0 
SHIF  1*1 
SHIF  1*2 
SHIF  1*3 
SHIF  1** 
SHIF  1*5 
SHIF  1*6 
SHIF  1*7 
SHIF  1*8 
SHIF  1*9 
SHIF  150 
SHIF  151 
SHIF  152 
SHIF  153 
SHIF  15* 
SHIF  155 
SHIF  156 
SHIF  157 
SHIF  158 
SHIF  159 
SHIF  160 
SHIF  161 
SHIF  162 
SHIF  163 
5HIF  16* 
SHIF  165 
SHIF  *66 
SHIF  167 
SHIF  168 
SHIF  169 
SHIF  170 
SHIF  171 
SHIF  172 
SHIF  173 
SHIF  17* 
SHIF  173 
SHIF  176 
SHIF  177 
SHIF  178 
SHIF  179 


9t 


102  Mc>*e*NPtmNPm  snrr  iso 

NP(1)»0  SHIP  111 

GO  T»  110  SHIP  182 

101  mJOUT-JMJ  103*103*105  SHIP  1*3 

10S  tRROR—105  SHIP  184 

GO  TO  7734  SH1F  185 

103  NE»NE*NP(J0UT1  SHIP  186 

110  NP(J0UT42l"NP(J0UT+2)*NPf J0UT+1)  SHIP  187 

NPI J0UT*1 1*0  SHIP  188 

NPIJ0UTl*O  SHIP  189 

120  RETURN  SHIP  190 

END  SHIP  191 
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*!BETC  SLID  LIST. DECK .M94/2 

SLID 

0 

subrbut i nc  slide 

SLID 

1 

C 

P.  ElUSSER  TECHNICAL  0PERATI0NS 

RESEARCH 

SR  SLIDE 

SLID 

2 

C 

26  EEB  67 

SLID 

3 

C 

SLID 

4 

C**»$USROUT!NE  SLICE 

WVES  THE  C0NTENTS 

OF  THE  RIGHT  BUFFER  ZONE 

SLID 

5 

c»«*int#  the  left  buffer  zone  and  zeros 

OUT  THE 

REMAINING 

SLIO 

6 

C***ENTRI£$  IN  THE  MAP  ARRAY. 

SLID 

7 

C»<»*.NiX  ■  number  be 

OUTPUT  GRID  PRINTS  IN  THE 

ZONE  ITSELF 

SLID 

8 

c**t 

C BUNT  INC  IN  TH“ 

X  DfR£( TI0N. 

SLID 

9 

c*»« 

NBZX2-  NUMBER  OF  OUTPUT  GRID  POINTS  IN  THE 

BUFFER  ZONE 

SLID 

10 

c»*« 

COUNTING  IN  THE  X  DIRECTION. 

slid 

11 

c*#* 

NX*  NUMBER  0E 

OUTPUT  GRID  P0INTS 

IN  THE  ZONE  ITSELF  PLUS 

SLID 

12 

C«« 

CNE  BUFFER  ZBNE 

COUNTING  IN  THE  X 

DIRECTION 

• 

SLID 

13 

c*»* 

NXMAP*  NUMBER  0E 

OUTPUT  GRID  POINTS  IN  Th£ 

ENTIRE  MAP  .COUNTING 

SLID 

14 

C»#«TWt>  BUFFER  ZONES) 

IN  THE  X  DIRECTION. 

SLID 

15 

c*** 

nymap.  number  0F 

0UTPUT  GRID  POINTS  IN  ENTIRE  MAP  COUNTING 

sue 

16 

C»»* 

IN  Y  DIRECTION. 

SLID 

IT 

c*.« 

FMAP*  VAP  STBRAGE 

SLID 

18 

r 

SLIO 

19 

c 

•SLID 

20 

c 

SLID 

21 

C0MM0N  /SET1/ 

SLIO 

22 

1  DIAM 

.0ETI0(12) .IRISE 

•  IEXEC 

•  I S IN  •  ISOU  • 

SLID 

23 

?  SO 

.  SPAR  >  SSAM 

.  TME 

.  TMP1  .  TMP2  . 

SLID 

24 

3  T2M 

.  U  .  VPR 

>  W 

.  HBURST  .  SCLOMB  . 

SLID 

25 

4  TI0I4PJ 

•  RMIN  .  10 I STR 

.  SPAR1 

.  MBTAPE  .  F SUM  . 

SLID 

26 

5  SPaRa 

.  SPAR5  .  SPAR6 

.  S?AR 7 

.  SPAR8  .  SPAR9 

SLID 

27 

COMMON  /SET  3/ 

SLID 

28 

1  BZ 

,BZ2 

>BZZ 

.3222 

SLID 

29 

?  .DELTAX 

.nGX 

»OjY 

.OIFC0N 

SLID 

30 

3  «DIPADJ 

.FMASI5CC) 

.FMASS120 

01  .IC118) 

SLID 

31 

4  .ICON 

» 1 C  TR 

.  IH 

.131(18) 

SLID 

32 

5  *IP 

. I  POUT 

•  ittubj 

•  IV 

SLID 

33 

6  »JC(1B) 

•  JIN 

•  JiiUT 

.JPUJT 

Slid 

34 

7  *KTR ( 500 ) 

.XTAPE 

.LAST 

.MAPRUN 

SLID 

35 

6  .MARRAY 

•MIN 

•MXREQ 

SLID 

36 

9  «N 

•  NA 

.NBZX 

.N3ZX2 

SLID 

37 

1  .N3ZY 

»NCL 

•  NE 

.NF 

SLID 

38 

2  «NI  J 

•  NMAP 

•  NKAX 

tNOX 

SLID 

39 

3  .NPI21) 

»nrec 

•  NS 

.ntapes 

SLID 

40 

4  .NTaPET 

.ntask 

.NXMAP 

.NYHAP 

SLID 

41 

5  »  YV  |  N 

»PS 1 500 ) 

.PSIZEI200)  .PACT i 200 ) 

SLID 

42 

6  .R0PART 

•SVC200) 

•  T  T  500  J 

.T1 

SLID 

43 

7  »T2 

.TLIMJT 

.X  5  500 ) 

♦  XF 

SLID 

44 

8  .X0 

»XMA.V 

.  XM I N 

.XNMAP 

SLID 

45 

9  .XI 

*X2 

»X3 

.X4 

slid 

46 

1  »Y(5CC J 

.YF 

.YO 

*  YMAX 

SLID 

47 

COMMON  /SET4/ 

0M'AP(40(>0) 

SLID 

48 

c 

SLID 

49 

c  •* 

•SLID 

50 

c  •• 

•*LJD 

51 

c 

SLID 

52 

K*NBZX2*N0X 

sue 

53 

NX*N 

SLID 

54 

N«0 

SUC 

35 

00  303  J*l. NYMAP 

SLID 

36 

00  301  K*1  .NB2X2 

SUO 

57 

NN*N*K 

SLID 

50 

MM*M*K 

SLIO 

>9 

301 


C*** 


SHIFT  BUFFER  ZONE 

SLID 

60 

0MAP(NN)«0MAPIMM> 

SLID 

61 

L»N*NBZ'X2*1 

SLID 

62 

LMN-L+NX-l 

SLIO 

63 

00  502  I-L.LMN 

SLID 

64 

ZER0  BUT  WHATEVER  IS  LEFT 

SLID 

65 

!  0MAPUI-C.O 

SLID 

66 

M»M*NXMAP 

SLID 

67 

l  N»N*NXMAP 

SLID 

68 

RETURN 

SLID 

69 

END 

SLID 

70 
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•  IBTTC  If*  trST.0fCKrtm/2 

ZER 

0 

submutine  zerb 

ZER 

1 

C  Zi  FEB  «7 

ZER 

2 

C  P.  FLUSSER  '  TECHNICAL  fePERATlBNS 

RESEARCH  SR  ZERB 

ZER 

3 

C***  THIS  SUBRfUTJNE  MAKES  RfcBM  FBR  NIJ  PARTICLES 

TB  BE  WRITTEN  IN 

ZER 

4 

C***  THE  PARTICLE  ARRAY  BY  M0VJNG  THE 

PARTICLE  PARAMETERS  INTB  THE 

>- 

o 

3 

CD 

ZER 

r 

C***«F  THE  ARRAY. 

ZER 

6 

C***  jjm  t«P  CBUNTER 

ZER 

7 

C»«*  JB"  BBTTBM  CBUNTER 

ZER 

8 

C 

ZER 

9 

*zer 

10 

CBMMBN  /SET1/ 

ZER 

li 

1  01 AH  *D£  T 1 0 1 1 2 ) • IRI SE 

.  I  EXEC  . 

I SIN  .  IseuT 

• 

ZER 

12 

2  SO  •  SPAR  .  5 SAM 

.  TME  . 

TMP1  .  TMP2 

t 

ZER 

13 

3  T2M  .  U  .  VPR 

»  W  t 

H3URST  .  5CL0H8 

» 

ZER 

14 

*  TI0I40I  •  RUIN  *  I0ISTR 

.  spari  . 

MBTAPE  .  F SUM 

9 

ZER 

15 

5  SPAR4  .  SPARS  •  SPAR6 

.  SPAR7  . 

SPARS  •  GPAR9 

ZER 

16 

CB^MfN  /SET3/ 

ZER 

17 

1  B l  .?Z2 

•  B2Z 

•  BZ22 

ZER 

18 

2  »OELTAX  .OGX 

•  DGY 

.DIFC0N 

ZER 

19 

3  .OIFAOJ  .FMAS ( 5C0 ) 

. FMASSt 200 ) 

•  ICUSl 

ZER 

20 

4  » ICPN  .ICTR 

•  IH 

. I0T( 18  1 

ZER 

21 

3  .IP  . JPBUT 

» ITT  f i8 ) 

•  IV 

ZER 

22 

6  »  JC (18)  .JIN 

*  JBUT 

.JP3UT 

ZER 

23 

7  *KTR  1 500 )  .XTAPE 

.LAST 

•MAPRUN 

ZER 

24 

8  .MARRAY  .MIN 

•MXREO 

ZER 

25 

9  «N  >NA 

.NBZX 

.NBZX 2 

ZER 

26 

1  .NBZY  ,NCL 

•  NE 

.NF 

ZER 

27 

2  .NIJ  .NHAr 

.NMAX 

•  N3X 

ZER 

26 

3  «  NP ( 2 1 )  .NREO 

.NS 

•NTAPES 

ZER 

29 

4  .NTaPET  .ntask 

.NXHAP 

.NYMAP 

ZER 

30 

«  .YHIN  .PSISOO) 

.PSlZtUOOI 

.PACT  I  200) 

ZER 

31 

6  .R0BART  »SV( 200 ) 

.7(500) 

»T1 

ZER 

32 

7  »T2  .TL1MIT 

. X 1 300 ) 

.XF 

ZER 

33 

8  «X3  .XMAX 

.XM!N 

»  XNMAP 

ZER 

34 

9  .XI  .X2 

»X3 

*X4 

ZER 

35 

1  .Y(5?0>  .YF 

.  Y2 

•  Y  MAX 

ZER 

36 

C2MMPN  /SET4/  avAP(*COOI 

ZER 

37 

DATA  PR8GRM/6H  ZFRe  / 

ZER 

38 

JT.1 

ZER 

39 

JB*HARr,AY 

ZER 

40 

C*»*  IS  ThERF  a  H3lr  IN  The  '.0P...  IF 

YES.  g<:  UN. 

IF  NjJT.  CHECK 

ttom. 

ZER 

41 

2  IFIKTR(J)  13.4.5 

ZER 

42 

3  IRRSR  — 3 

ZER 

43 

7734  CALL  ERR?R(PRBGRM.IRRCR.153L'T) 

ZER 

44 

a  n*JT*i 

ZEP 

45 

C*«*  ARE  WE  03NE... 

ZER 

46 

IFI JT-NIJ)?.?.1C 

ZER 

47 

c*»*  is  there  m?le  in  the  eeTTen...  if 

YES.  M3 VE  A 

PARTICLE  FK0N  THE 

ZER 

48 

C •••  T3P.  IF  N?Y .  TRY  ANO  FI  NO  6NE . 

ZER 

A? 

5  IFIXT7I JB; 16*Y.S 

ZER 

SO 

4  |RR*R»-6 

ZER 

51 

32  Te  7734 

ZER 

32 

<;•••  «avE  particle  parahpters  and  decrement  bbttzm 

cbunter. 

ZER 

33 

?  XIJP1-XUT1 

ZER 

34 

Y< J8)«Y(JT) 

ZER 

33 

TIjBI-TIJT) 

ZER 

36 

PS i JB ) »PS I JT 1 

ZER 

37 

FMAS  I 35) -FMAS! JT > 

ZER 

38 

KYRI JB)»KTR(JT> 

ZER 

59 

H 

0 

0 


JB-JB-1 

ce  re  k 

8 

C*«»  ARE  THERE  EN3UGH  H0LES.... 
mje-N!J>9.9.$ 

9  IRR0R--9 
G0  T0  7734 
1C  RETURN 
ENO 
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ZER  62 
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ZER  65 
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THE  PRINTED  OUTPUT  OF  THE  PROGRAM 


An  example  of  the  output  is  shown  in  the  following  pages.  This  consists  of 
two  parts:  (1)  a  listing  of  the  particles  Impacted  tape  (IPOUT)  which  is  optional 
(see  p,  49),  and  (2)  a  map.  The  map  has  been  hand -contoured  to  shew  the  limit 
of  particle  deposition  (outside  contour  representing  0  R/h)  and  the  10  R/h  line 
(inside  contour).  Laige  maps  are  divided  into  a  number  of  vertical  strips  and  the 
strips  are  numbered  in  increasing  sequence  for  ease  of  identification  and  assembly 
into  a  complete  map.  This  sequencing  is  accomplished  automatically  no  matter 
how  large  the  map  is.  The  maximum,  strip  width  is  fixed  by  the  geometry  of  the 
printer,  but  the  program  automatically  ascertains  that  the  total  number  of  grid 
points  across  all  strips?  is  that  which  the  user  specified.  In  the  vertical  direction 
on  the  paper  (north-south)  all  required  data  points  are  simply  printed  in  their 
correct  positions  and.  therefore,  no  further  identification  is  required.  It  should 
be  .‘oted,  for  the  sake  of  correctly  interpreting  the  output  map,  that  the  coordi¬ 
nates  of  the  lowest,  leftmost  map  data  point  are  not  the  minimum  coordinate  pair 
put  in  by  the  user;  they  are  the  minimum  coordinate  pair  incremented  by  one  x 
and  y  grid  interval.  Thus,  this  lower  left  point  is  the  first  data  point  falling  one 
full  grid  interval  within  the  area  of  interest  to  the  user. 
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OUTPUT  PROCESSOR  MODULE 
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TECHNICAL  OPERATIONS  RESEARCH, INC. 
BURLINGTON,  MASS. 


I 


LISTING  OF  GROUNDED  PARTICLES 


TRANSPORT  IDENTIFICATION 


I 


FIFTH  LARGE  SCALE  TEST  OF  THE  DcLFIC  MODEL.  1  FEB.  1967,  TBANSPOST 


block  i 


I 

f 

I 


B 

P 


NO.  OF  PARTICLES  IN  THIS  BLOCK  IS  150 


X 

Y 

T 

PS 

MASS 

999813.7031 

999866.8594 

109.7602 

661.0439 

0.0090 

S',9972»2656 

999950.7734 

153.9715 

661.0439 

0.0C90 

999970.7891 

999951.6797 

150.3999 

661.0439 

0.0090 

1000027.9297 

999938,4609 

238.7078 

661.0439 

C.0090 

999662.468b 

999698,2500 

217.6443 

661.0439 

0.0090 

999663.429/ 

999698.6719 

218,0461 

661.0439 

0.0090 

999877 .9453 

999762.4375 

270,3820 

661.0439 

0.0090 

999956.3516 

999877.6875 

258.7040 

6ul .0439 

0,0090 

999836.9375 

999713.0938 

304,3261 

661.04J9 

0.0090 

999806.2500 

999705.2266 

290.9941 

661.0439 

0.0090 

100015H.4_.-3l 

999878.6328 

330,9649 

661.0439 

0.0090 

1000141.9453 

999877.0000 

320,445s 

661.0439 

0.0090 

1000103, -'*141 

999840.0781 

355.6223 

661.0439 

0.0090 

1000U  4.9844 

999879.2813 

335.1282 

661.0439 

0.0090 

1000040.3928 

999782.6016 

365.5172 

461.0439 

0.0090 

1000086.0625 

999809.9453 

400.1892 

661.0439 

0,0090 

1000C64.4609 

999804.8984 

387.9202 

661.0439 

0.0090 

1000341.3359 

999952.9609 

417.7250 

661.0439 

0.0090 

1000336.1484 

999951.8359 

414,3079 

661.0439 

0,0090 

1000511.2344 

1000045.2578 

441.2204 

661.0439 

0.0090 

100C5 1 3.3906 

1000045.8125 

442,6173 

661.0439 

0.0090 

1000350.3750 

999944.4766 

467,4481 

661.0439 

0.0090 

1000370.6328 

999960.4609 

460,4454 

661.0439 

0.0090 

1000618.8516 

1000088.9922 

431.9517 

661.0439 

0.0090 

,000“’? . ’““1 

1000206.9531 

504.9391 

661.0439 

0.0090 

1000823.4531 

1000202.7/34 

501.5282 

661.0439 

0.0090 

[1 

B 

S 


101 


1000922.6719 

1000252.3125 

521.0816 

661.0439 

0.0090 

1001043.8829 

1000316.1406 

541.8622 

661,0439 

0.0090 

1001034.0156 

1000311.8759 

538.6753 

661.0439 

0.0090 

1001116.6328 

1000350.0391 

558.5450 

661.0439 

0.0090 

1001293.7109 

1C00514.4766 

559.4243 

661.0439 

0.0022 

1001293.7109 

1000177. 1ST* 

559.4243 

661,0439 

0.0022 

1000956.4141 

1O0O177.U75 

559.4243 

661.0439 

0.0022 

1000956.4141 

1000514*4766 

559.4242. 

661.0439 

0.0022 

999779.9453 

999855.3594 

108.1759 

540.6692 

0.0090 

999662.8516 

999700-3594 

215.4246 

540.6692 

0.0090 

9999*0.2422 

999855.7188 

300.7138 

540.6692 

0.0090 

999942.0234 

999860.6953 

292.7125 

540.6692 

0.0090 

999942.3906 

999861.7266 

291.0587 

540.6692 

0.0090 

999835.2656 

999764.7969 

303.9042 

540.4692 

0.0090 

1000310.9297 

999921.4219 

396.4778 

540.6692 

0.0090 

1000287.4766 

999921.4141 

379.9917 

540.6692 

0.0090 

1000281.7109 

999921.9766 

375.6066 

540.6692 

0.0090 

999970.1563 

999686.3359 

436.9803 

540.6692 

0.0090 

1000015.8828 

999719.2344 

428.3916 

540.6692 

0.0090 

1000466.8594 

999978.7656 

481.6200 

540.6692 

0.0090 

1000286.0391 

999888.2109 

455.6174 

540.6692 

0.0090 

1000518.9219 

1000003.6953 

499.9505 

540.6692 

0.0090 

1000493.5938 

999994.2266 

490.0412 

540.6692 

0.0090 

1000362.9453 

999905.4375 

526.5606 

540.6692 

0.0090 

1000376.5038 

999914.9297 

522.92C3 

540.6692 

0.0090 

1000880.8281 

1000192. *253 

555.3261 

540.6692 

0.0090 

1000883.1719 

1000190.8828 

554.7974 

540.6692 

0.0090 

1000631.6953 

» 030043. 5469 

548.0274 

540.6692 

0.0090 

1001035.7188 

1000251.0547 

580.5498 

540.6692 

0.0090 

1001186.9922 

1000330.3047 

615.6166 

540. 6692 

0.0090 

1001174.6406 

1000325.7813 

611.5813 

540.6692 

0.0090 

1001587.5938 

1000627.9844 

643.7514 

540.6692 

0.0090 

1001523.6016 

1000625.4375 

643.1794 

540.6692 

Q.0090 

1001816.2266 

1000778.6563 

675-0872 

540,6692 

0.0090 

1001998.7031 

1000954.3125 

673.0941 

540.6692 

0.0022 

1001998.7031 

1000617.0156 

676.0941 

540.6692 

0.0022 

1001661.4063 

1000617.0156 

676.0941 

540.6692 

0.0022 

1001661.4063 

1000954.3125 

676.0941 

540.6692 

0.0022 

1002215.5547 

1001094. 7}09 

707.4751 

540.6692 

0.0022 

1002215.5547 

1000757.4141 

707.4751 

540.6692 

0.0022 

1001878.2656 

1000757.4141 

707.47 51 

540,6692 

0.0022 

1001878.2656 

1001094.7109 

707,4751 

540.6692 

0.0022 

1002200.6953 

1001074.4219 

732.1252 

540.6692 

0.0022 

1002200.6953 

1000737.1250 

702.1252 

540.6692 

0.0022 

10016e3.3984 

1000737.1250 

702.1252 

540.6692 

0.0022 

1001863.3984 

1001074.4219 

702.1252 

540.6692 

0.0022 

1002505.4766 

1001082.7578 

733.0195 

540.6692 

0.0022 

1002505.4766 

1000745.4609 

733.0195 

540.6692 

0.0022 

1002034,2031 

1000847.0000 

731.0921 

540.6692 

0.0022 

1002036.2031 

1001184.2969 

731.0921 

540.6692 

0.0022 

1002408.0/81 

1001090.6641 

722.2159 

5*0.6692 

0.0022 

1002408,0781 

1000753.3672 

722.2159 

540.6692 

0.0022 

1001982.2656 

1000816.6563 

720.8784 

540.6692 

0,3022 

.1001982.2656 

1001153.9531 

720.8784 

540.6692 

0.0022 

1002876.6484 

1001161.4609 

740.9310 

540.6692 

0.0022 

1002876.6484 

‘.000824.1641 

749.9310 

540.6692 

0.0022 

1002440.8984 

1000950.7969 

751.4804 

540.6492 

0,0022 

1002440.8984 

1001288.0938 

751 ,4804 

540.6692 

0.002*2 

999890.5156 

999902.7109 

196,0140 

455.9091 

0.0090 

999681. 0859 

999708.1719 

230.6680 

455.9091 

0.0090 

999935.2891 

999831.8516 

286.7701 

455.9091 

0.0090 

999932.1406 

999842.4141 

335.8081 

455.9091 

0.0090 

999719.0313 

999599.5781 

383.2466 

455.9091 

0.0090 

1000218.9531 

999755.1172 

510.1834 

455.9091 

0.0090 

1000217.0000 

999761.1094 

503.5328 

455.9091 

0.0090 

1000216.7422 

999761.9063 

502.6511 

455.9091 

0.0090 

1000218.5078 

999756.4844 

508.6675 

455.9091 

0.0090 

1000436.6328 

999919.3047 

565.0691 

455.9091 

0.0090 

1000513.9531 

999961.4063 

544.4145 

455.9091 

0.0090 

1000672.2578 

1000006.0469 

606.9872 

455.90C1 

0.0090 

1000541.6172 

999931.5469 

597.7697 

455,9091 

0.0090 

999993.9063 

1000001.7734 

3.5909 

12473.7008 

0.0090 

999998.6875 

1000001*9922 

3.94* 

12473.7008 

0.0090 

999998.6953 

1000001.9844 

3.0341 

12473.7008 

0.0090 

999998.6484 

1000002.0469 

4.1363 

12473.7008 

0.0090 

99999P.3594 

1000002.3516 

4.7463 

12473.7008 

0.0090 

999998.1641 

1000002.5703 

5.2631 

12473. 70C8 

0.0090 

999990.0156 

1000024.9922 

53.2560 

661.0439 

0.0090 

1000379.0938 

1000004.4219 

395.6474 

850.1669 

0.0090 

1000510.0625 

1000086.4375 

379.1820 

850.1669 

0.0090 

1000373.4609 

1000010.7031 

363.5162 

850.1669 

0.0090 

1000488.1719 

1000076.0469 

369.2923 

850.1669 

0.0090 

1000351.4375 

1000003.2734 

352.2446 

850.1669 

0.0090 

1000355.7344 

1000004.7266 

354.4451 

850.1669 

0.0090 

1000119.0391 

999873.4531 

335.7248 

850.1669 

0.0090 

1000097.4609 

999865.6875 

318.4868 

850.1669 

0.0090 

1000080.7734 

999856.7344 

323.2555 

850.1669 

0.0090 

1000085.4844 

999863.3672 

304.1440 

850.1669 

0.0090 

1000102.7578 

999890.5781 

286.5992 

850.1669 

0.0090 

1000101.1016 

999877.9297 

295.9191 

850.1669 

0.0090 

959932.3828 

999785.7656 

275.2522 

850.1669 

0.0090 

999875.5547 

999765.7109 

256.5722 

850.1669 

0.0090 

999896.9644 

999789.9844 

239.9431 

850.1669 

0.0090 

999867.7578 

999763.5469 

253.3261 

850.1669 

0.0090 

999916.6250 

999805.0859 

233.3254 

850.1669 

0.0090 

999904. 5000 

999791.4219 

243.5866 

850.1669 

0.0090 

999979.8047 

999875.2422 

220.1690 

850.1669 

0.0090 

1000025.3984 

999947.6094 

199.6659 

850.1669 

0.0090 

999994.2656 

999924.6172 

206.2278 

850.1669 

0.0090 

100CC44.4453 

999967.6719 

184.0057 

850.1669 

0.0090 

1000025.2109 

999948.9531 

196,4965 

850.1669 

0.0090 

999792.9219 

999790.8984 

166.1074 

850.1669 

0.0090 

999772.4453 

999774.8750 

179.2275 

850.1669 

0.0090 

999925.6094 

999920.3984 

140.4025 

850.1669 

0.0090 

999899.8359 

999885.5625 

158.2701 

850.1669 

0.0090 

1000013.5703 

99999'. 6094 

110.8231 

850.1669 

0.0090 

999973.7578 

999956.3516 

126.5207 

850.1669 

0.0090 

1000006.9766 

999980.9219 

117.1562 

850.1669 

0.0090 

999994.1484 

1000020.0234 

50.5690 

850.1669 

0.0090 

999981.1250 

1000025.0703 

26.4409 

850,1669 

0,0090 

1000102.8203 

999930.7500 

229.3430 

1219.1012 

0.0090 

999953.0234 

999835.1875 

219.4342 

1219.1012 

0.0090 

999974.7344 

999852.7344 

209.8566 

1219.1012 

0.0090 

999948.1873 

999843.1016 

201.0213 

1219.1012 

0.0090 

999967.6172 

999850.2969 

207.1759 

1219.1012 

0.0090 

999964.4453 

999854.3984 

196,6471 

1219.1012 

0.0090 

999942.1328 

999848.2031 

186.7100 

1219.1012 

0.0090 

999995.4375 

999904.6328 

177.6155 

1219,1012 

0.0090 
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BLOCK 


2 


NO.  OF  PARTICLES  IN  THIS  BLOCK  IS  1*4 


X 

1000944.4844 

1001193.3281 

1001217.7422 

1001494.0938 

1001388.8594 

1001401.7031 

1002030.1797 

1002055.3828 

1002563.2969 

1002850.9609 

1002850.9609 

1002308.2031 

1002308.2031 

1003109.2031 

1003109.2031 
1002732.5234 
1002732.5234 
1003090.6719 
1003090.6719 
1002707.1719 
1002707.1719 

1003648.9375 

1003648.9375 

1003311.6406 

1003311.6406 

1003535.8281 

1003535.8281 
1003198.5391 
1003198.5391 
1003899.5156 
1003899.5156 

1003562.2266 

1003562.2266 
1003923.1875 
1003923.1875 
1003585.8984 
1003585.8984 
1003926.3672 
1003926.3672 
1003589.0781 
1003589.0781 
1003876.7266 
1003876.7266 
1003539.4297 
1003539.4297 
1003767.3906 
1003767.3904 


V 

999931.8594 

1000307.7031 

1000318.8281 

1000445.9844 

1000401.9141 

1000408.2031 

1000862.6484 

1000875.4219 

1000975.6406 

1001049.3047 
1000712.0078 

1000932.3047 
1001269.6016 
1001158.5859 
1000821.2891 
1000,13.1290 

1001250.4219 

1001138.2344 

1000800.9375 
1000910.0469 
1001247.3438 
1001413.1484 
1001075.8516 
1001075.8516 
1001413.1484 

1001338.2266 
1001000.9297 
1001000.9297 

1001338.2266 

1001505.2969 
1001168.0000 
1001168.0000 

1001505.2969 
1001508.5313 

1001171.2344 

1001171.2344 
1001508.5313 
1001489.6797 
1001)52.3906 
1001152.3906 
1001489.6797 
1001457.6953 
1001120.4063 
1001120.4063 
1001457.6953 
1001385.8672 
1001048.5781 


T 

PS 

MASS 

£.03.1177 

455.9091 

0.0090 

653.7841 

455.9091 

0.0090 

660.8842 

455.9091 

0.0090 

709.7954 

435.9091 

0.0090 

699.7327 

455.9091 

0.0090 

700.7300 

455.9091 

0.0090 

742.8521 

455.9091 

0.0090 

746.3553 

455.9091 

0.0090 

788.0953 

455.9091 

0.0090 

792.0753 

455.9091 

0.0022 

792.0753 

455.9091 

0.0022 

793.1015 

455.9091 

0.0022 

793.1015 

455.9091 

0.0022 

831.9199 

455.9091 

0.0022 

831.9199 

455.9091 

0.0022 

833*4435 

455.9091 

0.0022 

833.4435 

45, .9091 

0.0022 

826.9933 

455.  T'M 

0.0022 

826.9933 

455.9091 

0.0022 

830.0415 

455.9091 

0.0022 

830.0415 

455.9091 

0.0022 

866.0B82 

455.9091 

0.0022 

866.0882 

455.9091 

0.0022 

866.0882 

455.9091 

0.0022 

866.0882 

455.9091 

0.0022 

853.9816 

455.9091 

0.0022 

853.9816 

455.9091 

0.0022 

853.9816 

455.9091 

0.0022 

853.9816 

455.9091 

0.0022 

894.5464 

455.9091 

0.0022 

894.5464 

455.9091 

0.0022 

894.5464 

455.9091 

0.0022 

894,5464 

455.9091 

0.0022 

896.4286 

455.9091 

0.0022 

896.4286 

455.9091 

0.0022 

896.4286 

455.9091 

0.0022 

896.4286 

455.9091 

0.0022 

933.6617 

455.9091 

0.0022 

933.6617 

455.9091 

0.0022 

933.6617 

455.9091 

0.0022 

933.6617 

455.9091 

0.0022 

927.9183 

455.9091 

0.0022 

927.9183 

455.9091 

0.0022 

927.9183 

455.9091 

0.0022 

927.9183 

455.9091 

0.0022 

965.8685 

455.9091 

0.0022 

965.8685 

455.9091 

0.0022 
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«**-»  uiil  PUT  PROCESSOR  lUtNT  I  PICA T ION  **** 

FIFTH  LAAGfc  SCALe  TcST  UP  OcLFIC.  CHECK  TUrtPKlNa.  OPP 


***»  INITIAL  CoN-JII  IONS  (FIREbALLI  iOtNT  IF  ic  AT  IUN  **** 

PIPTr)  LAKGt  SCALE  TcsT  OP  THt  UcLPIl  HUU'.L.  i  Fee.  Iia7.  IN1T.  CUNo. 


****  CLUUO  A I  St  iOENTIFICAl IUN  **»* 

FIFTH  cilit  AwAlE  TeSI  up  Mt  UtLFIC  MUDcL.  i  Ffcb.  19a?.  CLUUU  RISt 


Pan  T ICLc  Sc f  cAPANSiUN  IUENTIFILaI ION  •*** 

FIFTH  LARGi.  SUct  leaf  OF  THc  OtLPiC  MUOEL.  I  Fed.  i9o7.  PSE 


IaaNSPURT  IOcNT IfIlATIoN 

FIFTH  LARGE  SCALE  TeST  OF  THt  OkLFIL  HOUtt.  A  Fto.  19a/.  TRANSPOST 
****  HlNJ  IJtNTIFUAllUN  **** 

THIS  Hi iO  FicLJ  OSes  THfc  FRENCHMAN  FLATS  ANO  koAO  d  SlATlONS.  l/3G/a7 


***»  TUPuGRAPHY  lOkNllFICATION  **** 


»***  OlHtR  INPOIS  •*** 

»**♦  lit  Control  VARIABLE  ARRAY.  IcIJI.  has  GIVcN  THE  FULL Jdi No  VALUES 
1  i  i  n  j  a  7  d  9  10  II  te  ii  In  i»  1  u  0 


Tic  FJLLUhING  LOGICAL  TAPES  Arc  AVAILABLE  ICR  SORTING. 

I  i  a  -u  -o-O  -0  -o  *0  ~o  -0  ~0  -0  -o  -u  “u  -0  -0 
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MIA  X4*X  VMlN  VMAX 
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OfcUA  A 
1SJ0.3 


I-te  uJHTAOL  VAHlAOLt  AAAAV •  iUlJI.  HAS  «6H  UVtN  TH6  fULlJalH*  tfAkJcS. 

1  i  3  4  s  o  7  a  9  10  11  12  U  l*  *1  i  it  0 

•MOUlil  bdUOHNfcSS  r*t.TUH  0.500 


itJj.il  >n»Ha£«  i 

TV«*C  2  Tl  •  1.0  12  •  2.0  HASCMN 

HAHJcJ  CM  OKI 0  Ult-tVALS  OOX  *  1500.0  Oo»»  1<5>.J 

FHc  JJIPOT  fKfcSfci?*  110.0  IS  A 
T4J-ul*e  e  AUXHAf  4»1> 


(He  J-jAHinv  f>tu5t*rtu  !S 

JjSt  4«lt  MOHIAtltcd  U  lint  1  *00*. 

bi4JJ4d  ItM  SS  Ld>-»ffcO  A I  X  »  luJOOOO.Q  »  7  •  toOOJOj.o 


OALTA  V 
1500.0 
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1014250. 

101)000. 

10117)0. 

1010530. 

ICO  ¥2*0. 

100*000. 

1C  06  7  50, 

1C055-70. 

1006  250. 

1005000. 

.act /so. 

1000  500. 

996250. 


1016250. 

1013000. 

1011750. 

lOlusoO. 

1009250. 

1oO#6«3. 

Io0»7  50. 

100»5b0. 

1006250. 

loOJaOC. 

1001756. 

1400500. 

669250. 
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APPENDIX 


PRIMARY  ARRAYS  OF  THE  OUTPUT  PROCESSOR 


OMAP(  ) 

X(J) 

Y(J) 

T(J) 

PS(J) 

FMAS(J) 

KTR(J) 

PSIZEfl) 

PACT  (I) 

F  MASS,!) 

SV(I) 


One -dimensional  map  array  used  to  store  the  part  of  the  map 
currently  being  prepared  for  printing.  Explicit  index  conversion 
is  used  to  store  and  retrieve  two-dimensional  map  data  from  this 
one -dimensional  storage  array. 

X  coordinate  of  the  Jth  central  particle  description  that  is  cur¬ 
rently  in  memory.  The  arrays  X,  Y,  T,  PS,  FMAS,  and  KTR  are 
spoken  of  as  the  particle  arrays. 

Y  coordinate,  t.ee  X(J). 

Impact  time  associated  with  the  Jth  central  particle  description. 

Particle  diameter  (microns)  associated  with  the  Jth  cenrral  par¬ 
ticle  description. 

Mass  per  unit  area  (inks  system)  associated  with  the  Jth  cloud 
subdivision  at  the  time  of  its  definition. 

Class  indicator  number  for  the  Jth  central  particle  description. 
KTR(J)  =  0  indicates  that  the  Jth  particle  description  is  not  in  use. 
KTR(J)  >  0  indicates  the  number  of  the  map  zone  or  buffer  zone 
into  which  the  Jth  central  particle  has  fallen.  See  NP(K,. 

Central  particle  size  (microns,  d'ameter;  of  the  Ith  size  range 
of  the  tabulated  particle  size  V3  property  arrays.  The  uier 
arrays  in  this  jet  are  PACT,  FMASS,  and  SV. 

Minimum  particle  size  of  the  Ith  particle  size  range.  See 
PSIZEfl). 

Fraction  of  the  total  particulate  mass  of  the  cloud  that  is  repre¬ 
sented  by  the  Ith  oarticle  size  range. 

Surface  to  volume  ratio  of  the  Ith  particle  size  range  (as  required 
by  the  Freiling  radial  distribution  model). 


ill 


NP(K) 


CRID(  ) 
OPi:>(  ) 
PSEID(  ) 
TOPID(  ) 
WIDj  ) 
DETID(  ) 
TID(  ) 
IC<  ) 
!OT{  ) 


ITT{  ) 
JMAP(  ) 


Count  of  the  number  of  central  particles  that  are  in  the  particle 
arrays  and  are  known  to  belong  in  the  Kth  numbered  zone  or 
buffer  zone  of  the  desired  mop.  These  zones  and  buffer  zones 
are  numbered  sequentially  from  left  to  right  starting  with  the 
buffer  zone  just  to  the  right  of  the  first  map  zone  to  be  printed 
(Zone  1  of  Figure  1).  Thus  NP(K)  for  odd  K  are  counts  for  buffer 
zones  and  for  even  K  they  are  counts  for  interbuffer  map  zones. 

Cloud  Rise  run  identifier . 

Output  Processor  run  iaentifier. 

CloudRise-Transportlnterface  run  identifier. 

Topography  data  identifier. 

Wind  data  identifier. 

Initial  Conditions  run  identifier. 

Transport  run  identifier. 

Control  integer  array. 

An  array  containing  the  logical  numbers  of  those  tape  units  that 
are  availat-3  for  use  by  the  Output  Processor  in  sorting  grounded 
particle  de.  riptions. 

Permanent  copy  of  the  original  state  of  the  array  IOT(  ). 

Array  internal  to  subroutine  MAP  used  to  store  temporarily  a 
one -line  group  of  integers  for  use  in  the  printing  of  two-line  E  or 
F  format  maps. 
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H  +  T2,  assuming  all  particles  have  arrived  by  H  +  T1  hours;  (8)  same  as 
7  but  integrated  to  'ufinlty;  (9)  concentration  of  an  individual  mass  chain 


Unclassified 

Security  CU»*ific*tioci 


Unclassified 

Security  CUssifiemrioa 


INSTRUCTIONS 


1.  ORIGINATING  ACTIVITY:  Enter  the  name  sad 
»f  the  contractor,  subcontractor,  grantee,  Department  ol 
Dtlnw  activity  or  other  organintioe  (corporate  author) 
issuing  the  report. 

la.  REPORT  SECURITY  CLASSIFICATION:  Eater  the  over- 
ell  nacarity  classification  cl  the  report.  Indicate  whether 
"Raatricted  Data'*  ia  included.  Markieg  ia  to  he  ia  accord¬ 
ance  with  appropriate  eeewity  regulations. 

24.  CROUP:  Automatic  downgrading  ia  a  pacified  ia  OoD 
Directive  5100. 10  and  Armed  Force*  Industrie;  Manual. 

Enter  the  (roup  number.  Alao,  when  applicable,  ahow  that 
optional  markings  have  been  ueed  for  Group  3  end  Group  4 
aa  authorised. 

3.  REPORT  TlTl-E:  Enter  thu  complete  report  title  ia  all 
capital  letter*.  Titles  ia  alt  cases  shoald  be  nncleasified. 

If  a  meaningful  title  cannot  be  selected  without  classifica¬ 
tion,  ahow  title  eUssilicafion  ia  all  capitals  ia  parencheiit 
immediately  following  the  title. 

4.  DESCRIPTIVE  NOTES:  If  appropriate,  enter  the  type  of 
report,  e.g,,  interim,  progress,  summary,  enauel,  or  Ijoal. 

Give  the  inclusive  sates  when  a  specific  reporting  period  ia 
covered. 

5.  AUTHOR(S):  Enter  the  uantefe)  of  author!*)  as  ihows  oa 
or  in  the  report.  Eater  lent  atuae,  first  name,  middle  initial. 

If  military,  ahow  rank  and  branch  of  service.  The  nfiae  of 
the  principal  author  ia  aa  absolute  minimum  requirement, 

6.  REPORT  DATE:  Enter  the  date  of  the  report  a*  day, 
month,  year,  or  month,  year.  If  mor'  than  one  date  appear* 
on  the  report,  use  dele  of  publication. 

7a.  TOTAL  NUMBER  OF  PACES:  The  total  page  count 
should  follow  normal  pagination  procedures,  i.e„  eater  the 
number  of  pages  containing  information* 

76.  NUMBER  OF  REFERENCES:  Eater  the  total  number  of 
rnfaroace*  cited  ia  the  ropori. 

8a.  CONTRACT  OR  GRANT  NUMBER:  If  appropriate,  eater 
the  applicable  number  of  the  contract  or  grant  under  which 
the  report  was  written. 

86.  Sc,  8  8d.  PROJECT  NUMBER:  Enter  the  appropriate 
military  department  identification,  aaeh  as  project  number, 
aubf-Uct  number,  system  numbers,  teak  number,  etc. 

9o.  OMCNATOR’S  REPORT  NUMBER(S):  Enter  the  oifi- 
cial  report  number  by  which  the  document  will  be  identified 
and  controlled  by  the  originating  activity.  This  number  mart 
be  aaiqae  to  this  report. 

96.  OTHER  REPORT  NUM3ER1S):  II  the  report  has  been 
assigned  any  other  report  numbers  (either  by  the  originator 
or  by  the  epomorX  also  eater  this  numbers). 


10.  AVAIL ABIUTY/LEMTATION  NOTICES:  Enter  any  limi¬ 
tations  on  farther  dissemination  of  the  report,  ether  than  timer 
imposed  by  security  classification,  using  standard  statements 
sack  a a: 

(1)  "Qualified  requesters  may  obtain  copiae  of  this 
report  from  DDC.” 

(2)  "Foreign  announcement  end  diseemi notion  of  this 
report  by  DDC  is  not  authorised." 

(3)  ‘*U.  S.  Government  agencies  may  obtain  copies  of 
■hie  report  directly  from  DDC.  Other  qoalified  DDC 
users  shall  request  through 


(4)  ’  U.  S.  military  agencies  may  oblnia  copies  of  this 
report  directly  from  DDC.  Other  qualified  users 
shell  request  through 


(S)  "All  distribution  of  this  report  in  controlled.  Quail* 
find  DDC  users  shall  request  through 


if  the  report  ha*  been^furniahed  to  the  Office  of  Technical 
Service#,  Department  of  Commerce,  for  eels  to  the  public,  indi¬ 
cate  this  fact  tad  enter  the  price,  if  known. 

11.  SUPPLEMENTARY  NOTES:  Use  for  additional  explana¬ 
tory  notes. 

12.  SPONSORING  MILITARY  ACTIVITY:  Enter  the  inure  of 
the  departmental  project  office  or  laboratory  sponsoring  (pay- 
tag  for )  the  research  and  development.  Include  address. 

13.  ABSTRACT:  Enter  an  abstract  giving  a  brief  and  factual 
summary  ol  the  document  indicative  or  the  report,  even 
though  It  may  also  upper-  elsewhere  in  the  body  of  the  tech, 
nice:  report.  If  additional  space  ie  required,  a  continuation 
sheet  shall  be  attached. 

It  ie  highly  desirable  that  the  abstract  of  classified  re¬ 
port*  be  unclassified.  Each  paragraph  of  tha  abstract  ahull 
and  with  aa  indication  of  the  military  security  classification 
of  tha  information  in  the  paragraph,  represented  as  (TS).  (SX 
(Cx  or  (UX 

There  is  no  limitation  on  the  length  of  the  abstract.  How¬ 
ever,  the  suggested  length  ia  from  130  to  225  words. 

14  KEY  WORDS:  Key  words  are  technically  meaningful  tenet 
or  nhort  phrases  that  characterise  a  —port  end  may  be  ueed  ns 
index  entries  for  cataloging  the  report.  Key  words  must  be 
selected  so  that  no  security  clnsslf  lent  ion  is  required,  identi¬ 
fiers.  such  ts  equipment  model  designation,  trade  name,  mill-, 
tery  project  code  name,  geographic  location,  may  be  used  as 
key  words  but  will  he  followed  by  an  indication  of  technical 
context.  The  assignment  of  links,  rules,  end  weights  is 
optional. 


Unclassified 

Security  Gexalficntioo 


DD  Form  1473 
Abstract  (cont'd.) 


(curies/m2);  (10)  time  of  arrival;  (11)  time  of  cessation;  (12)  smallest  par¬ 
ticle  deposited;  (13)  largest  particle  deposited;  (14)  mass  deposited  by  par¬ 
ticles  in  the  size  range  SI  to  S2;  (15)  H  +  1  hour  "normalized*  exposure 
rate  resulting  from  particles  in  the  size  range  SI  to  S2;  and  (16)  die  number 
of  cloud  (model)  subdivisions  affecting  each  map  grid  point.  The  user  is 
free  to  specify  any  limiting  coordinates  and  scale  factors  for  the  map  display 
that  will  be  produced  and  can  also  cause  the  resulting  map,  or  maps,  to  be 
recorded  on  magnetic  tape  for  further  processing. 


